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ABSTRACT 


This document contains a detailed description of the computer code for the 
trajectory processor (#TRAJ) of the HP-9825A High Fidelity Relative Motion 
ogram (HFRMP). The #TRAJ processor Is a 12-degrees-of- freedom trajectory 
integrator (6 degrees of freedom for each of two vehicles) which can be used 
to generate digital and graphical data describing the relative motion of the 
Soace Shuttle Orb iter and a free-flying cylindrical payload. The description 
of the processor includes a listing of the code, coding standards and conven- 
tions, detailed flow charts, and discussions of the computational logic. 
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dft^/dt 

p 

Keplerian angular velocity magnitude of payload CG, jr x v|/r 
dw^/dt 
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FOREWORD 


This document contains a detailed description of the trajectory processor 
(#TRAJ) of the High Fidelity Relative Motion Program (HFRMP) Version 03T(1229/ 
10N0V79). The HFRMP is coded in Hewlett-Packard Language (HPL) for the 
HP-9825A Desktop Calculator with memory option 002 (23,228 bytes of read/ 
write memory) and the following peripheral devices and read-only memory (ROM) 
modules: 


Required Peripheral Devices 

HP-9885M Flexible Disk Drive 

HP-9872A Plotter 

HP-9866B Thermal Line Printer 

or 

HP-9871 A Character Impact Printer 

Required ROM Modules 

HP-9821 OA String-Advanced Programming 

HP-9821 1A Matrix Programming 

HP-9821 7A 9885M Flexible Disk Drive 

HP-9821 6A 9872A Plotter-General I/O-Extended I/O 

or 

HP-9821 5A 9872A Plotter-General I/O 

Applicable HPL programming information is contained, along with operating in 
structlons for the calculator and its peripheral devices, In the following 
HP-9825A manuals published by the Hewlett-Packard Calculator Products Divi- 
sion: 
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HP Part No. 


Title of Manual 


09825-90000 

09825-90020 

09825-90021 

09825-90022 

09825-90024 

09825-90026 

09871-90000 

09885-90000 


Operating and Programming 
String Variable Programming 
Advanced Programming 
Matrix Programming 
General I/O Programming 
9872A Plotter Programming 
9871A Printer... Operating Manual 
Disk Programming 


A basic understanding of the contents of the cited manuals is required for a 
full comprehension of this program document. 


1 . INTRODUCTION 


The HP-9825A High Fidelity Relative Motion Program (HFRMP) Is a disk- 
resident flight design software system specifically Intended for the analysis 
of Space Transportation System (STS) proximity operations. The mainline logic 
flow + is illustrated in Figure 1. The code Is divided Into several primary 
modules that are designated "processors" and identified as such (according to 
HFRMP convention) using the pound-sign (#) symbol as the first character In 
their identification strings (i.e., names). Some processors, being too large 
to fit in the available read/write memory area, are broken down into subsidiary 
links which are identified by the use of the percent (%) symbol as the first 
character in their identification strings. Subsidiary links overlay each other 
in a shared memory area, as illustrated in Figure 2. 

Processors (such as #TRAJ) are called into read/write memory from the 
program disk by use of the HPL get command, which erases all variables (numer- 
ical data) that may have been stored in memory. Information transfer between 
processors is accomplished by means of disk-resident data files that are 
stored and retrieved, as necessary, by appropriate commands coded into the 
processors. This includes user-supplied input data, which can be manipulated 
(edited/saved/recalled) by the user with the aid of interactive input proces- 
sors (#DBED and #FPED in Figure 1). 

From the preceding remarks it can be seen that every HFRMP processor rep- 
resents, in a limited sense, a stand-alone program. Therefore it is possible. 


4 * 

i.e., that which is invoked by a regular production run of the program. Not 
shown are a number of special processors that are used for program mainte- 
nance. 
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TEST I : 
TEST 2 : 
TC5T3 : 
TEST H : 


( «HMMP ^ 
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SPECIAL OUTPUT PROCESSOR 


1 'N 
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MORE FLIGHT SEGMENTS? 

OIO A PAYLOAD SRM BURN OCCUR? 

PLUME IMPINGEMENT DAMAGE ANALYSIS REQUESTED? 
START A NEW RUN? 

Figure 1. HFRMP Mainline Logic Flow 
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for Instance, to Interrupt the normal HFRMP execution sequence after Input 
processing has been completed and then later to resume execution by comnandlng 
get "#TRAJ M from the calculator keyboard and pressing the (1jU3) button. By the 
same token. It Is not necessary to execute the HFRMP particle Impact damage 
Integrator (#PIDI ) Immediately after executing ITRAJ. The vehicle ephemerls 
data that #PIDI requires Is stored on the program disk by ITRAJ, where It Is 
preserved Indefinitely until written over by a subsequent execution of ITRAJ. 
Therefore, a user can execute the trajectory processor, shut off the calculator 
and remove his program disk, and then Initiate a particle Impact damage analy- 
sis of the subject trajectory at any later date simply by Inserting his disk, 
commanding get "flPIDI 11 from the keyboard, and pressing the CffUtD button. 

The remainder of this document will be concerned only with the #TRAJ pro- 
cessor, which contains the HFRMP dynamic models, the trajectory/ attitude Inte- 
gration logic, and the basic output computations. Detailed program documen- 
tation for the HFRMP Input processors and the #PIDI output processor are to be 
published under separate cover. (Reference 1 contains operating Instructions 
for the HFRMP, with emphasis on Input operations. Reference 2 contains a def- 
inition of the basic equations used In the #PIDI processor.) 
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1.1 GENERAL DESCRIPTION OF #TRAJ 


The #TRAJ processor Is a 12-degrees-of- freedom trajectory Integrator (6 
degrees of freedom for each of two vehicles) which generates digital and graph- 
ical data to describe the relative motion of the Space Shuttle Orblter and a 
free-flying payload. These data are obtained by differencing the geocentric 
states of the individual vehicles, computed to a numerical precision of 12 
decimal digits, with respect to an oblate earth whose gravitational model 
includes the second harmonic coefficient (Jg). The state of the two-vehicle 
system is computed as a function of time by means of a fourth-order Runge- 
Kutta numerical integration scheme which uses unit quaternions (versors) to 
define both the rotational and the translational states of each vehicle. The 
derivation of the versor-based translational state vector is contained In 
Reference 3, which is reproduced for the most part (with some changes in 
coordinate system definition) in Appendix A and Sections 1.2 and 1.3 of this 
report. 

The payload is modeled geometrically as a cylinder whose length and diam- 
eter are specified by program input. The Orblter and the payload are treated 
as rigid bodies whose individual mass properties (gross weight, moments and 
products of inertia, and center of gravity location) are specified by program 
input and are assumed to remain constant 1 during the HFRMP run. Gravity 
gradient torque is included In the rotational equations of motion for both 
vehicles. At the user's option, aerodynamic torque and drag can also be 


+ An exception is made when integrating the trajectory of an upper stage during 
a major translational maneuver such as a solid rocket motor burn. In such 
a case, the payload gross weight is decremented during the burn to reflect 
the consumption of propellant. 
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Included for either or both vehicles. Aerodynamic forces normal to the rela- 
tive wind vector are Ignored. 

The atmosphere Is assumed to rotate with the earth, and Is modeled as a 
function of geodetic altitude by a curve fit of the 1962 Standard density pro- 
file (Reference 4). The curve fit Is valid down to a minimum altitude of 
approximately 400,000 feet. The density profile can be modified by a program 
Input factor (K p ) to account for the major effects of solar activity. Aerody- 
namic drag and moment coefficients for the Orblter are computed as curve-fit 
functions of Its attitude with respect to the relative wind vector (References 
5-7). 

Aerodynamical ly, the payload Is modeled as a flat plate whose size and 
shape are determined by projecting Its cylindrical outline onto a plane normal 
to the relative wind vector. The payload drag coefficient Is assumed to be 
2.0, based on Its projected frontal area. Aerodynamic effects on the Orblter 
and the payload can be modified (or cancelled entirely) by means of Input fac- 
tors (Ka and ka) which are applied uniformly to all aerodynamic forces and 
torques that are computed for the specified vehicle. The possible effects of 
aerodynamic shadowing (of one vehicle by the other) are not accounted for In 
the Internal calculations of #TRAJ. 

Several options are provided for defining the Initial state of the two- 
vehicle system. The translational state of the Orbiter can be described 
either In terms of osculating orbit elements referenced to the Mean of 1950.0 
(M50) geocentric equatorial frame, or in terms of Invariant orbit elements 
(Reference 8) measured in the Mean of Launch Oate (MLD) equatorial frame. The 
initial attitude of the Orbiter Is defined by pitch, yaw, and roll angles 
(taken in that order) referenced to the rotating Shuttle-centered local- 
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vertical (SLV) coordinate system. The Orblter's angular velocity, measured 
In terms of rate components about Its body axes, can be defined with respect 
to either the M50 (Inertial) frame or the SLV (rotating) frame. 

The Initial translational state of the payload Is defined by rectangular 
position and velocity components which are measured relative to the Orblter’s 
center of gravity (CG). At the user's option, these components can be mea- 
sured In the SLV coordinate system, or In the Shuttle body (SBY) coordinates 
system. The payload's initial pitch, yaw, and roll angles can be referenced 
either to the payload-centered local-vertical (PLV) system, or to the SBY 
system. The payload's angular rate components about Its body axes can be 
defined relative to the M50, the PLV, or the SBY frai;:i. 

The initial state of the system can be advanced through up to 40 flight 

profile segments, each of arbitrary length, which are defined by the user in 

a pric' 1 execution of the #FPED processor/ At the beginning of any segment 

the user may command the application of an impulsive (l.e., instantaneous) 

* 

Increment to the angular rate of either or both vehicles. In this regard, 
the user may specify a particular rate Increment (INCR), a desired rate with 
respect to inertial space (IR), or a desired rate with respect to the local- 
vertical frame of the vehicle In question (LVR). In all cases, the components 
of the desired rate or rate increment are measured about the body axes of the 


All input data, including the flight profile definition, are saved in disk 
files whence they can be recalled (and edited. If necessary) for use In 
subsequent runs. 

Commanded angular velocity impulses, and the linear velocity impulses which 
they may induce as a result of RCS translational cross coupling, are the only 
types of state variable discontinuity that are permitted by #TRAJ. These are 
allowed only at the beginning of a flight profile segment. 
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vehicle In question. 

After applying the specified angular velocity Increment (if any) at the 
beginning of the flight segment, #TRAJ then (for each vehicle independently 
of the other) 

1. allows the attitude to drift (D) under the influence of Inertia and 
natural torques, or 

2. performs inertial rate hold (IRH) control (i.e., maintains a constant 
angular velocity relative to inertial space), or 

3. performs local -vertical rate hold (LVRH) control (i.e., maintains a 
constant angular velocity relative to the rotating local -vertical 
frame of the vehicle being controlled) 

for the duration of the segment, depending on the user's specifications for 
that segment. 

When the IRH or the LVRH attitude-maintenance option is specified for the 
Orbiter, a simplified RCS/DAP model (Reference 9) is used to compute average 
values for the propellant consumption rates and translational cross-coupling 
accelerations that result from the intermittent thruster firings which are 
required to apply the necessary control torques. The model takes into account 
the mass properties of the Orbiter, the electrical width (an integer multiple 
of the DAP cycle time) anc the effective width (the duration of steady-state 
acceleration) of the RC3 thruster pulses, and the width of the attitude dead- 
band about each of the Orb iter's body axes. Deadbands can be changed from 
segment to segment in the flight profile, as can the selection of primary or 
vernier thrusters and the mode of cross-coupling compensation. Translational 
cross-coupling accelerations are integrated along with those produced by 
gravity, aerodynami cs , and commended translational thrust. They are reflected 
in the output data by the flight path deviations they produce. Propellant 


consumption rates are also integrated (but not subtracted from the Orblter 
gross weight), and the accumulated expenditures are tabulated, along with other 
data, at user-specified time intervals. Pronellant consumption is broken down 
according to source (forward, aft left, or aft right tank) and function (trans- 
lational or rotational control). 

When the IRH or the LVRH attitude-maintenance option is specified for the 
payload, the magnitudes of the necessary control torques are integrated and the 
accumulated rotational impulse (measured in pound-foot-seconds) in the positive 
and negative direction about each body axis is printed along with the other 
digital output data. Since no specific method of implementation is modeled, 
it is not possible to compute propellant consumption rates or cross-coupling 
effects that may result from payload attitude control. 

Translational thrust acceleration of either or both vehicles can be 
commanded at the beginning of any flight profile segment. Payload transla- 
tional thrust is always applied in the direction of the payload's +X body axis 
and is assumed to be directed through the CG. Once initiated, payload thrust 
acceleration continues until all of the rocket motor propellant is consumed, 
as determined by a table of flow rates versus burn time. 

Translational acceleration of the Orbiter is initiated by commanding 
ignition of either or both of the OMS engines (L, R, or L+R) and/or by firing 
primary RCS thrusters to produce thrust nominally in the positive or negative 
directions of the Orbiter body axes (+X, -X, +Y, -Y, +Z, or -Z). Once initi- 
ated, Orbiter translational acceleration is applied continuously at the 
nominal steady-state level, throughout the duration of the flight profile 
segment. 
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1.2 STATE VECTOR DEFINITIONS 


There are many numerical methods which can be used to integrate the dif- 
ferential equations that govern the motion of a rigid body such as the Orbiter 
or a free-flying payload. A common requirement of such methods is the defi- 
nition of a state vector + along with a set of equations which, given the val- 
ues of the state variables at any particular time, permit the calculation of 
the first derivative of every state variable with respect to time. 

A fourth-order Runge-Kutta numerical integration technique is used in the 
#TRAJ processor. However, the state vector definitions and the derivative 
calculations (which comprise the bulk of the computations involved in state 
vector propagation) are not influenced directly by the integration logic 
except in the sense that the integrator determines how frequently the deriv- 
atives have to be evaluated. In #TRAJ, the integration logic and the deriva- 
tive calculations are isolated from each other in separate subroutines, thus 
making it easy to change the numerical integration technique if that should 
become desirable. 

The purpose of this section of the report is to describe the state vari- 
able that have been adopted for use in #TRAJ. The associated derivative 
calculations, described in terms of external forces and torques that act on 
the body under consideration, will be covered in Section 1.3. The detailed 
calculation of these forces and torques will be described in Section 2.3, and 
the Runge-Kutta integration logic will be described in Section 7.2. 


4 * 

A set of variables that define the position and linear velocity of the body's 
CG (in the case of translational motion) with respect to a chosen system of 
reference, and/or the body's attitude and angular velocity (in the case of 
rotational motion about the CG). 
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The #TRAJ processor makes use of versors (unit quaternions) as Internal 
state variables to describe the translational as well as the rotational 
motions of the Orblter and the payload. (The HFRMP user never sees the state 
vectors in quaternion form; input and output data are defined in terms of 
Euler angles.) The computational advantages of using versors to describe 
rotational motion are well known. Certain advantages also accrue from their 
use to describe the translational motion of a satellite, as explained in 
Reference 3. Versors are also used extensively in #TRAJ for coordinate trans- 
formations in general, and a working knowledge of their characteristics and 
rules of manipulation is necessary for a good understanding cf the remainder 
of this report. A careful _tudy of Appendix A is recommended at this point 
for the reader who is not accustomed to the use of quaternions. Even for 
those who are familiar with quaternion applications, it probably will be help- 
ful to review the conventions and the system of notation defined in Appendix 
A. All readers should examine the coordinate system definitions in Appendix B. 

1.2.1 Rotational State Vector 

The rotational state of the Shuttle at any given time t is defined by the 
—IB =; 

variables [q jg ,a) B where Qjb is a versor defines the orientation of the 

Shuttle body-fixed (SBY) coordinate system B with respect to the geocentric 
mean-of-launch-date (MLD) equatorial inertial coordinate system I. The vector 

T n 

symbol oi B represents the angular velocity of the Orb iter. The double super- 
script IB identifies the angular velocity as that of system B with respect to 
system I, and the subscript B indicates that the vector is resolved into 
components along the axes of the B system. The rotational state of the payload 
is similarly defined by the variables whose definitions are the 

same in all respects except that the lower-case symbol b is used to identify 
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the body-fixed coordinate system of the payload. 

1.2.2 Translational State Vector 

The translational state of the Orbiter is defined by the variables 
[qjQ,R,R,fi K ]. The symbol qjQ represents a versor that defines the orientation 
of the Shuttle local -vertical (SLV) coordinate system G with respect to the 
inertial system I. The symbol R represents the distance from the center of 
the earth to the Shuttle CG, and R stands for its derivative dR/dt. The 
symbol represents the scalar magnitude of the instantaneous geocentric 
angular velocity of the Shuttle CG. The subscript K in this case does not 
designate a coordinate system; it merely identifies the angular velocity in 
question as that which is associated with Kepler's law of equal areas. 

The translational state of the payload is similarly composed of the vari- 
ables [q Ig ,r,r,w K ], which are defined in the same manner as those described in 
the preceding paragraph except for the use of lower-case symbols to repre- 
sent the scalar quantities and the coordinate-system designator (g) of the 
payload local -vertical (PLV) coordinate system. 
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1.3 BASIC EQUATIONS OF MOTION 


The purpose of this section is to define the equations for calculating 
the derivatives of the #TRAJ internal state variables, assuming the forces and 
torques acting on the bodies are known. Equations for calculating the forces 
and torques will be given in Section 2.3. 

The only equations that will be defined explicitly in this section are 
those that apply to the motion of the Shuttle. The payload equations are 
identical in all important respects, and can be obtained simply by substituting 
the lower-case symbols b, g, r, r, and for their upper-case counterparts 
B, G, R, R, and 

1.3.1 Rotational Motion 

^ TD 

What is required here are the first derivatives of q JB and u> B with re- 
spect to time. The derivative of the orientation versor is equal to one-half 

T n 

of its quaternion product with , i.e., 



TD 

The derivative of w B is given by the equation 

W B = ( L b " “B x CID b W B )’ (2) 

where the vector Lg represents the external torque acting on the Shuttle. 

The symbol [I]g represents the Shuttle's inertia tensor (a 3x3 matrix whose 
elements consist of the moments and negative products of inertia, referenced 
to the B system), and [I]^ 1 represents its inverse. For the purpose of comput- 
ing the product of a matrix and a vector, the vector is treated as a column 
matrix, e.g. , 
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mama atom 


-IB . 
W B " 


1.3.2 Translational Motion 

We seek now to define the first derivatives of the translational state 
variables q"j G , ft, and with respect to time. (The derivative of the state 
variable R is no problem, because it is equal to R). By definition, the 
Shuttle's local -vertical coordinate system 6 rotates in such a manner that 

R G = - k R (3) 

and 

(RxV) g = - j|RxV| = - j R 2 n K (4) 

at every instant of time, where R and V represent the geocentric position and 
inertial velocity vectors of the Shuttle CG. 

An expression for the inertial velocity can be obtained by differentiat- 
ing Equation (3) with respect to time. This yields 

V G = - k R - x k R, (5) 

where the vector 

— IG _ 2 IG , • IG . l IG ic\ 

W G " 1 “Gl + ^ W G2 + ^ W G3 6 

(whose components have yet to be evaluated) represents the total angular 
velocity of coordinate system G with respect to the inertial system I. We 


IB 

W B1 

IB 

W B2 

IB 

“B3 
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.V 1 '•** ’ ^ 


now differentiate Equation (5) with respect to time and obtain the expression 
/T 6 = - k R - 2 x k R 

- wq G x k R - wg G x (“q G X k r) (7) 

for the linear acceleration of the Shuttle CG, where the vector 


-4G 

“G 


IG 


IG 


_ „• * VJ , • , | 

- » “61 + J “C9 + k 


G2 


•IG 

°G3 


( 8 ) 


represents the angular acceleration of coordinate system G. 

T p J_T p 

We now proceed to evaluate the components of tog and . First, we 

Tp 

substitute the expression for wg from Equation (6) into (5), and expand the 
vector cross product. This results in the equation 


V- 


i R 


IG 

W G2 


+ j R 


«4 G - k R. 


(9) 


Using the expression tor Rg and Vg from Equations (3) and (9), we form the 
cross product 

(RxV)g = 1 R 2 a,g G + j R 2 ug G . (10) 

Equating the vector components of Equations (4) and (10), we obtain 

<4? = 0 (ID 

and 

W G2 ‘ ’ »K <' 2 > 


which lead to 


\ 


and 


•16 _ • 
“G2 “ ’ 


( 14 ) 


Substituting (11), (12), (13), and (14) Into Equations (6), (8), and (9), we 
obtain 


-16 . : . . r 16 

' W G - 0 + k 


-IG . ;• r -ig 

W G " ' J + k “G3* 


(15) 

(16) 


and 


V G * 1 ' R n K ‘ k R* 


(17) 


T P * T p 

Using the expressions for wg and wg from Equations (15) and (16), the 
cross products appearing in Equation (7) can now be expanded to obtain 

= i (2 R + R a K ) + j R Wg^ + k (R - R). (18) 

From Newton's second law, we have 


Fg = M Ag, 


(19) 


where M represents the mass of the Shuttle and the vector 


Ff: = ^ + J + K 


G2 G3 


( 20 ) 


represents the external force acting on it. Combining Equations (18), (19), 
and (20), and equating vector components, we obtain the equations 


R - R uj - F g3 /M 


( 21 ) 
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h K « (F 61 /M - 2 R n K )/R 


( 22 ) 


and 

“G3 * *W (M R V- (23) 

Equations (21) and (22) define the derivatives of the translational state 
variables R and The derivative of the remaining state variable qj G Is 
given by 

q IG - (q IG » « g )/2. (24) 

where the angular velocity vftctor 

4 G = - j + k F G2 /(M R n K ) (25) 

is obtained by substituting Equation (23) into (15). 

Equations (21), (22), (24), and (25) are valid for all states of motion 
except those where |Rx?| = 0, In which case the geocentric flight path would 
be vertical, and the orientation of the X and Y axes in the local -vertical 
coordinate system would be Indeterminate. 
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1.4 CONVENTIONS AND CODING STANDARDS 


HP-9825 programs are coded In the Hewlett-Packard Language (HPL), which 
Is functionally similar to BASIC but unique In terms of Its syntax. To 
facilitate correlation with the #TRAJ code (Appendix C), some of the syntacti- 
cal features of HPL have been carried over Into the flow charts (Sections 2-7). 

For instance, consider a computation which Involves forming the product 
of two simple variables A and B, where the resulting value Is to be assigned 
to a third variable C. In FORTRAN code this operation would be represented by 
writing C » A ♦ B , and in most flow charts by writing C « AB . In HPL code, 
and in the flow charts that follow, the same operation would be described by 
writing AB •» C . (Perhaps the most succinct verbalization of the HPL value- 
assignment symbol is obtained by substituting the phrase "goes to" for the 
right- running arrow 

Another feature of HPL code and the flow charts, which may be confusing 
to FORTRAN programmers, is the use of square brackets to enclose the index or 
indices of a dimensioned variable. The Jth element of the HPL dimensioned 
variable A[*] Is designated by writing A[J], whereas In FORTRAN it would be 
written A(J), In HPL code, the square brackets are always used to enclose 
indexing or dimensioning information about an array (dimensioned) variable, 
and for no other purpose . This Is necessary in HPL to distinguish, for 
instance, the simple (scalar) variable P from the array variable P[*] whose 
dimensions might have been declared by a statement dim P[5,10] , which Is 
analagous to the FORTRAN statement DIMENSION P(5,10) . In a FORTRAN program 
where such a dimension statement exists, a reference to the variable P is an 
implicit reference to P(l,l), l.e., the first element in the P array. However, 
this is not the case in HPL, where the simple variable P has no connection 
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whatever with the array variable P[*] or any of Its elements. For Instance, 
the sequence of HPL statements 2 ■» P; 3 -» PCI ,1]; P PC1.1] •» X would cause the 
number 6 to be assigned to the variable X. 

1.4.1 General Arrangement of Program Files 

The first line of code (line 0) In every HFRMP program file (processor or 
processor link) contains, In the form of a label, the name under which the file 
Is stored on the program disk, followed by the date and time of Its most recent 
revision. The executable part of this line contains the statement gto "RUN" , 
where "RUN 11 Is the label of the line where the main logic flow begins. The 
next t.'n lines (1 and 2) contain statements that facilitate an automatic list- 
ing of the complete program, and which are never executed in a normal produc- 
tion run. 

All subprograms are located between line 2 and the "RUN" label, and they 
are arranged in the following general order: 

1. Function subprograms 

2. Utility subroutines 

3. Special subroutines. 

The rationale for this order is to make the program listings as Insensitive as 
possible to corrections or revisions of the logic. Function subprograms and 
utility subroutines are revised much less frequently than the main (driver) 
logic in a given link. Therefore the code In the top part of each link tends 
to remain comparatively stable. 

To the greatest extent that Is practicable, the code In each subroutine 
and in the driver logic of each Individual link Is arranged so that execution 
control flows from top to bottom, i.e., from smaller line numbers to greater 


ones. 




1.4.2 Branching Conventions 


For ease of program maintenance, branching to absolute line numbers Is 
avoided. Relative addressing (e.g,, gto + 4) Is used whenever practicable; 
otherwise, the general rule Is to branch to a symbolic label (e.g., gto "RUN") . 
The jmp statement is used only In the rare Instances where It Is necessary to 
compute a relative address (line number) at execution time. 

To avoid confusion regarding the flow of execution control (and also 
generally to permit the transfer of date through an argument list), the 
HPl cfl statement is used for branching to a subroutine In lieu of the gsb 

sta tenant. 

1.4.3 Program Labels 

No lower-case alphabetic symbols are used In program labels (character 
strings used to identify branching targets, Including the names of subpro* 
grams). The purpose of this convention Is to provide a contrast between such 
labels and the commands and standard functions of HPL, all of which are spelled 
with lower-case alphabetic symbols. 

1 .4.4 Input/Output Mnemonics 

Although numerical option codes are sometimes used internally by #TRAJ, 
it has been adopted as a policy that HFRMP users should not have to memorize 
(or look up) obscure numerical codes in order specify logical options In the 
Input data files. Rather, such options are specified by the user in the form 
of mnemonic character strings which, if necessary, are converted to numerical 
codes Internally. 

Mnemonic character strings are also used to identify the HFRMP digital 
output data. With regard to all input/output mnemonics, a determined effort 
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has been put forth to make them as meaningful and as consistent as possible. 


1.4.5 Data-Reqlster Protocol 

Because of the limitations of the HPL syntax, It has been 'ound necessary 
to adopt a comparatively rigid protocol to govern the use of named variables 
and of the numbered r-reglsters. Appendix D contains storage allocation Infor- 
mation for the r-reglsters. 

1.4. 5.1 Simple Variables 

The simple (scalar) variable names, which are limited In HPL to single 
upper-case alphabetic characters (A thru Z), are used primarily as loop 
counters and for the temporary storage of Intermediate computational results. 
Eight of the 26 simple variables are designated as volatile registers, which 
means that they are least rigidly controlled, and that their contents are most 
susceptible to frequent change. These are the registers Identified by the 
characters H,I,J,K and W,X,Y,Z. 

1.4. 5. 2 r-Reglsters 

Because of the limited number of usable names for variables, and also 
because array names cannot be passed through the argument lists HPL sub- 
programs, most of the data that would normally be assigned unique names or 
stored In Individual arrays (e.g., as In a FORTRAN program) are stored Instead 
In the numbered r-reglsters. The correlation between the r-numbers and the 
logical symbols, shown In Appendix D, Is extremely critical to the understand- 
ing of the #TRAJ code. 

When there Is a need for a utility subprogram to perform a standard compu- 
tation involving one or more logical arrays (such as vectors and quaternions). 
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the address of each logical array (i.e., the r-number of its first element) is 
passed through the argument list (in lieu of an array name, which is not per- 
mitted in HPL). For instance, suppose the components of a vector Wg reside 
in the registers rl through r3, the components of another vector Vg reside In 
r4 through r6, and it is desired to compute the vector ffg « Rg x Vg and store 
it in r7 through j Symbolically, in the flow charts, this would be written 
Rg x Vg -> Hg, or possibly as c// 'CRSP' (Rg,Vg,Hg). In the HPL code, it would 
appear as cl/ 'CRSP 1 (1 ,4,7) , where 'CRSP' is the name of the vector cross- 
product routine, and the numbers 1,4,7 are the addresses of the vectors in- 
volved in the operation. 

The first 19 r-registers (rO through rl8) are similar to the simple-vari- 
able registers in that they are used mainly to store the intermediate results 
of array computations. The first ten r-registers (rO through r9) are analo- 
gous to the simple variables H,I,J,K and W,X,Y,Z in that they are designated 
as volatile, i.e., most frequently used. Preferential use of the lower-num- 
bered registers tends to minimize code-storage requirements, simply because 
they entail writing out fewer digits (each of which requires one byte of code 
storage) tc: identify the registers. That is to say, the statement c // 'CRSP ' 
(101 ,104,107) requires 6 more bytes of code storage space than the statement 
c // 'CRSP' (1 ,4,7) . 

1.4. 5. 3 Restrictions on Utility Subprograms 

During the execution of an HPL subprogram, the HP-9825 operating system 
allocates temporary storage (only for the duration of subprogram execution) to 
numbered p-variables. These variables are numbered sequentially beginning 
with pO, into which the operating system loads the number of parameters that 
appear in the argument list of the calling routine's reference to the sub- 
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program. If there should be three arguments, then at the beginning of sub- 
program execution pO would contain the number 3. The registers pi, p2, and p3 
would contain the values assigned to the arguments of the subprogram by the 
calling routine. The subprogram can make use of as many higher-numbered p-reg- 
isters (p4, p5, etc.) as it may need for temporary storage of intermediate 
computational results. All of the p-registers are de-allocated (in effect, 
erased) when execution control is returned to the calling routine. 

The operating system also permits the subprogram to reference (get values 
from or store values into) any of the registers accessible to the main program 
(this includes the simple variables, array variables, string variables, and 
r-registers). However, to prevent inadvertent modification of the contents of 
registers that the calling routine may be using, certain conventions have been 
adopted that limit the access of utility subprograms to global variables (i.e., 
any other than the p-variables). 

In general, utility (general-purpose) subprograms are not allowed to make 
literal reference to &ny r-register. They are allowed to make logical refer- 
ence to such registers by means of addresses that may be passed to them 
through their argument lists by the calling routine. For example, the charac- 
ter strings rO or r25 are not permitted in the code of a utility subprogram. 
However, the character strings r£l_ and r(p3+3) are permitted, where pi and p3 
represent addresses (in this case, the first and third numbers in the argument 
list) of logical arrays that are passed to the subprogram by the calling 
routine. 

Function (as distinguished from subroutine) subprograms are not allowed to 
make reference to any simple variable. They are required to use p-variables 
for any temporary storage they might need. The same restriction applies to 


any utility subroutine whose logical argument list contains only scalars and 
vectors . 

Utility subroutines that perform quaternion and matrix operations are 
permitted to use the volatile simple variables (H,I,J,K and W,X,Y,Z), but no 
others. Special subroutines such as ' DERVIS ' (Section 2.3), whose calculations 
generally are more complicated than those of utility subroutines, have no 
such general restrictions on their access to storage registers. 


2. BASE LINK (#TRAJ) 


The "#TRAJ" program file is the base link of the #TRAJ processor, contain- 
ing all of the general-purpose subprograms, i.e., all that are used by two or 
more subsidiary links. 

2.1 FUNCTION SUBPROGRAMS 

2.1.1 * ANG1 1 

The function subprogram 'ANG1' converts an input angle to its equivalent 
value in the range of 0 to 2ir. 

2. 1.1.1 Arguments 

pi = Input angle, measured in radians. 

2. 1.1. 2 Example of Usage 

The instruction 'ANGl 1 (-3 tt) •» A in the calling routine would cause the 
value of +ir to be assigned to the variable A. 

2. 1.1. 3 Computations 

The fractional part of pl/2ir is multiplied by 2ir. If the result is nega- 
tive, 2tt is added. A flow chart is shown in Figure 3. 
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The function subprogram ' ANG2 ' converts an input angle to its equivalent 
value in the range of -ir to +u. 

2. 1.2.1 Arguments 

pi = Input angle, measured in radians. 

2. 1.2. 2 Example of Usage 

The instruction 1 ANG2 1 -» A in the calling routine would cause the 

value of -ir to be assigned to the variable A. 

2. 1.2. 3 Computations 

The fractional part of pl/2ir is multiplied by 2ir. If the result is great 
er than ir, 2-n is subtracted. A flow chart is shown in Figure 4. 







2.1.3 'ATNl '/'ATN2 1 


This multiple entrypoint function subprogram computes the arctangent of 
the ratio of two input quantities, whose individual signs determine the quad- 
rant of the output angle. The ' ATN2 ' entrypoint is analogous to the FORTRAN 
ATAN2 function, in that the output angle lies in the range of -n to +*. When 
the 'ATNl' entrypoint is used, the output angle lies in the range of 0 to 2*. 

2.1.3. 1 Examples of Usage 

The instruction 'ATNl 1 (-1 ,1) •» A would cause the variable A to be as- 
signed the value of 7 tt/ 4 , while the instruction 'ATN2'(-1 ,1) -» A would cause 
it to be assigned the value of -u/4. 

2. 1.3. 2 Computations 

Tests are made to avoid division by zero. Depending on whether pi or p2 
has the greater magnitude, a base value of the output angle is computed from 
the expression tt/2 - atn(p2/abs(pl)) or the expression atn(abs(pl/p2)) , where 
atn and abs_ represent the HPL arctangent and absolute value functions. Tests 
on the signs of pi and p2 then are made so that the base value can be rotated 
into the proper quadrant. If pi * p2 * 0, the output angle is set equal to 
zero. A flow chart is shown in Figure 5. 
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2.1.4 'HM.S' 

The 'HM.S' function converts time from seconds to hours* minutes, and 
seconds. The output Is a single nwnber In the form of HHM.SS, where HH 
denotes the hours digits, MM denotes the minutes digits, and SS denotes the 
seconds digits. 

2. 1.4.1 Arguments 

pi * Input time, measured In seconds. 

2. 1.4. 2 Example of Usage 

The Instruction 'HM. $'(36385. 8) -*• H In the calling routine would cause 
the variable H to be assigned the value 1006.258 (10 hours, 6 minutes, and 
25.8 seconds). 

2. 1.4. 3 Computations 
See Figure 6. 
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2.1.5 * SECS * 

The 'SECS' function Is the Inverse of the 'HM.S' function described In 
Section 2.1.4. It converts time from hours, minutes, and seconds (expressed 
In the HHMM.SS format) Into seconds. 

2. 1.5.1 Arguments 

pi » Input time, In the HHMM.SS format. 

2. 1.5. 2 Example of Usage 

The Instruction 'SECS' (1006.258) •» S In the calling routine would cause 
the variable S to be assigned the value 36385.8 (the number of seconds In 10 
hours, 6 minutes, and 25.8 seconds). 

2. 1.5. 3 Computations 
See Figure 7. 
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2.1.6 *001? 1 

The 'DOTP* function subprogram computes the dot (scalar) product of two 
vectors. 

2. 1.6.1 Arguments 

pi * Address of first vector. 
p2 ■ Address of second vector. 

2. 1.6. 2 Example of Usage 

If vector F resides In registers rl, r2, r3, and vector F resides In 
registers r4, r5, r6, then the Instruction 'DOTP' 0.4) -» D In the calling 
routine would cause the value of F • F to be assigned to the variable 0. 

2. 1.6. 3 Computations 
See Figure 8. 
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2.2 UTILITY SUBROUTINES 


2.2.1 'SXV* 

The subroutine 'SXV' calculates the product of a scalar with a vector. 

2. 2. 1.1 Arguments 

pi = Value of scalar. 

p2 = Address of input vector. 

p3 = Address of output vector. 

2. 2. 1.2 Examples of Usage 

The instruction c/i 'SXV' (6. 8,0,3) In the calling routine would cause the 
values 6.8 rO, 6.8 rl , 6.8 r2 to be stored In r3, r4, r5. If vector IT resides 
In rll, rl2, rl3, and IT resides In r7, r8, r9, the Instruction c/l 'SXV'CDOTP 1 
(11, 7), 11 >7) would cause the value of (A • B) J to be stored In r7, r8, r9. 

2. 2. 1.3 Computations 
See Figure 9. 
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Figure 9. 'SXV Subroutine Logic Flow 



Depending on which entrypolnt Is used, this subroutine will calculate 
either the sum or the difference of two vectors. 

2. 2. 2.1 Arguments 

pi ■ Address of first Input vector. 
p2 * Address of second Input vector. 
p3 = Address of output vector. 

2. 2. 2. 2 Example of Usage 

If vector A resides in rl, r2, r3 and vector IT resides In r4, r5, r6, the 
Instruction cj£ 'VADD*(1 ,4,7) in the calling routine would cause the value of 
A" + IT to be stored In r7, r8, r9. The Instruction c/i 'VSUB'O ,4,4) would 
cause the value of A - F to be stored in r4, r5, r6. The Instruction 
ell 'VADD' (4,4,4) would cause the value of 2F to replace F In r4, r5, r6. 

The Instruction c U 'VSUB'O ,1 ,1) would cause zeros to be stored In rl, r2, r3. 

2. 2. 2.3 Computations 


See Figure 10. 








2.2.3 * CRSP 1 


The 'CRSP* subroutine computes the cross product of two Input vectors. 

2.2.3. 1 Arguments 

pi * Address of first Input vector. 
p2 * Address of second Input vector. 
p3 = Address of output vector. 

2. 2. 3. 2 Example of Usage 

If the vector J resides in rl, r2, r3 and the vector IT resides In r4, r5, 
r6, the Instruction c // 'CRSP'O .4.7) in the calling routine would cause the 
value of X x F to be stored in r7, r8, r9. The Instruction c 1/ 'CRSP 1 (4,1 ,4) 
would cause the value of IT x ft to replace IT in r4, r5, r6. 

2. 2. 3. 3 Computations 
See Figure 11. 
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Figure 11. • CRSP 1 Subroutine Logic Flow 



2.2.4 *ROT '/ * IROT* 


This multiple entrypolnt subroutine Is used for transforming the com- 
ponents of a vector between two Cartesian coordinate systems that are angularly 
displaced from one another by a rotation about their common origin. Let the 
two coordinate systems be designated F and G, and let the versor (unit quater- 
nion) q F6 define the orientation of G with respect to F (see Appendix A). 

The 'ROT' entrypolnt performs the computation symbolized by ^p G o Vp o q’pg Vq. 
The symbol Vp represents the value of a vector 7 In coordinate system F, and V G 
represents the value of the same vector In system G. The IROT entrypolnt 
performs the Inverse transformation, symbolized by q FG o V fi o ?FG ’ V 

2. 2. 4.1 Arguments 

pi = Address of input vector. 

p2 * Address of versor. 

p3 * Address of output vector. 

2. 2. 4. 2 Examples of Usage 

If the versor c^q resides in rO, rl , r2, r3 and the vector Vp resides 
In r4, r5, r6, the Instruction cl/ 'ROT '(4,0,7) In the calling routine would 
cause the vector V G to be stored in r7, r8, r9. A subsequent Instruction 
ell * IROT ' (7,0,7) would cause a copy of the vector Vp to replace Vq In r7, r8, 
r9. The Instruction ell 'ROT' (4,0,1) would cause V G to be stored In rl, r2, 
r3 (thus destroying the last 3 components of the versor q^). 

2. 2. 4. 3 Computations 

The quaternion operation q FG o o q pG -*• ? G Is equivalent to the matrix 
operation [R] Vp -*■ V G , where 
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(HH+II-JJ-KK) 

2(IJ-HK) 

2(IK+HJ) 


2(IJ+HK) 

(HH-II+JJ-KK) 

2(JK-HI) 


2(IK-HJ) 

2( JK+HI ) 
(HH-II-JJ+KIO 


and where q FG ■ H + 1 I + j J + k K. The Inverse transformation 
q F G o o cfp G •*> Vp Is equivalent to the matrix operation [R] T Vg ■ 
[R] T Is the transpose of [R]. By Inspecting Equation (26), It can 
[R]^ Is obtained simply by reversing the sign of H, as reflected 1i 


( 26 ) 


■ Vp, where 
be seen that 
i the flow 


chart that Is shown In Figure 12. 


tor'll, 


rf>i4x 
r(^i4i>» v 

r(j>i4aV»S 

r(|»+0*i 
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Figure 12. * ROT V I ROT ’ Subroutine Logic Flow 






2.2.5 *QDOT */ *QXQV *QXQC */ *QCXQ 1 


This multiple entrypolnt subroutine performs a variety of related quater- 
nion operations. As In Section 2.2.4, let qpg represent a versor (unit 
quaternion) which defines the orientation of a Cartesian coordinate system 
6 with respect to another system F that has the same origin. Let the vector 
Mg’ represent the angular velocity (measured In system G) of system G with 
respect to system F. The 'QDOT 1 entrypolnt Is used to perform the operation 
symbolized by (q F g o wg )/2 q F g, where q F6 represents the derivative of q pG 
with respect to time. 

Let F and F represent any two quaternions (not necessarily versors), and 
let F and 5 represent their conjugates. The entrypolnts 'QXQ', 'QXQC', and 
' QCXQ * (respectively) perform the operations symbolized by F o $ •* F, 

F o ? + F, and F o F F. 

2.2. 5.1 Arguments 
'QDOT* Entrypolnt : 

pi * Address of the orientation versor qj-g. 

rp 

p2 ■ Address of the angular velocity vector JZ . 
p3 * Address of the output derivative qpg. 

'QXQ 7 'QXQC 7 'QCXQ' Entrypolnts : 

pi * Address of the first Input quaternion F. 
p2 > Address of the second Input quaternion F. 
p3 * Address of the output quaternion F. 

2. 2. 5. 2 Examples of Usage 

Assume the variable A has the value of 25, that the versor q^ resides In 

rp 

r50, r51, r52, r53, and that the vector uig resides In r54, r55, r56. Then 
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a 

the statement c// l 0D0T l (A+25.A+29.A+125) will cause the value of qpg to be 

stored In r150, rl51, rl52, r!53. If the versor q^ resides In rO, rl, r2, r3, 

the statement c U , QXQ'(A+25.0.4) will cause q fH to be stored In r4, r5, r6, 

r7. If this were followed by c/l 'QXQC' (0.4.8) , the versor 

qgH o » q 6H o q H p » q G p would be stored In r8, r9, rlO, rll. Then the 

statement c/l 'QCXQ' (8.0.0) would cause a copy of the versor 

%F 0 ?6H * ^FG 0 ^GH * ^FH ( alrea< *y residing In r4, r5, r6, r7) to be stored 

In rO, rl, r2, r3. 

2. 2. 5. 3 Computations 

The computational sequence depicted by the flow chart In Figure 13 results 
from a straightforward application of the rules of quaternion algebra that are 
defined In Appendix A. 
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Figure 13. 'QDOT'/W/'QXQC'/'QCXQ' Subroutine Logic Flow 
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2.2.6 ‘Q2317'Q213' 


The digits 1,2,3 In the *Q231 ' and 'Q213' entrypolnt names Identify coor- 
dinate-system rotations about the X,Y,Z axes (respectively) In two different 
Euler-angle sequences. Generally, rotation about the X axis Is designated as 
roll, rotation about the Y axis is designated as pitch, and rotation about the 
Z axis Is designated as yaw. Given an orientation versor qpg, the 'Q231' entry- 
point- Is used to calculate the values of the pitch, yaw, roll angles (taken In 
that sequence) that are required to rotate the axes of coordinate system F Into 
coincidence with those of system G. The 'Q213' entrypolnt Is used to calculate 
the corresponding Euler angle set for a pitch, roll, yaw sequence. 

2.2.6. 1 Arguments 

pi • Address of the orientation versor. 

p2 • Address of the output array of Euler angles (measured In radians, and 
stored In the same order as the logical rotation sequence). 

2. 2.6. 2 Example of Usage 

If the orientation versor under consideration resides In rO, rl, r2, r3, 
the statement cil 1 Q231 (0,4) would cause the first rotation (pitch) angle to 
be stored In r4, the second rotation (yaw) angle to be stored In r5, and the 
third rotation (roll) angle to be stored In r6. The statement c// *0213(0,0) 
would cause the rotation angles of a pitch, roll, yaw sequence to be stored In 
rO, rl, r2 (thus wiping out the first three conponents of the orientation 
versor) . 


2. 2. 6. 3 Computations 

The output of the 'Q2317'Q213' subroutine Is an ordered set of Euler angles 


49 


a, 6, y. For a pitch, yaw, roll sequence the Euler angles are related to the 

- AAA 

orientation versor q^ s H + i I + j J + k K by the equations (see Reference 

10) 


H a C C 0 C - S S c S , 
a 3 Y o 3 y 

(27) 

I ■ c c 0 s + s s 0 c , 

a 6 Y a 0 y 

(28) 

J = S C. C + C S 0 S, 
a 0 y a 3 Y 

(29) 

and 


K * C S 0 C - S C c S , 
a 3 Y a 3 Y 

(30) 

where C a = cos k a, S o = sin h a, C g * cos % $, etc. 
(30) for a, 3, y results in 

Solving equations (27) • 

H (a + y) = tan* 1 [(J + I )/ (H + K)], 

(31) 

H ( a - Y ) ■ tan" 1 [(J - I )/ (H - K)L 

(32) 

3 « tan" 1 [2 (HK + IJ)/^ + Y*], 

(33) 

where 


X = H 2 - I 2 + J 2 - K 2 , 

(34) 

and 


Y * 2 (JK - HI). 

(35) 


As reflected by the flow chart shown in Figure 14, Equations (31) - (35) are 
used to compute the Euler angles when the 'Q231' entrypoint is used. 

For a pitch, yaw, roll sequence, the orientation versor and the Euler 
angles are related by the equations 


The Euler-angle solutions In this case are 


>s (o - Y ) ■ tan' 1 [(J - K)/(H + I)], (40) 

>5 (a + Y ) » tan' 1 [(J + K)/(H - I)], (41) 

6 = tan' 1 [2 (HI - JK)/ ^777], (42) 

where 

X = H 2 - I 2 + J 2 - K 2 (43) 

and 

Y = 2 (IJ + HK). (44) 


It can be seen that Equations (31) - (35) become Identical with (40) - (44) 
when 

(a) -K is substituted for I, 

(b) I is substituted for K, and 

(c) the sign of y is reversed in the first set of equations. 

This characteristic of the two Euler sequences allows common logic to be used 
for the major portion of the computations, as illustrated in Figure 14. 
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Figure 14. *Q231 ’/ *Q213* Subroutine Logic Flow 







2.2.7 *0313* 


Given an orientation versor q FG , the *Q313 ' subroutine calculates a set of 
Euler angles (a.e.y) that corresponds to a rotation of coordinate system F 
Into coincidence with system G, assuming F Is rotated sequentially about Its 
Z, X, Z axes. 

2.2.7. 1 Arguments 

pi = Address of the orientation versor. 

p2 = Address of the output array of Euler angles (measured in radians, and 
stored In the same order as the logical rotation sequence). 

2. 2. 7. 2 Example of Usage 

Suppose that coordinate system F Is a geocentric Inertial system defined 
such that the first point in Aries lies on the Xp axis, and the earth's 
angular momentum vector lies along the Zp axis. Let G be another geocentric 
system defined such that the Instantaneous position of some satellite lies 
on the X 6 axis, and such that the angular momentum vector of the satellite 
CG lies along the Zq axis. Let the versor qpg reside in rO, rl, r2, r3. Then 
the statement c // *Q31 3 1 (0,4) would cause the right ascension of the ascending 
node of the satellite orbit to be stored in r4, the inclination of the orbit 
to be stored In r5, and the argument of latitude corresponding to the satellite 
position to be stored in r6. 

2. 2. 7. 3 Computations 

BBI * A 

The Euler angles are related to the versor q FG = H + i I + j J + k K by 
the equations 
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H = COS *5 (a + y) cos H 6* (45) 

I » cos h (a - y) sin h 6. (46) 

J ■ sin H (a - y) sin H 3* (47) 

and 

K = sin H (a + y) COS H 3 (48) 

(see Reference 10). It follows that 

J5S = tan" 1 (vTI+JJ, \/RH+i(l(), (49) 

h (a + y) * tan 1 (K,H), (50) 

and 

H (o - y) = tan 1 (J,I). (51) 


A flow chart of the 'Q313' logic Is shown In Figure 15. 
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Figure 15. *Q31 3 1 Subroutine Logic Flow 



2.2.8 


The '231Q'/'213Q' subroutine performs the Inverse function of the 'Q23V/ 
'Q213' subroutine described In Section 2.2.6. That Is to say, given an order- 
ed set of Euler angles a, P, y, it calculates the associated orientation versor 
q"pg. The ' 231Q ' entrypolnt Is used when the rotation sequence Is pitch, yaw, 
roll, and the * 21 3Q ' entrypolnt Is used when the sequence Is pitch, roll, yaw. 

2. 2. 8.1 Arguments 

pi = Address of the input array of Euler angles (measured in radians, and 
stored In the same order as the logical rotation sequence). 

p2 = Address of the orientation versor. 

2. 2. 8. 2 Example of Usage 

If rl , r2, r3 contain the pitch, yaw, and roll angles that define a 
certain orientation of the Orblter's body axes B relative to its local -vertical 
frame G, then the statement c U * 231Q 1 ( 1 ,0) In the calling routine would cause 
the orientation versor q^g to be stored In rO, rl , r2, r3 (thus destroying the 
Input values of the Euler angles, plus whatever value may have been stored In 
rO before calling the subroutine). If the variable A has been assigned the 
value 100, and If r4, r5, r6 contain the pitch, roll, yaw angles that define 
the orientation of system G with respect to another system F, the statement 
c£l l 213Q'(4.A+10) would cause the versor q F g to be stored in rllO, rill, rl 1 2 , 
rl 1 3 . 


2. 2. 8. 3 Computations 

The computational sequence defined by the flow chart In Figure 16 Is a 
straightforward Implementation of Equations (27) - (30) and (36) - (39) from 
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Figure 16. '231Q7'213Q' Subroutine Logic Flow 







2.2.9 '3130' 

The '313Q' subroutine performs the Inverse function of the 'Q313' subroutine 
described In Section 2.2.7. Given an ordered set of angles a, $, y that 
defines the orientation of a coordinate system G with respect to another system 
F In terms of sequential rotations of F about Its Z, X, Z axes, the '313Q' sub- 
routine computes the associated versor qp G 

2.2.9. 1 Arguments 

pi = Address of the input array of Euler angles (measured In radians, and 
stored In the same order as the rotation sequence). 
p2 * Address of orientation versor. 

2. 2. 9. 2 Example of Usage 

Let the F and G coordinate systems be defined as In Section 2.2.7. Let 
the right ascension of the satellite orbit reside In r4, the inclination of the 
orbit In r5, and the satellite's argument of latitude in r6. Then the state- 
ment cjU l 313Q l (4,0) would cause the orientation versor q FG to be stored in 
rO, rl, r2, r3. 

2. 2. 9. 3 Computations 

The computational sequence defined by the flow chart in Figure 17 Is a 
straightforward Implementation of Equations (45) - (48) from Section 2.2.7. 
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Figure 17. 1 31 3Q ' Subroutine Logic Flow 
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2.2.10 ' 1MATQ * 


Given the Inverse 


M' 1 • [r] t ■ ft] 


11 

T 12 

T 13 

21 

T 22 

T 23 

31 

T 32 

T 33 


( 52 ) 


of the coordinate transformation matrix [R] defined by Equation (26) In Section 
2.2.4, the * IMATQ 1 subroutine computes the associated orientation versor 

q FG * q 0 + i Pi + i q 2 + * q 3* (53) 


2.2.10.1 Arguments 

pi * Address of the Inverse transformation matrix, stored (columnwise) In 
the order T^, T 21 , T 31 , T 12 , T 22 , T 32 , T 13> T 23 , T 33 - 

p2 3 Address of the output versor. 

2.2.10.2 Example of Usage 

Let F represent an ECI coordinate system, defined such that the first 
point In Arles lies on the Xp axis, and the earth's angular momentum vector Is 
aligned with the Zp axis. Let G represent another geocentric coordinate 
system defined such that the axis Is aligned with the geocentric position 
vector ft of a satellite, and such that the Zg axis Is aligned with ft x 7 , where 
7 Is the geocentric Inertial velocity vector of the satellite. Suppose that 
ft and 7 are known In terms of their components In F, and that it is desired to 
evaluate cfpg. 

The inverse transformation matrix can be expressed as 
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where Hp X » ?p/|lTp| . Up* • (ftp x Vp)/|ftp x ?p| , and U^ Y ■ U^ 2 x Ujl x . If ftp re- 
sides In rl3, rl4, rl5, the statement ctt •SXV'nx/C D0TP'n3J3)) J3.0) will 
cause oS x to be stored In rO, rl, r2. If Vp resides In rl6, r!7, r18, the 
statement c // 'CRSP’(13.16.6): c jj •SM , C\Aj\ , DOTP'l6.6)).6.e) will cause U® 2 to 
be stored In r6, r7, r8, and the statement cil ’CRSP'tS.O.S) will then cause 
Il£ Y to be stored In r3, r4, r5, thus completing the columnwise storage of [Tj 
In rO through r8. The statement c // * IMATQ ' (0 .9) will then cause q FG to be 
stored In r9, rlO, rll, rl2. 


2.2.10.3 Computations 


By comparing Equations (26) and (52), It can be seen that [T] can be 
expressed. In terms of the versor components appearing In Equation (53), as 


[T] 


(q 0 +q l 

2(q,V‘)(fts ) 

2(q, q 3 -q 0 q 2 ) 


2(qiq 2 -q 0 q3 ) 

<q§-qf + q|-q|) 

2(q 2 q3+q 0 qi) 


2(q 1 q 3 + q 0 q 2 ) 

2(q 2 q 3 -q 0 qi ) 

(qJ-qi-qf+qj) 


(54) 


2 2 2 2 

Taking cognizance of the facts that q Q + q^ + q 2 + qj = 1 » and that one (but 
only one) of the versor components can be given an arbitrary sign, four 
different solutions for the values of q^k » 0,1 ,2,3) can be found In terms of 
the Tjj(1,j * 1,2,3). These solutions are given by the equations 


q Q * >/l + + T 2 2 + T 33^ 2 


(54a) 


q l s ^ T 32 " T 23^ /4 q 0 


(54b) 


q 2 “ ^ T 13 " T 31* /4 q 0 

(54c) 

q 3 * ^ T 21 " T 12) /4 q 0 

(54d) 

q l " T n - T 22 - V 2 

(55a) 

q 2 * ^21 + ^12^ 4 q l 

(55b) 

q 3 “ ^ T 13 + T 31^ 4 q l 

(55c) 

q 0 * ^32 * T 23^ 4 q l 

(55d) 

V ^ “ *11 + T 22 “ T 33^ 2 

(56a) 

q 3 * ^ T 32 + T 23^ 4 q 2 

(56b) 

q 0 * ^ T 13 ' T 31^ 4 q 2 

(56c) 

q l * (T 21 + T ]2 )/4 q 2 

(56d) 

"3 ‘ S * 'll - t 22 + V 2 

(57a) 

q 0 ‘ * f 2) * T 12^ /4 q 3 

(57b) 


* ^ T 13 + T 31^* ^3 


(57c) 


C 


• * 


^ ^2 " ^ T 32 + T 23 ^* <, 3 ‘ 

Since any versor component (or as many as three of them) may turn out to 

be zero, numerical tests are necessary to determine which one of the alternate 

computational sequences defined by the preceding equations can be used In a 

2 2 2 2 

given case. It Is known from the versor Identity < 1 q ♦ dj + ^ ♦ <13 - 1 that 
at least one component must have a magnitude as great as *s, and that none can 
* have a magnitude greater than unity. As defined by the flow chart shown In 

Figure 18, the 'IMATQ' logic tests the magnitude of each versor component In 
sequence until one is found that has a magnitude at least as great as As 

< soon as a component q N Is found that satisfies the inequality 4 q^ >_ 1, it Is 

arbitrarily given a positive sign, and the remainder of the components are 
calculated by using the equations that contain 4 q N In the denominator. 





m 





i 
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Figure 18. * IMATQ 1 Subroutine Logic Flow 
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2.3 STATE VECTOR DERIVATIVES SUBROUTINE ('DERIVS') 


The 'DERIVS' subroutine computes the first derivative (with respect to time) 
of each quantity contained In an extended array of state variables that In- 
cludes, In addition to the translational and rotational state vectors defined 
in Section 1.2, such things as propellant expenditures and control torque 
Integrals. The memory allocation table contained in Appendix D lists the 
specific quantities that are included In the Shuttle state variable array 
(registers r25~r49) and the payload state variable array (registers r50-r74). 

Containing mathematical models of the orbital flight environment and of 
all vehicle systems affecting the dynamics of the Shuttle and the payload, the 
'DERIVS' subroutine represents the heart of the #TRAJ processor. 

2.3.1 Input Data 

2. 3. 1.1 Argument List 

The 'DERIVS' argument list consists of a single parameter (pi) which is 
assigned a value of 0 or 1 by the calling routine according to whether state 
variable derivatives are to be computed for the Shuttle or the payload. 


2. 3. 1.2 r-Registers 


The contents of the following r-registers (see Appendix D ) are used as 
input quantities in the 'DERIVS' calculations: 


Shuttle Computations 


(Pi 


Compi 


r23, r24 
r25-r39 
r75-r86 
r91-r99, r!24 


Payload Computations 
(Pi = 1) 

r23, r24 

r50-r64, r74 

rlOO-rlll 

r!16-rl24 
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2. 3. 1.3 Array Variables 

The appropriate 12x12 RCS/DAP response matrix (see Section 3) must reside 
In the HPL array variable R[*] (represented symbolically as [R] In the following 
flow charts) before calling ' DERI VS ' for Shuttle derivatives (i.e., with pi * 0). 

When SRM thrusting is to be simulated, the name of the disk file containing 
the appropriate SRM thrust table must have been assigned to file no. 3 and the 
array variable A[*] (see Figure 19d for contents) must have been initialized 
before calling DERIVS for payload derivatives (I.e. , with pi = 1). The contents 
of the A[*] array (symbolically, [A]) are updated as required in 'DERIVS' by 
reading new thrust-profile coordinates from file no. 3 as the burn progresses. 

2.3.2 Output Data 

As indicated in Appendix D , the r-registers have been allocated so that 
the address of every state derivative can be found simply by adding 100 to the 
address of the corresponding state variable. Thus, Shuttle state derivatives 
are stored in r!25-rl49 when pi = 0, and payload state derivatives are stored 
in registers rl50-rl74 when pi = 1 . When pi = 0, the contents of r285 
(normally zero) will be set equal to 1 if 'DERIVS' finds the Shuttle RCS 
inadequate to maintain the commanded attitude. 

Sometimes the calling routine requires, in addition to the state variable 
derivatives, the angular velocity vector (wg G or ) of the local-vertical 
coordinate system. This vector resides in the volatile registers rl-r3 when 
execution control is returned from 'DERIVS' to the calling routine. Another 
vector sometimes needed by the calling routine is the linear acceleration 
vector (A G or ? g ), which resides in r7-r9 upon the return from 'DERIVS'. 
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2.3.3 Computations 


The computations performed by 'DERI VS' are defined by the flow chart shown 
in Figure 19a through 19g. For the most part, the computations are described 
in terms of logical symbols (as opposed to the HPL variable names and addresses 
that appear in the code.) The major portion of the computational logic applies 
equally to Shuttle and to payload calculations. To minimize confusion in 
following the flow of the logic, Shuttle symbols (e.g., fig rather than a^) are 
used almost exclusively throughout the flow chart, even in those sections of 
the logic that apply uniquely to the payload. The only exceptions to this 
general rule occur when there is no Shuttle-related equivalent of a payload- 
related variable (e.g., i^ + , / 2+ , etc.). 

Insofar as practicable, the order of computation in the HPL code (Section 
C.l) is arranged to follow that shown in the flow chart. Even so. Appendix D , 
which shows the correlation between logical symbols and r-register numbers, is 
indispensable to the understanding of the code. 

Following is a summary of all the data registers that will or may be 
modified by an execution of ' DERI VS 1 : 

Shuttle Computations Payload Computations 

(Pi ■ 0) (Pi s D 


r0-rl8 

r87-r90 

rl25-rl49 


r0-rl8 
r 1 1 2-rl 1 5 
rl50-rl74 


r285 


A[5]-A[8] 


A-Z 


A-Z 
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Figure 19b. * DERI VS * Subroutine Logic Flow 
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Figure 19c. 'DERIVS' Subroutine Logic Flow (cont.) 
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2.4 MAIN LOGIC FLOW 


The "main logic" of the #TRAJ link consists of nothing more than an 
Instruction that causes the %RMAT link (Section 3) to be appended to #TRAJ. 
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3. RESPONSE MATRIX COMPUTATION LINK (XRMAT) 


The function of the XRMAT link Is to compute matrices that describe the 
response of the Shuttle to OAP/RCS commands, and to store the matrices In 
designated disk files for subsequent use during trajectory Integration. 

3.1 INPUT DATA 

3.1.1 Jet Force Table 

Basic Inputs to the response matrix computations Include the data shown In 
Table 1, which reside In the disk files named "$JFT" and "SJFTM". The "$JFTM H 
file contains the mnemonic Identification symbols (shown In the second column 
of Table 1) for the Orblter's 44 RCS jets, as defined In Figure 4. 2. 2. -2 of 
Reference 11. The "$JFT M file contains the body-fixed thrust components and 
the station coordinates of the thrust application point (which together define 
a torque component normal to the thrust line) for each thruster, along with a 
factor C (In the last column of Table 1) whose product with the thrust vector 
defines a component of torque parallel to the thrust line. This latter com- 
ponent of torque (parallel to the thrust line) Is the result of RCS jet plume 
Impingement on the exterior surfaces of the Orblter. Table 1 results from 
combining the basic RCS jet data shown In Table 2 with the plume Impingement 
data shown In Table 3, as explained In Reference 9. 

3.1.2 Jet Select Tables 

For each flight profile segment, the HFRMP user must specify which one of 
three basic combinations of thrusters Is to be used for attitude and/or trans- 
lational control of the Shuttle. The available options are designated V 
(vernier jets), P (primary jets), and PZI (primary jets with +Z thrusters 
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Table 1. MS Thruster Data, Disk Files "$JFT" and "SOFTM" 




306.72 

306.72 

306.72 

362.67 

364.71 

362.67 

364.71 


360.62 

350.93 



1515.41 

1528.43 

1541.45 

1554.47 

1515:41 

1528.43 

1541.45 

1554.47 

1516.36 

1529.23 

1542.10 

1516.36 

1529.23 

1542.10 

15 ) 0.34 

1526.50 

1542.65 

15 ) 0.34 

1526.50 

1542.65 

324.35 

324.35 

1565.00 

1565.00 

1565.00 

1565.00 


• 14.65 

• 69.50 

- 71.65 


71.65 

14.39 

0.00 

- 14.39 

61.42 

• 61.42 

66.23 

• 66.23 

137.00 

124.00 
• 137.00 
- 124.00 
• 149.91 
• 149.91 
- 149.91 
- 149 . 9 ) 

149.91 

149.91 

149.91 

149.91 

• 115.43 

- 1 ) 5.43 

- 115.43 

115.43 

1 ) 5 . 43 . 

115.43 

- 104.34 

- 103.22 

- 102.11 

104.34 

103.22 

102.11 

59.70 

- 59.70 

149 .C 7 

- 149.87 

118.00 
• 118.00 


392.96 

394.45 
392.96 

373.73 
359.25 

373.73 
359.25 

413.46 
414.53 
4 ) 3.46 

356.95 

356.95 

358.44 

358.44 

473.06 

473.06 

473.06 

473.06 

452.87 

452.88 

452.88 

452.88 

452.87 

452.88 

452.88 

452.88 

481.95 

481.95 

481.95 

481 .95 

481 .95 

481.95 

431.12 
434.78 

438.44 

431.12 
434.78 

438.44 

350.12 

350.12 

459.00 

459.00 

455.44 

455.44 


.0000 
• 0.2421 
- 0.2759 
- 0.3098 
- 0.3437 
0 . 242 ) 
0 . 27 S 9 
0.3098 
0.3437 
- 0.1586 
- 0.0648 
0.0290 
0.1586 
0.0648 
- 0.0290 
0.1103 
- 0.0959 
- 0.3320 
• 0.1103 
0.0959 

































o 


Table 2. Basic RCS Thrust Data (Without Plume Impingement) 



liiil 


mm 

III 

mm 

sYa 

(IN) 

il 

(IN) 

ML 

(IN) 

1 

F2F 

-879.4 

mn 

m 

306.72 

14.65 

392.96 


F3F 

•879.5 



306.72 

0.00 

394.45 


F1F 

-879.4 

26.2 

MHXl 

306.72 

-14.65 

392.96 


Fit 

•26.3 

873.6 

18.2 

362.67 

-69.50 

373.73 



-21.0 

870.3 

0.5 

364.71 

-71.65 

359.25 



*26.3 

-873.6 

18.2 

362.67 

69.50 

373.73 



•21.0 

-870.3 

0.5 

364.71 

71.65 

359.25 

8 

F2U 

-32.3 

-11.7 

874.4 

350.93 

14.39 

413.46 

9 

F3U 

-31.9 


873.5 

350.92 

0.00 

414.53 

10 

F1U 

-32.3 

11.7 

874.4 

350.93 

-14.39 

413.46 

11 

F2D 

•28.0 

-616.4 

-639.5 

333.84 

61.42 

356.95 

12 

FID 

-28.0 

616.4 

-639.5 

333.84 

-61.42 

356.95 

13 

F4D 

-24.8 

-612.6 

•639.4 

348.44 

66.23 

358.44 

14 


-24.8 

612.6 

-639.4 

348.44 

•66.23 

358.44 

IS 

R3A 

856.8 

0.0 

151.1 

1555.29 

137.00 

473.06 

16 

R1A 

856.8 

0.0 

151.1 

1555.25 

124.00 

473.06 

17 

L3A 

856.8 

0.0 

151.1 

1555.29 

-137.00 

473.06 

18 

L1A 

856.8 

0.0 

151.1 

1555.29 

-124.00 

473.06 

19 

L4L 

0.0 

870.5 

-22.4 

1516.00 

-149.87 

459.00 

20 

La 

0.0 

870.5 

•22.4 

1529.00 

-149.87 

459.00 

21 

L3L 

0.0 

870.5 

-22.4 

1542.00 

-149.87 

459.00 

22 

L1L 

0.0 

870.5 

•22.4 

1555.00 

-149.87 

459.00 

23 

R4R 

0.0 

-870.5 

-22.4 

1516.00 

149.87 

459. 0C 

24 

R2R 

0.0 

-870.5 

-22.4 

1529.00 

149.87 

459.00 

25 

R3R 

0.0 

-870.5 

-22.4 

1542.00 

149.87 

459.00 

26 

R1R 

0.0 

•870.5 

-22.4 

1555.00 

149.87 

459.00 

27 

L4U 

0.0 


870.0 

1516.00 

-132.00 

480.50 

28 

L2U 

0.0 

Hlfl 

870.0 

1529.00 

-132.00 

480.50 

29 

LIU 

0.0 

0.0 

870.0 

1542.00 

-132.00 

480.50 

30 

R4U 

0.0 

0.0 

870.0 

1516.00 

132.00 

480.50 

31 

R2U 

0.0 

HX1 

870.0 

1529.00 

132.00 

480.50 

32 

R1U 

0.0 


870.0 

1542.00 

132.00 

480.50 

33 

L40 

170.4 

291.8 

r801.7 

1516.00 

-111.95 

437.40 

34 


170.4 

291.8 

-801.7 

1529.00 

-111.00 

440.00 

35 


170.4 

291.8 

-801.7 

1542.00 

-110.06 

442.60 

36 

R4D 

170.4 

-291 .8 

-801.7 

1516.00 

111.95 

437.40 

37 

R20 

170.4 

-291 .8 

-801.7 

1529.00 

111.00 

440.00 

38 

R3D 

170.4 

-291.8 

-801.7 

1542.00 

110.06 

442.60 

39 

F5R 

-0.8 

-17.0 

-17.6 

324.35 

59.70 

350.12 

40 

F5L 

-0.8 

17. U 

-17.6 

324.35 

-59.70 

350.12 

41 

R5R 


-24.0 


1565.00 

149.87 

459.00 

42 

L5L 


24.0 

mm 

1565.00 

-149.87 

459.00 

43 

R5D 



-24.0 

1565.00 

118.00 

455.44 

44 



■a 

-24.0 

1565.00 

-118.00 

455.44 
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Table 3. Force and Moment Increments Due to Plume Impingement 
(CG @ STA 1076.7, BL 0, WL 375.0) 


THRUSTER 

ID 


A FORCE (LB) 


A MOMENT (FT-LB) 


-TT~ 

wtm 

F z 

3.34 

-1.66 

| 27.17 

3.34 

-1.66 

27.17 

3.34 

-1.66 

27.17 

3.34 

-1.66 

27.17 

0.09 

76.30 

2.21 

0.09 

76.30 

2.21 

0.09 

76.30 

2.21 

40.19 

26.59 

225.66 

40.19 

26.59 

225.66 

40.19 

26.59 

225.66 

3.34 

1.66 

27.17 

3.34 

1.66 

27.17 

3.34 

1.66 

27.17 

3.34 

1.66 

27.17 

0.09 

-76.30 

2.21 

0.09 

-76.30 

2.21 

0.09 

-76.30 

2.21 

40.19 

-26.59 

225.66 

40.19 

-26.59 

225.66 

40.19 

-26.59 

225.66 


-795.30 

-795.30 

1859.90 

1859.90 

1859.90 

-2475.50 

-2475.50 

-2475.50 

795.30 

795.30 

795.30 

795.30 

-1855.90 

-1859.90 

-1859.90 

2475.50 

2475.50 

2475.50 


762.42 

762.42 

762.42 

762.42 

94.21 

94.21 

94.21 


8469.00 

8469.00 

8469.00 


762.41 

762.41 
762.4; 
762,4: 


8469.00 

8469.00 

8469.00 


(Taken from Reference 13) 










Inhibited). Corresponding to each of these options Is a jet-select table 
(Tables 4-6) which Identifies the particular jet or combination of jets that 
In .0 be fired In response to each of the six RCS translational acceleration 
commands (+X, -X, +Y, -Y, +Z, -Z) and the six rotational acceleration commands 
(+ROL, -ROL, +PCH, -PCH, +YAW, -YAW). The jet select tables, which reside In 
the disk files named "$JSV", "$JSP", and "$JSPZI", are not routinely available 
for modification by the user. However, the HFRMP software system Includes a 
jet-select editing processor that makes It possible to update the tables or to 
provide additional options with little difficulty. 

As Indicated In Table 4 by the absence of any jet designations for the 
execution of translation commands, the V option (vernier jets) can be used only 
for rotational control. The P option (Table 5) and the PZI option (Table 6) 
can be used for translational and/or rotational control. 

In the PZI option, no jets are fired that would expel propellant directly 
upward with respect to the Orblter body. Translational acceleration In tre 
downward direction. If commanded. Is achieved (at a comparatively high propel- 
lant cost) by firing the +X and -X thrusters simultaneously. The cant angles 
of the +X and -X jet thrust lines produce a small net acceleration In the +Z 
(downward) direction. This option normally Is used only when the Orblter Is 
maneuvering In the near vicinity of a payload that must be protected from jet 
plume Impingement. 

3.1.3 Shuttle Mass Properties 

At %MAT execution time, the user-defined mass properties of the Shuttle* 


+ Wh1ch remain constant during HFRMP trajectory Integration. 
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reside In the disk file named "1*", where they are stored by the data-base 
editing processor #DBED. The "I* 11 file contains 32 scalar quantities, of 
which XW1AT makes use of the following: 


Item No. 
1 


2 

! XX 

3 

X YY 

4 

hi* 

5 

hi 

6 

hi ' 

7 

! xy j 

8 

STA' 

9 

BL • 

10 

WL 

• 

• 


32 

XRMAT 


Description 

Shuttle gross weight (1b) 


moments of 
Inertia 


products of 
Inertia 


(slug-ft*) 


station coordinates 
of CG (In) 


Item 32 In the "1*" file Is a flag that Is tested Immediately after entry Into 
XRMAT to determine whether a re-computation of the response matrices Is 
necessary. It Is set equal to zero during Initialization of the program disk, 
and thereafter (by the appropriate editing processor) when the contents of 
"1*", "$JFT", or any jet-select table are changed In any way. It Is set equal 
to 9 by XRMAT upon completion of the computation and storage of the response 
matrl ces . 


3.2 OUTPUT DATA 

The output of the XRMAT link consists of eight 12x12 matrices, each corre- 
sponding to a particular jet select table and a particular mode of RCS cross- 
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coupling compensation, and each stored In a separate disk file as Indicated In 
Table 7. The compensation modes available to the HFRMP user are designated 
NONE (no compensation), ROT (rotational compensation only), and FULL (rota- 
tional and translational compensation). In apy given flight profile segment, 
the HFRMP user may choose any jet-select option In combination with any com- 
pensation mode, except V with FULL. Full compensation Is Impossible with the 
vernier jets, and the compensation mode Is Internally defaulted to ROT when 
such a combination Is specified by the user. 

Tables 8-10 show typical response matrices for the P jet-select option 
and Increasing degrees of cross-coupling compensation. For purposes of Illus- 
tration, each matrix was transposed and then partitioned Into two 12x6 matrices 
so It could be printed conveniently on a single page. Each column of the 
(untransposed) response matrix corresponds to a particular translational or 
rotational command. Rows 1-3 contain the body-axis components of the steady- 
state linear acceleration, and rows 4-6 contain the components of angular accel- 
eration. Rows 7-12 contain RCS propellant consumption rates, broken down 
according to source (forward, aft left, or aft right tanks) and control function 
(translation or rotation). The computations that produce the results shown In 
Tables 8-10 are described by the flow chart shown In Figures 20a through 20g, 
wherein the symbols [U], [R], and [F] represent the uncompensated, rotatlonally 
compensated, and fully compensated response matrices. 

3.3 COMPUTATIONS 


It will be noted that six small rectangles, each enclosing a pair of 
acceleration components, run diagonally across the upper partition of the 
matrix shown In Figure 8. The components thus enclosed represent the uncompen- 
sated "principal response" of the Shuttle to each of the 12 translation and 
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Transposed and partitioned for purposes of illustration. 



0 

S 

u> 

1 

s 

M 

\ 

I 

£ 

iilisSSSSSsS 
• •••••••••*• 

OOf MOOWWmMWO 

5| 

5" 

0.04949 

0.30190 

2.62193 

1.66496 

0.24650 

0.13175 

3.42754 

3.51573 

3.13352 

3.12431 

0.62176 

3.72896 

si 

0.09246 
0.59096 
1.15363 
1.15903 
0.23295 
C. 18607 
0.41495 
0.40939 
6.21529 
3.10772 
3.12941 
3.13013 

I 

M 

N 

i 

3.10710 

0.00000 

0.00000 

3.10710 

3.10710 

6.21420 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

5| 

-* 

§ 

3.10710 

0.00000 

3.10710 

0.00000 

3.10710 

6.21420 

0.00009 

0.00000 

0.00000 

0.00000 

0.00000 

0.00000 

si 

0.00000 ' 
6.21420 
3.10710 
3.10710 
3.10710 
6.21420 
0.00000 
0.00000 
C. 00000 
0.00000 
0.00000 
0.OOCOC 

CM) 

+ *+•+! 
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Transpest? and partitioned for purposes of illustration. 




Purposes of illustration. 
















I 


( 


; ( 



r ‘ 




rotation coimands. All other components of acceleration represent "extraneous 
responses" , commonly referred to as "cross-coupling effects"* that arise 
from the canting of certain jet centerlines away from the body axes* uneven 
moment arms, and plume Impingement on the exterior surfaces of the Prblter. 

The details of the uncompensated response matrix calculations are defined In 
Figures 20b and 20c. 

The process of cross-coupling compensation can be explained in the follow- 
ing general terms. Assume that the jets activated by a given "primary command" 
P fire continuously for a long period of time at. The effects of extraneous 
accelerations are nullified by Intermittent firings of jets activated by 
"compensating conmands". Let the accumulated firing time of the jets activated 
Intermittently by any particular compensating command K be represented by stg. 
The ratio Y K ■ jt K /at Is referred to as the "duty cycle" of the (jets activated 
by the) compensating command. 

The uncompensated response to the primary command (a column In the uncom- 
pensated response matrix) we represent with the symbol [Ul p = [(U^ p ) (U 2 p ) ... 
(Uj 2 ,p)] T * The compensated response (a column In the compensated response 
matrix Is defined by the expression. 


K 

[C] p » [U] p + z y k [u] k * 


where the summation Includes all the compensating commands that are required 
to nullify the extraneous accelerations. 

Rotational compensation, which Involves the nullification of only the 
extraneous angular accelerations (see Table 9). Is carried out In two phases. 
The rotation conmands themselves are rotational ly compensated In the first 
phase (Figure 20d), which Is of necessity an Iterative process. Successive 







s 

i 


1 


90 



compensations introduce new extraneous angular accelerations, whose magnitudes 
must diminish progressively If convergence Is to be realized. Non-convergence 
of this process Implies that (given the Input mass properties and the jet- 
select table under consideration) the attitude of the Shuttle can not be 
controlled, and an appropriate warning message Is output. 

The translation commands are rotatlonally compensated In the second phase 
of rotational compensation (Figure 20e), This Is a non-iterative process that 
Is simplified by use of the results from the first phase. 

Full compensation Is achieved by nullifying the extraneous linear acceler- 
ations that remain In the rotatlonally compensated matrix. Again, the process 
Is carried out In two phases. Starting with the rotatlonally compensated 
matrix as an Input (which eliminates the need to worry about compensation 
commands producing extraneous angular accelerations) the first phase (Figure 
20f) Is devoted to the translational compensation of the translation commands. 
This Is an Iterative process that Is essentially Identical to that defined In 
Figure 20d. Likewise, the second phase of translational compensation (Figure 
20g) Is essentially Identical to that shown In Figure 20e. 

The computations are repeated for each of the three jet select tables, 
except that full compensation Is not attempted with the vernier jets. Appro- 
priate warning messages (if any) are stored in the desk files along with the 
response matrices. 

The last executable statement In the %RMAT link causes the code of the 
%TNIT link (Section 4) to be appended to that of the base link (#TRAJ), In the 
region of computer memory formerly occupied by the %RMAT code. The HPL get 
command Is used for this purpose, which causes all data registers used by %RMAT 
to be de-allocated (erased) before execution control Is handed over to %TNIT. 
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Figure 20a. fRMAT Logic Flow 







































Figure 20e. L.RMAT Logic Flow (cent.) 
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Figure 20f. %RMAT Logic Flow (cont.) 
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4. TRAJECTORY INITIALIZATION LINK (ZTNIT) 

The function of the %TNIT link Is to Initialize calculator memory In prep* 
aratlon for trajectory Integration. Elements of the user-supplied data base 
(Appendix E) are read from appropriate disk files (where they will have been 
stored by a previous execution of the #DBED processor) and processed as neces- 
sary to define the constants and the Initial state of the Shuttle/payload system. 

4.1 FUNCTION SUBPROGRAMS 


4.1.1 'JD* 

The ' JD' function subprogram computes the Julian day number corresponding 
to a given year, month, and day of the Gregorian (civil) calendar. 

4. 1.1.1 Argument List 


pi = year 
p2 8 month 
p3 8 day 


integers 


4. 1.1. 2 Example of Usage 

The Instruction 'JD'nSBO.A^) -» D would cause D to be assigned the value 
2444332 (the number of the Julian day commencing at Greenwich noon on 2 April 
1980). 


4. 1.1. 3 Computations 

Figure 21 Is a flow chart of the 'JD* computational procedure, which was 
taken from Reference 14. 
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Figure 21. 1 JD* Function Subprogram Logic Flow 
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4.1.2 * PLOTYP * 



Given the user-assigned mnemonic symbol which designates the desired type 
of data plot, 'PLOTYP* returns the appropriate numeric code for Internal use. 

4. 1.2.1 Argument List 


None. The mnemonic symbol must be assigned to the string variable P$ 
before 'PLOTYP' Is executed. 


4. 1.2. 2 Example of Usage 

If the character string "CPLV" resides In P$, the Instruction 'PLOTYP* ■» I 
will cause I to be assigned the value 2. 

4. 1.2. 3 Computations 
See Figure 22. 



Figure 22. ' PLOTYP ' Function Subprogram Logic Flow 
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4.1.3 ’UNIT* 


Given the user- as signed mnemonic symbol which designates the unit of 
distance for data plots, the 'UNIT' subprogram returns the appropriate conver- 
sion factor (the number of feet In the designated unit). 

4. 1.3.1 Argument List 

None. The mnemonic symbol must be assigned to the string variable P$ 
before 'UNIT 1 Is executed. 

4. 1.3. 2 Example of Usage 

The Instructions " NMI M - Pi: 'UNIT' * C would cause C to be assigned the 
value 6076.115. 

4. 1.3. 3 Computations 
See Figure 23. 


103 



04 








4.1.4 ‘MONTH 1 

Given the abbreviated name of a month, the 'MONTH' function subprogram 
returns the appropriate month number. 

4. 1.4.1 Argument List 

None. The abbreviated month name must be assigned to the string variable 
P$ before MONTH Is executed. 

4. 1.4. 2 Example of Usage 

If the character string " AUG" resides In P$, the Instruction 'MONTH' •> M 
will cause the number 8 to be assigned to the variable M. 

4. 1.4. 3 Computations 
See Figure 24. 






Given the size of a square matrix stored columnwise In a one-dimensional 
array, and the row and column Indices of one of Its elements, the 'C' function 
returns the relative address (In the one-dimensional array) of the matrix 
element. 

4.1. 5.1 Argument List 

pi 3 Row Index of matrix element. 

p2 = Column Index of matrix element. 

p3 3 Matrix size (number of rows 3 number of columns). 

4.1. 5. 2 Example of Usage 

Assume the 3x3 matrix [T] Is stored columnwise In registers r8 - rl6, as 
shown below: 


Relative 

Address 

Register 

No. 

Matrix 

Element 

0 

r8 

T l.l 

1 

r9 

T 2,l 

2 

rlO 

T 3,l 

3 

rll 

T 1 ,2 

4 

rl2 

T 2,2 

5 

r13 

T 3,2 

6 

rl4 

T 1 »3 

7 

rl5 

T 2,3 

8 

r!6 

T 3,3 


The Instruction r(8 + 'C'd, 2,3)) -» A would cause the value of 2 t0 be 
assigned to the variable A. 
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4.1.6 jr 

Given the row and column Indices of an element In a symmetric matrix 
whose upper triangular form* Is stored columnwise In a one-dimensional array* 
the *T' function returns the relative address (In the one-dimensional array) 
of the matrix element. 

4. 1.6.1 Argument List 

pi » Row Index of matrix element. 
p2 = Column Index of matrix element. 

4. 1.6. 2 Example of Usage 

Assume the upper triangular form 

*1,1 *1,2 X l,3 

*2,2 *2,3 

*3,3- 

of the symmetric 3x3 matrix [I] Is stored columnwise In registers r10-r!5 as 
Indicated below: 


Relative 

Address 

Register 

No. 

Matrix 

Element 

0 

rlO 

*1.1 

1 

rll 

*1 ,2 sI 2,l 

2 

rl2 

*2,2 

3 

rl3 

1 1 ,3 8 *3.1 

4 

rl4 

*2,3 “ *3,2 

5 

rl5 

*3,3 
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The Instructions 10 » TO, 2) •» H: rH X would cause the value of I 
to be assigned to the variable X. 



4. 1.6. 3 Computations 


See Figure 26. 





4.2 MATRIX DIAGONALIZATION SUBROUTINE ( *DIAG* ) 

Given a syirmetrlc matrix [M] of order C, the * DIAG ' subroutine solves the 
equation 

[R] CM] [R] T • [0] 

for [R] and [D], where [D] Is a diagonal matrix (one whose off-diagonal elements 
are all zero, within some specified tolerance). 

4.2.1 Argument List 

pi * Order (size) of the Input matrix [M]. 

p2 * Address of the Input matrix [M], stored columnwise In upper tri- 
angular form. 

p3 ■ Address of the diagonalized output matrix [D], stored columnwise 
In upper triangular form. 

p4 « Address of the square (pi x pi) output matrix [R], stored columnwise. 
p5 = Tolerance on the maximum squared value of Dj where l^j (this Input 
Is optional; If not supplied by calling routine, It will be calculat- 
ed by 'DIAG'). 

4.2.2 Example of Usage 

Assume that the upper triangular form of the Shuttle's Inertia tensor 
[I]g (referenced to the body coordinate system B) Is stored columnwise In 
registers rlO - rl5 as Illustrated In Section 4. 1.5. 2. Then the Instruction 
ett 'DIAG'OJO.IO.l ) would cause [I] g to be replaced by [I] p (the same Inertia 
tensor but now referenced to the coordinate system P, whose axes are the 
principal axes of Inertia). Upon return from 'DIAG', the registers rl - r9 
would contain the 3x3 coordinate transformation matrix [R] that satlsflos the 


equations 


[R] ftp * ftg 
and 

[R]' ! Aj • [R] t Xj . Ip, 
where ft Is any arbitrary vector. 

4.2.3 Computations 

The method used to find the matrices [D] and [R] Is of the type known as 
the Jacobi method (Reference 15). This consists of a series of matrix rotations 
of the form 

[D] • (DUjCAJj • • • ) [H] ( • • • [R] 2 T [R], T ) 

where each [R] k Is* selected so as to cause one pair of off-diagonal elements 
to be zero after the k th rotation. An Initial trigger level, 6, Is computed 
from 

6 * yToT/c 

where a Is the sum of the squares of all off-diagonal elements and C Is the 
order of [M]. The square of each off-diagonal element of [N] is compared to 
6. When an element Is encountered whose square Is larger than 6, the matrix 
[N] Is rotated so as to cause that off-diagonal element to be zero. 

For example, suppose that the square of the (1,j)th element of [M] Is 
greater than 6, then [M] Is rotated by [R] 1 to get [M]' wh»re 

M'[1.J] • M'[j,1] » 0 
The matrix [R]j has the form 
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RjCl.l] ■ R^CJJ] - cos 4 
RjClJ] ■ - sin 4 

RjCJ.1] ■ sin 4 

and all other elements of [R]j are Identical to the unit matrix [1]. Rotation 
with [R]j, 1.e. t 

W ■ W, CM3 CRD| 

gives rise to the elements 

M'[1,l] - MCI .13 cos 2 4 + 2M[1 ,j] sin 4 cos 4 + 

M'[j,j] * MCI .13 sin 2 4 - 2M[1,j] sin 4 cos 4 + 

M'[1.j] - M'[j, 1 ] » (M[j ,J] - MC1.13) sin 4 cos 

+ M[ 1 ,j] (cos 2 4 - sin 2 4 ) 

from which It can be seen that 

M'Cl.j] ■ M'Cj.l] ■ 0 

provided 4 Is selected such that 

tan 2 4 ■ 2M[1 ,J]/(M[1 ,1] - M[J.j]) 

After this rotation, scan of the off-diagonal elements continues using 
now, however, [M]'. It Is true that any one rotation will cause other off- 
diagonal elements which were zero to become non-zero, but the trend Is to re- 
duce all off-diagonal elements to small numbers. 

When the scan of all off-diagonal elements has been completed, the trigger 
level 6 Is reduced by dividing It again by C, and the process Is repeated 


M[j,j] sin 2 4 
M[j.j] cos 2 4 
4 

t 
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until < Is 1«ss than a pre-set tolerance. The programmer may Include this 
tolerance In the argument list or he may let ‘DIAS' set the tolerance for him. 
This process Is shown In detail In Figures 27(a) and 27(b). 
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Figure 27a. 'DIAG' Subroutine logic Flow 












4.3 MAIN LOGIC 


The i'TNIT computations are described by the flow chart contained In 
Figures 28a through 28g. The memory allocation table In Appendix D will 
have to be consulted for the purpose of correlating the symbols appearing In 
the flow chart with rhe r- register numbers appearing In the HPL code, which 
Is contained In Appendix C.3. 
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Figure 28a. %TNIT Logic Flow 
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Figure 2Bb 


55TNIT Logic Flow (cont.) 





Figure 28c. %TNIT Logic Flow (cont.) 
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Figure 28d. %TNIT Logic Flow (cont.) 
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Figure 28e. %TNIT Logic Flow (cont.) 

f 







25 








BC*1 


BtH -•[<] 

■mi 



Figure 28g. %TNIT Logic Flow (cont.) 










S. ORBITER GEOMETRY LINK (XSSVU) 

The function of this link of the #TRAJ processor Is to plot an end and/or 
a side view of the Shuttle Orblter profile. It Is executed only when the user 
specifies RSBY (Rectangular Shuttle Body-fixed coordinates) for the plot type. 
In the graphics data file (see Appendix E.3). In addition, It draws dashed 
lines that represent the nominal limits of visibility (through the overhead 
and aft windows) from the on-orbit pilot's control position. The code (Appen- 
dix C.4) Is straightforward, for the most part consisting of HPL £lt commands 
followed by the station coordinates (measured In Inches) of the points that 
define the end and side profiles of the Orblter. Typical Orblter profile plots 
are Illustrated In Appendix H. 
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6. FLIGHT SEGMENT INITIALIZATION LINK (XSNIT) 

The %SNIT link Is called Into memory and executed at the beginning of 
every flight profile segment. Each time It Is executed. It reads a new flight 
segment definition (see Appendix F) from file number 1, to which the name "7*" 
has been assigned In the %TNIT link. 

After each new flight segment definition Is read from disk, appropriate 
processing of the new Input data Is performed In preparation for state variable 
propagation through the ensuing flight segment. The propagation Itself (In- 
volving numerical Integration of the equations of motion) Is performed by the 
%PR0P link (Section 7). Upon completion of the segment Initialization process, 
execution control Is passed to %PR0P by means of a chain Instruction that 
causes the %SNIT code to be replaced by that of %PR0P. 

Trajectory Integration Is terminated (l.e., %PR0P Is not called Into 
memory) when an end-of-flle mark Is encountered during the attempt to read a 
new flight segment definition. If a payload solid rocket motor (SRM) burn has 
been simulated, the Particle Impact Damage Integrator Processor (#PIDI) Is 
called Into memory Immediately after the termination of trajectory Integration. 
Otherwise, the user Is asked whether he wants to start a new run. If the 
answer Is yes (( ^ONTINUE )) , the Data Base Editor Processor (#DBED) Is called 
Into memory. If the answer Is no, ((©, ( CONTINU^ ) , program execution Is 
terminated by a st£ command. 

6.1 FUNCTION SUBPROGRAMS 

6.1.1 'KAM* 

Given the user-supplied mnemonic symbol that designates an attitude- 
maintenance option, 'KAM* returns the appropriate numeric code for Internal use. 
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6.1. 1.1 Argument List 

None. The mnemonic symbol must be assigned to the string variable P$ 
before 'KAN 1 Is executed. 

6. 1.1. 2 Example of Usage 

If the character string "LVRH" resides In P$, the Instruction 'HAM* K 
will cause K to be assigned the value 2. 

6. 1.1. 3 Computations 
See Figure 29. 
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6.1.2 * JSEL* 

Given the mnemonic symbol that designates a Shuttle RCS jet-select table, 
'JSEL* returns the appropriate numeric code for Internal use. 

6. 1.2.1 Argument List 

None. The mnemonic symbol must be assigned to the string variable P$ 
before 'JSEL' Is executed. 

6. 1.2. 2 Example of Usage 

If the character string " P" resides In P$ # the Instruction *JSEL' •» J 
will cause J to be assigned the value 1. 

6. 1.2. 3 Computations 
See Figure 30. 
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Figure 30 ' JSEL 1 Function Subprogram Logic Flow 
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6.1.3 * COMP 1 


Given the mnemonic symbol that designates a Shuttle RCS cross-coupling 
compensation option, 'COMP* returns the appropriate numeric code for Internal 
use. 


6. 1.4.1 Argument Hr t 

None. The mnemonic symbol must be assigned to the string variable P$ be- 
fore 'COMP' Is executed. 

6. 1.3. 2 Example of Usage 

If the character string "NONE" resides In P$, the Instruction 'COMP' C 
will cause C to be assigned the value zero (0). 

6. 1.3. 3 Computations 
See Figure 31. 
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6.1.4 'KTOMS' 


Given the mnemonic symbol that designates which (if any) of the Shuttle 
OMS engines are to be fired firing a given flight profile segment, 'KTOMS' 
returns the appropriate numeric code for Internal use. 

6. 1.4.1 Argument List 

None. The mnemonic symbol must be assigned to the string variable P$ 
before 'KTOMS' is executed. 

6. 1.4. 2 Example of Usage 

If the character string " L+R" resides in P$, the instruction 'KTOMS' -» K 
would cause K to be assigned the value 3. 

6. 1.4. 3 Computations 
See Figure 32. 
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6.1.5 ' KTRCS ' 

Given the mnemonic symbol that designates a Shuttle RCS thrust command, 
- 'KTRCS' returns the appropriate numeric code for Internal use. 

6. 1.5.1 Argument List 

None. The mnemonic symbol must be assigned to the string variable P$ 
before 'KTRCS' Is executed. 



6. 1.5. 2 Example of Usage 

If the character string "+R0L" resides In P$, the Instruction 'KTRCS' -» K 
would cause K to be assigned the value 7. 

6. 1.5. 3 Computations 
See Figure 33. 







6.1.6 'DATYP 1 


HFRMP Input data files contain three general types of data: non-integer 
numbers. Integers, and mnemonic symbols (character strings). Given the Item 
number of an entry In a flight segment definition (see Appendix F), 'DATYP' 
returns a numeric code which Identifies the type of data so that It can be 
printed In the appropriate format. 'DATYP' also loads the appropriate 
Identification text for the data Item Into the string variable B$, and the 
description of Its unit of measurement Into the string variable U$. 

6. 1.6.1 Argument List 

pi * Item number In the flight segment definition (see Appendix F). 

6. 1.6. 2 Example of Usage 

The instruction 'DATYP* (7) -» I would cause I to be assigned the value 1, 

the character string "SS XB RATE OR INCR " to be loaded 

Into B$, and the character string "DEG/SEC " to be loaded Into U$. 

6. 1.6. 3 Computations 

The 'DATYP' code, which Is straightforward and has no effect on trajectory 
computations, can be found In Appendix C.5. 
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6.2 SEGMENT DEFINITION LISTING SUBROUTINE ( 'SLIST* ) 

The 'SLIST ' subroutine Is used to list each flight profile segment defi- 
nition on the output line printer. Immediately after It Is read from the disk. 
In the format that Is Illustrated In Appendix F. The *SLIST' code, which Is 
straightforward and has no effect on trajectory computations, can be found In 
Appendix C.5. 
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The %SNIT computations are described by the flow chart contained In Figures 
34a through 34f. The memory allocation table In Appendix D will have to be 
consulted to correlate the logical symbols appearing In the flow chart with the 
r- register numbers that appear In the HPL code, which Is contained In Appendix 



Figure 34a. %SNIT Logic Flow 
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Figure 34b. %SNIT Logic Flow 
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Figure 34e. %SNIT Logic Flow (cont.) 
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Figure 34f. %S,NIT Logic Flow (cont.) 







7. STATE PROPAGATION LINK (%PROP) 

The function of the %PROP link Is to propagate the state of the Shuttle/ 
payload system through a single flight profile segment. Digital output data 
describing the state of the system are printed at the beginning and the end of 
the segment, and at user-specified Intervals between these two points. If the 
user has specified "RSBY" (Rectangular Shuttle Body-fixed coordinates) or "CPtV" 
(Curvilinear Payload Local Vertical coordinates) for the PLOT TYPE (in the 
graphics data file. Appendix E-3), then the appropriate graphical data are 
also plotted at each data output point. Upon reaching the end of the flight 
profile segment, %PR0P passes execution control back to %SNIT. 

7.1 ROTATED-ELLIPSE PLOTTING SUBROUTINE ('RELIP') 


( When the user selects the "RSBY" plot type, %PR0P draws one or two pictures 

of the payload (depending on how many views are specified) at each data output 
point. These pictures represent orthogonal projections of the payload's 
cylindrical outline into the X-Z and/or the Y-Z planes of the Shuttle's body 
axes. The pictures are composed of straight lines (corresponding to the sides 
of the cylinder) and ellipses and elliptical arcs (corresponding to canted 
views of the circles that represent the ends of the cylinder). 

The ellipses and elliptical arcs are drawn by using the HPL ofs Instruction 
tu move the origin of plotter coordinates to the appropriate point on the 
( HP-9872A plotting surface, and then by calling 'RELIP' to draw an ellipse (or 

portion thereof) about the offset origin. The geometry of the ellipse, refer- 
enced to the offset origin of plotter coordinates, Is shown in Figure 35. 
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pi * a = Semi -major axis of ellipse. ; 

p2 * b = Semi-minor axis of ellipse. j 

3 

p3 ■ Ep ■ Eccentric anomaly of first point on elliptical arc. ^ 

p4 * aE = E L - Ep, where E L = eccentric anomaly of last point on elliptical i 

arc. ; 

p5 = Number of chords to be drawn between the first and last points, for | 

the purpose of approximating the true arc. 
p6 = cos e , where e * ellipse rotation angle, measured from X pL Q T axis 
p7 * sin 0 to semi -major axis. 
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7.1.2 Examples of Usage 


Suppose that It Is necessary to draw the upper half of the ellipse shown 
In Figure 35, centered on the Shuttle body-fixed coordinates X g « 20, Z g « -50. 
We will assume that the relationship between the HP-9872A plot coordinates 
(Xp L0 T» y plot) and the shutt ^ e coordinates (X g , Z g ) has already been defined 
by the execution of an HPL scl_ instruction. We will further assume that the 
lengths of the semi -major and semi -mi nor axes, respectively, reside In regis- 
ters A and B. The values of cose and sine we assume to reside In registers C 
and S, respectively. The necessary Instructions for drawing the desired arc 
are then pen; ofs 20, -50; c // 'RELIP'(A,B,0,7r,N,C.S) t where N is the number of 
chords that are to be used to approximate the true shape of the semi -ell Ipse. 

The 'RELIP' subroutine can be used to draw figures other than elliptical 
arcs. For instance, the instruction pen; c H 'RELIP' (R,R,0,2tt, 40,1 ,0) would 
cause a circle of radius R (approximated by 40 chords of equal length) to be 
drawn about the origin of coordinates. If this were followed by the Instruc- 
tions pen; cil 'RELIP' (R,R,tt/2,27t, 3,1 ,0) an equilateral triangle would then be 
inscribed within the circle. 
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7.1.3 Computations 


The coordinates of any particular point P, on the ellipse shown In Figure 
35, are given by the equations 

XpiQT = a cosE cose - b slnE sine 
and 

Yrlot a a cosE sine + b slnE cose. 

As Indicated by the flow chart In Figure 36, these equations are embedded In a 
loop and evaluated at regular Intervals of E. The elliptical arc is approxi- 
mated by drawing straight lines (chords) between the points thus defined. 

Since a great number of chords may be required to obtain an accurate 
approximation of the true arc, the trigonometric identities 

sin (E + 6E) * sinE cos 6E + cosE sin 6E 

and 

cos (E + 6E) * cosE cos 6E - slnE sin $E 

are used within the loop. This avoids repetitive references to the si£ and cos 
functions, which could result in excessive execution time requirements. 

It should be noted that the * RELIP ' subroutine makes use of the volatile 
simple variables (H,I,J,K and W,X,Y,Z); therefore, values assigned to those 
registers by the calling routine will be lost upon execution of 'RELIP'. 
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Figure 36. 'RELIP' Subroutine Logic Flow 
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7.2 FOURTH-ORDER RUNGE-KUTTA INTEGRATION SUBROUTINE ( ' RK4 * ) 

The function of the ' RK4 ' subroutine Is to propagate the state of the 
Shuttle/payload system across one Integration time step, using the fourth-order 
method of Runge-Kutta. 

7.2.1 Input Data 

7. 2. 1.1 Argument List 

pi ■ h = Value of time step. 

7. 2. 1.2 Others 

The 'RK4' routine calls 'DERIVS'; therefore, all the 'DERIVS' Input data 
listed in Sections 2. 3. 1.2 and 2. 3. 1.3 must have been defined before calling 


{ 

7.2.2 Output Data 

* RK4 ' updates the contents of the extended array of state variables (r25- 
r74), which will be referred to symbolically In this section as [X], It should 

/ 

1 be noted that time Is one of the state variables In the HFRMP, and Is Integrat- 

ed just like any other state variable. It should also be noted that, although 
the contents of the derivatives array ( rl 25-rl 74) will change as a result of 
executing 'RK4', the values that reside there upon return from ' RK4* do not 
represent the true derivatives at the end of the time step. To obtain the 
true derivatives. It Is necessary to call 'DERI VS' again after executing 
' RK4 ' . 

7.2.3 Example of Usage 

{ The Instruction cf£ 'RK4* (300) would cause the state of the system to be 

advanced 300 seconds. 


7.2.4 Computations 


Let [X] represent the array of first derivatives (with respect to time) 
of the state variables In the array [X]. Let [X] p represent the state of the 
system at time t n , and [X] n+ ] represent the state at Vl* where 

h * Vl - ‘n 

Is some relatively small time Interval. Since [X] 8 f([X]), then according to 
the fourth-order method of Runge-Kutta we can write 

[X] n+1 8 [X] n + ([k] 1 + 2[k] 2 + 2[k] 3 + k[4])/6, 

where 

DO] = h f([x] n ) 

[k] 2 8 h f([X] n + >s [k]]) 

[k] 3 8 h f([X] n + Js [k] 2 ) 
and 

[k] 4 8 h f([X] n + [k] 3 ). 

The error Introduced Into the system state by a single fourth-order RK 
Integration step Is on the order of 

i. [x ] hi . 

d t 5 5! 

As a rule of thumb, It has been found In using the HFRMP that the Integration 
stepslze should always satisfy the relationship 

h <_ 40 degrees/ oi max . 
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The symbol w max represents the angular velocity magnitude, measured In degrees 
per second, of whichever one of the four state-variable reference coordinate 
systems (B,b,G, or g) Is rotating most rapidly with respect to Inertial space. 
That Is to say, no reference coordinate system should ever be allowed to rotate 
more than 40 degrees during a single Integration step. In some cases It may be 
necessary to reduce the single-step rotation limit to 20 or even 10 degrees to 
achieve the desired Integration accuracy. 

A flow chart of the 'RK4' subroutine Is shown In Figure 37. The [Y] array 
is stored In registers rl75-r224, and the [Z] array Is stored In registers 
r225-r274. Each quaternion In the [X] array (only) Is normalized, each time 
that array Is updated, by means of the computational sequence 

\/^q + qf + 4 + 4 * m 

V m ** % 

q^/m ■+ q<| 
q 2 /m -*• q 2 
q 3 /m q 3 . 

The quaternion normalization procedure Is never applied to the [Y] and the [Z] 
arrays; to do so would introduce a systematic Integration error. 














7.3 MAIN LOGIC 


The XPROP computations are described by the flow chart contained In Figures 
38 a through 38h. The memory allocation table In Appendix D will have to be 
consulted to correlate the logical symbols appearing In the flow chart with 
the r- register numbers that appear In the HPL code, which Is contained In 
Appendix C.6. 



Figure 38a. XPROP Logic Flow 
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Figure 38b. %PR0P Logic Flow (cont.) 
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Figure 38c. %PR0P Logic Flow (cont.) 
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Figure 38d. %PR0P Logic Flow (cont.) 




































ACA/ANLE STATE OF SHUTTLE /PAYLOAD SYSTEM TO NEXT PRINT/PLOT POINT 



Till* PRlWT/PLOT INTERVAL 


Till ■ t AT END OF FLIGHT PROFILE 
SEGMENT 


TU]«MM ALLOWABLE IHTI&RATlOM 

STEF «IU 

TCH1» INTEGRATION STEP SHE 

TC93* NUMBER OF STEPS REQUIRED TO 

REACH NERT PRINT /PLOT POINT 


TLIT* PRIN1/PLOT COUNTER 


Figure 38h. %PR0P Logic Flow (cont.) 


166 












REFERENCES 


1. S. M. Klndall , "User's Guide for the HP-9825A High Fidelity Relative Notion 

Program (HFRMP) TRW IOC No. 78:2511.4-33, 16 October 1978. 

2. S. W. Wilson, "Equations for Calculating Orblter Surface Erosion and Break- 

age Rates In IUS and SSUS SRM Exhaust Plumes," TRW Report No. 

2841 5-H007-R0-O0 , 6 June 1978. 

3. S. W. Wilson, "A Quaternion Method for the Numerical Integration of Satel- 

lite Orbits," TRW Report No. 99815-HOOl-RO-OO, 30 December 1977. 

4. United States Committee on Extension to the Standard Atmosphere (COESA), 

U.S. Standard Atmosphere, 1962, U. S. Government Printing Office, 
Washington, D. C., December 196?. 

5. M. D. Zuteck, "Improved Shuttle On-Orbit Aerodynamic Moment Coefficient 

Models," TRW IOC No. 77:2511.4-37, 1 August 1977. 

6. S. W. Wilson, "Shuttle On-Orbit Aerodynamic Moment Equations," TRW IOC No. 

77:2511 .1-3C, 28 April 1977. 

7. M. D. Zuteck, "Shuttle On-Orbit Aerodynamic Drag Model Improvements," 

TRW IOC No. 76:2511.4-60, 3 September 1976. 

8. E. C. Lineberry, "I variant Orbital Elements for Use In the Desertion of 

Motion About an Oblate Earth," NASA/ JSC IN 74-FM-84, 4 December 1974. 

9. S. W. Wilson, "Engineering Description of the OMS/RCS/DAP Models Used in the 

High Fidelity Relative Motion Program (HFRMP)," TRW Report No. 

2841 5-H009-R0-00 , 6 October 19/8. 

10. D. M. Henderson, "Euler Angles, Quaternions, and Transformation Matrices - 

Working Relationships," NASA/JSC Memorandum 77-FM-37, July 1977. 

11. Lyndon 8. Johnson Space Center, "Shuttle Operational Data Book, Volume I, 

Shuttle Systems Performance and Constraints Data," NASA/JSC Document No. 
JSC-08934 (Vol. I), Revision A, updated 13 June 1978. 

12. E. W. Purcell and W. B. Cowan, "Relating Geodetic Latitude and Altitude to 

Geocentric Latitude and Radius Vector," ARS Journal , July 1961. 

13. B. G. Jackson, "Update of On-Orbit Plume Impingement Values," JSC Memo- 

randum No. EX32/7901-4, 17 January 1979. 

14. H. F. Fllegel and T. C. Van Flandern, "A Machine Algorithm for Processing 

Calendar Dates," Communications of the ACM , Vol. 11, No. 10, p. 657, 
October 1968. 


1 


167 




J 








A Ralston and H. S. Wllf (editors). Mathematical Models for Digital 
Computers . Vol. 1, Chapter 7, John Wiley and Sons, New York, 1£$2. 

C. A. Denham, "Trajectory Prediction Parameters for AAP, Space Station/ 
Space Shuttle, and Interplanetary Mission," MSC Internal Note No. 
70-FM-33, 5 March 1970. 

Nautical Almanac Offices of the United Kingdom and the United States of 
America, Explanatory Supplement to the Ephemerls , p. 30, Her Majesty's 
Stationery office, London, 1961. 


168 






APPENDIX A: 


BASIC QUATERNION OPERATIONS 


A.1 QUATERNION ALGEBRA 


Reference A-l defines a quaternion to be a mathematical quantity of the 

form 

Q 0 + IQ, + JQ 2 + kQ 3 (A-l) 

A A 

where Qq, Q-j , Q 2 » and Q 3 are real numbers and where the products of 1, j, and 
k are governed by the following conventions: 


lol * joj * kok e -1 , 

(A-2) 

A A A A A 

jok » -ko j » 1 , 

(A-3) 

A A A A A 

kol « -lok = j, 

(A-4) 

and 


A A A A A 

loj « -jo 1 * k. 

(A-5) 

A A 

In addition to having the properties of Imaginary numbers, the quantities 1, j, 

A 

and k also have the properties of unit vectors that are aligned with orthongo- 
nal coordinate axes In a three dimensional space. 

The circular symbol (o) which Is used to denote the quaternion product Is 
adopted from the notation of Reference A-2. This symbol should not be confused 
with the dot (•) that Is used to denote the scalar product of vector algebra, 
which Is governed by the following conventions: 

A A A A A A 

1*1 8 j*j = k*k = 1 , 

(A-6) 

A A A A 

j«k s -k*j » 0, 

(A-7) 

A A A A 

k»i * -1*k * 0, 

(A-8) 


A-2 


and 


A A A A 

1*j « -j«1 - 0. (A-9) 

It Is worthwhile to observe also the similarities and differences between the 
quaternion product and the vector cross product, which Is governed by the 
conventions 

i x i • j x j * k x k * 0, (A-10) 

A A AAA 

j X k » -k X j - 1, (A-ll) 

k x 1 • -? x k - 3. (A-12) 

and 

1 x J ■ -J x 1 ■ k. (A-13) 

A quaternion can be thought of as having a scalar part Qq, and a vector 

part 

« 10, + jQ 2 +kQ 3 . (A-14) 

Sometimes, then, It Is convenient to express Equation (1) In Its equivalent 
form 


7 8 Q 0 + (A-15) 

The sum of two quaternions Is defined by 

M - (Pq+Qq) + <( p i*0i) + 3(P 2 + Q 2 ) + MP 3 +Q 3 ), (A-16) 

and their product by 


A- 3 


(Wo - fill • P 2^2 - W 


( 

PoQ • 

♦ HPoQ, * P,0 0 * p 2«3 ’ W 

{ 

* J(P 0 <l2 ’ e 2% * P 3«1 * P l»3> 

♦ k(P 0 q 3 ♦ P 3 Q 0 ♦ p,Q 2 - PjQ,). (A- 17) 

♦ 

Equation (A-17) results from applying the distributive law of algebra along 
with the conventions defined by Equations (A-2) through (A-5). It should be 
* noted that. In general, quaternion multiplication Is not commutative (l.e., 

f PotST) * Except for the commutative property of multiplication, quaternions 
satisfy all the requirements for the definition of a field. 

; The quaternion product of a scalar and a quaternion ^s commutative, and Is 

given by 

Sotf - $>S « SQq + 1SQ ] + jSQ 2 + kSQ 3 . (A- 18) 

which follows from (A-17) when the scalar Is treated as a quaternion whose 
vector part Is zero. In a similar vein, the product of a quaternion and a 
vector If formed by treating the vector as a quaternion whose scalar part Is 
zero. This results In 

SoV • (-Q,V, - Q 2 V 2 - Q 3 V 3 ) 

* <<Q„ V 1 * Q 2 V 3 - «3 V 2 J 
'♦ 3<V2 * Vi * <W 


+ k(Q 0 v 3 * 0,V 2 - q 2 v,) 


A-4 


(A-19) 


and 


Vo$ * * VgQg " V 3 Q 3 ) 

+ ♦ v 2 Q 3 - v 3 q 2 ) 

+ 3(V 2 Q 0 ♦ v 3 Q] - v 1 Q 3 ) 

+ k(V 3 Q 0 ♦ Vfl 2 • V 2 Q,). (A-20) 

By examining Equation (A-17), It Is seen that the quaternion product can 
be expressed In the form 

Fotf » P 0 Q 0 - F»ff + PqIT + Q q F + F x (A-21 ) 

which leads to a relationship, 

$>F * Fotf - 2 F x (A-22) 

that Is sometimes useful. Equation (A-22) shows that quaternion multiplication 
Is commutative whenever the vector parts of the two quaternions are parallel 
to each other. 

The conjugate of the quaternion $ * Qq ♦ IQj + 3 q 2 ♦ kQ 3 Is defined by 
5 « Q 0 - IQ, - 3 q 2 - kQ 3 . (A-23) 

The norm of a quaternion Is defined by the product 

So? « ?od * Qq * Q* ♦ <)f + Q 3 . (A24) 


A-5 


A. 2 COORDINATE TRANSFORMATION 


Quaternions that have a norm of unity exhibit properties that make them 
very useful for transforming coordinates from one Cartesian system to another 
having a different orientation. Although the distinction Is not made In some 
of the relevant literature, a quaternion of this special class Is more pre- 
cisely known as a versor. In any event, It Is well to remember that when the 
word "quaternion" is used in relation to coordinate transformation or the 
orientation of rigid bodies, almost always It refers to a quaternion whose 
norm is equal to one (1.0). A convention of using lower-case alphabetic 
symbols to designate versors (unit quaternions) has been adopted in this 
report. In other words, use of the symbology 

q = q 0 + ^ + Jq 2 + *“) 3 (a-25) 

implies that 

qjj + q l + q l + q 3 = (A " 26) 

According to one of Euler's theorems, any two Cartesian coordinate systems 
F and G that have a common origin can be brought into coincidence by rotating 
F through some angle a about a single fixed axis. That is to say, the angular 
displacement of any Cartesian system G with respect to another Cartesian 
system F can be described in terms of a rotation about a single fixed axis 
which is usually referred to as the Euler axis. The orientation of the Euler 
axis can be defined by a unit vector 

d = d F = d G = id 1 + jd 2 + kd 3 (A-27) 

with components d-j , d 2 , d^ which not only satisfy the equation 


(A-28) 




d 1 + d 2 + d 3 a 1 • 


but which also have identical values in F and 6 . Assuming that the angle of 

A 

rotation is restricted to lie in the range 0 a < 2tt, the unit vector d can 
be taken to define both the orientation of the Euler axis and the sense 
(direction) of rotation according to the customary right-hand ccrew convention. 

It can be shown (e.g., see Reference A-2 and A-3) that any vector V 
having the form 


V F = 1 V F1 + jV F2 + kV F3 


(A-29) 


in coordinate system F can be transformed to its equivalent form 


V G . iv G , * JV G2 ♦ kv G3 


(A-30) 


in system G by use of the equation 


r G * ?F r- ■' r F 0 a‘fg* 


(*-31 ) 


where 


q FG = cos { x 2 a) + d sin (» a a). 


(A-32) 


It should be noted that in this system of notation the unit vectors i, j, and k 
are not associated with a particular set of reference axes. Instead, they are 
understood to be aligned with the axes of whatever coordinate system is desig- 
nated by the alphabetic subscript attached to the vector symbol. In this 
connection, a vector symbol having no subscript represents an intrinsic physi- 
cal value that exists independently of the system of reference. For instance, 
the symbol V might_represent the inertial velocity of one body with respect _ 
to another, while V p represents just one of many possible quantifications of V: 
namely, that resulting from the projection of V onto the axes of coordinate 
system F. 


A-7 


The double subscript FG in Equation (A-32) identifies the versor as being 


the orientation versor which defines the angular displacement of coordinate 
i-y tern G with respect to system F . 

The Inverse displacement (of F with respect to G) is defined by 

?GF = §fg * cos ^ °) ‘ d sin & “)» (A-33) 

and the inverse transformation of coordinates by 

Fp = q^p o Vq o qgp. (A-34) 

It is easily seen that the equations 

V G = ^GF 0 V F 0 ^GF (A-35) 

and 

Vp = qpg o V G o <Tp G (A-36) 

are equivalent to (A-31 ) and (A-34). 

We now consider a third coordinate system H whose angular displacement 
relative to G is defined by the orientation versor 

q GH = cos (h e) + e sin (> 2 3 ). (A-37) 

Application of the coordinate transformation law yields 


An alternate rotation that produces the same result is defined by 

qp G = - cos (' 2 a) - d sin (* 2 a), which represents a rotation through an angle 

of 2n - a in the opposite (-d) direction. However, there is no need to con- 
sider the alternate rotation in the present discussion. 




" ?6H 0 0 ^GH 

= *GH 0 (*FG 0 Vp o q F6 ) o q GH 

= ^GH 0 ^FG* 0 ^F 0 ^FG 0 ^GH^’ (A-38) 

whence It follows that 


^FH = ^FG 0 ^GH* (A-39) 

The above result can be extended to any number of successive rotations. For 
instance, suppose that the body-fixed frame B of a spacecraft Is displaced from 
an inertial reference frame I by rotating the spacecraft first through a pitch 
angle e about Y-axis, then through a yaw angle ♦ about its body-fixed Z-axIs, 
and finally through a roll angle + about its X-axis. The total displacement is 
defined by the quaternion product 

q IB = [cos (h e) + j sin (H e)] o [cos Os *) + k sin (h *)] 0 


[cos (>s ♦) + 1 sin Os *) , (A-40) 

which, after carrying out the indicated multiplications, reduces to 

^IB s (c e % c $ _ s e V 


+ *< c e S s * + s e s * V 


+ j(S. C. C. + C s, SJ 

0 \p <f> 0 ip 9 


+ k ^ c e s * c , " s e c * V 


(A-41) 


A-9 






where C 0 * cos (h e), S fl » sin (h e), C. ■ cos (% 0 ), etc. 

Hie angles e, 4 represent one of many Euler angle sets that can be 

used to define the orientation of one Cartesian system relative to another. 

* 

Any relative orientation can be described by three Euler angles representing 
successive rotations about particular coordinates axes, taken In a specified 
sequence; the only restriction being that the second axis of rotation must not 
coincide with the first or the third. Reference A -4 contains a useful 
compendium of the relationships between versors , transfer ::'/ .on matrices, and 
all of the possible Euler angle sets. 


£ 

Not to be confused with the single angle of rotation about the Euler axis 
which produces the same result. 


A-10 


A. 3 TIME DERIVATIVE OF AN ORIENTATION VERSOR 


If the coordinate system B Is rotating with angular velocity fi relative 
to system I, the time derivative of q JB Is given by 

q IB = ** q IB 0 n B ( A-42 ) 

or, alternatively, by 

q IB * Oj o q IB . (A-43) 

The consistency of Equations (A-42) and (A-43) with each other can be readily 
verified by substituting the relationship 

= ^IB 0 ®l 0 ^IB (A-44) 

into (A-42), which yields 

QlB = ^ (djg 0 q IB ) o Oj 0 q^ B (A— 45 ) 

One of the major advantages of using a versor to define the orientation 
of a rotating coordinate system lies in the fact that it has a finite deriva- 
tive at every possible orientation (assuming of course that the angular velocity 
H is finite), thereby facilitating numerical integration of the differential 
equations that govern the rotational motion of the system. Such is not the 
case when Euler angles are used. No matter what rotation sequence is chosen, 
it is possible for the derivatives of two of the angles to approach infinity 
in the vicinty of certain critical orientations. 

The problem of infinite derivatives can be avoided by defining the 
orientation with direction cosines (i.e., the elements of a coordinate trans- 
formation matrix); however, this requires the integration of nine real 
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variables as compared to four In the case of a versor. The Integration of 
direction cosines Is further complicated by the necessity of maintaining the 
normality and orthogonality of the transformation matrix, as defined by six 
different equations among the direction cosines. In the case of versor Inte- 
gration, only a single ancillary condition of this nature Is of concern: the 

maintenance or normality as defined by Equation (A-26). 
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B.l SHUTTLE STATION (STRUCTURAL) COORDINATES 

Shuttle Station coordinates are referenced to a Cartesian system that Is 
fixed to the Orblter structure. The origin of the system lies In the plane of 
symmetry at a point 400 Inches below the cargo bay centerline and 576 Inches 
forward of the aft face of the forward bulkhead of the cargo bay. The X-axis 
Is parallel to the cargo bay centerline and positive In the aft direction, the 
Y-axis Is normal to the Orblter's plane of symmetry and positive to starboard, 
and the Z-axIs completes a right handed orthogonal system. The X coordinate Is 
referred to as the station (STA), the Y coordinate Is referred to as the buttock 
line (BL), and the Z coordinate Is referred to as the water line (WL). This 
system is shown In Figure B1. 
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B.2 PAYLOAD STATION (STRUCTURAL) COORDINATES 

The Payload Station reference system Is defined such that Its coordinate 
axes are parallel to the Shuttle Station coordinate axes when the payload is 
stowed In the Orblter payload bay. The origin of the Payload Station coordinate 
system Is located In the center of the front face of the payload cylinder. This 
system Is shown In Figure B2. Coordinates In this system are always measured in 
inches. 



B.3 SHUTTLE BODY COORDINATES 


Shuttle Body coordinates are referenced to a Cartesian system that Is 
fixed relative to the Orblter structure, with Its origin at the Orblter CG. 

The X-axis Is parallel to the cargo bay centerline and positive In the forward 
direction, the Y-axis Is normal to the Orblter's plane of symmetry and positive 
to starboard, and the Z-ax1s completes a right handed orthogonal system. These 
axes are parallel to the Shuttle Station axes; however, the positive directions 
of the X and Z axes are reversed. This system Is shown In Figure B3. 




B.4 PAYLOAD BODY COORDINATES 

Payload Body coordinates are referenced to a Cartesian system that Is fixed 
relative to the payload structure, with Its origin at the payload CG. The X- 
axis Is parallel to the payload longitudinal axis and positive In the Orblter's 
forward dlrctlon when the payload Is stowed In the Orblter's cargo bay. The Y- 
axis is normal to the Orbiter's plane of symmetry when the Payload Is stowed In 
the cargo bay, and the Z-axIs completes a right-handed orthogonal system. This 
system Is shown In Figure B4. 


* 
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B.5 RECTANGULAR LOCAL VERTICAL COORDINATES 


Rectangular Local Vertical coordinates are referenced to a rotating 
Cartesian system centered at the CG of an orbiting vehicle. The direction of 
the Y-axis Is opposite to the orbital angular momentum vector of the vehicle's 
CG with respect to the center of the earth, the Z-axIs points toward the center 
of the earth, and the X-axis completes a right hand orthogonal system. This 
system Is shown In Figure B5. 



B.6 CURVILINEAR LOCAL VERTICAL COORDINATES 


Curvilinear Local Vertical coordinates are referenced to the CG of an 
orbiting vehicle. As shown In Figure B6, the X and Y components of position 
are measured along the surface of an Imaginary earth-centered sphere that 
pa<.oe.' through the CG of the vehicle. The Z component Is measured normal to 
the sphere, positive In the direction of the earth's center. These coordinates 
are essentially Identical with those defined In Section B.5 when the distance 
from the vehicle CG Is small. 
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B.7 MEAN OF 1950.0 COORDINATES 


Mean of 1950.0 coordinates are referenced to a non-rotating earth-centered 
Cartesian system that Is defined by the orientation of the earth's mean equator 
at the beginning of the Bessellan year 1950.0 (l.e., the Julian Date 
2433282.423357). The Z-axIs points In the direction of the earth's angular 
momentum vector, the X-axis Is aligned with the Intersection of the equatorial 
and ecliptic planes (positive In the direction of the sun as seen from the 
earth at the time of the vernal equinox), and the Y-axis completes a right hand 
orothogonal system. This system Is shown In Figure B7. 
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X-Y plane Is Earth's equator of epoch 


X Is directed toward the mean vernal equinox of 
epoch 

Z Is directed along Earth's mean rotational axis 
of epoch and Is positive north 

Y completes a right handed system 

Epoch Is the beginning of Bessellan year 1950 
for the Mean of 1950.0 system or the date of launch 
for the Mean of Launch system 


Figure B7. Mean of 1950.0 and Mean of Launch Date Coordinates 
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B.8 MEAN OF LAUNCH DATE COORDINATES 


The Mean of Launch Date system of reference Is also a non-rotating earth- 
centered Cartesian system defined In exactly the same manner as the Mean of 
1950.0 system, except that the orientation of the equatorial plane Is fixed at 
the date of launch. This system Is shown In Figure B7. 
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C.1 Base Link (#TRAJ) 


0* " NTR0J< 1 223/ 1 0MOV79) " t ft e “RUN" 

It vtb 6 . limit *0 

2t t#t "kRHAT“»0»l 
31 " ANGI " 1 2cf rc<Pl '2#>4 p 2 
4t if p2<012**p2+p 2 
3t r«t p2 

01 "flNG2"*2ffrc<Pl/2if>*p2 
71 if p2>*Ip2~2*+p 2 
8* ret p2 

9* "ATNrtEc+pSIlto *2 
101 "ATH2"*0‘*p 5 
111 Plt2*P3lP2t24p4 
12c if p3*p4«0l«to +6 
13i if p3>p4l«to +3 

141 atn<abi<Pl^p2>>^P3l if p2<0Jc-p3+p 3 
13* *to *2 

16* c/2-atn<p2/«bi<Pl>>*p3 
17* if p1<01 p3-p3-»p 3 
10* r«t p3 

19* “HM.S’ *abi<Pl>-»p2 
20* lnt'»2' , 00>*p3f int<p3''00> , H»4 
21* 100p4*p3-0Op4+'»2-€0p3)/100+p 3 
22* rtt i>3*fn<»l> 

23* “SECS M *ab#«.»1^100-»p2 
24* Int'’p2>+P3l 100f rc<P2>-»P4 
23* 0O<00»3* i nt < p 4 ) ) ♦ 1 00f rc < p4 ) +p 3 
20* r»t p3i«n<Pl> 

27* "D0TP“* rtt rplrp2+r<Pl*l>r(p2M>*r<Pl*2>r<p2*2> 

20* "$XV“i&lrp2-*r»3lPlr<p2*l>*r'.c3*l>lPl r<p2+2>+r<»3+2'l rtt 
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:?* “VADir*i-»4ifto *2 



) 




ITRW Cent. 

Ml WHM 

91< r»IMnl«n>t 

991 r(»m)Mr<»2M)*r(»9«l> 

391 r(»ltl)«Mr(»l«l)*r(i3*t)tr«l 

341 ”C*$r'tr<Pl4|>r<p242>-r<pi42>r<»S4l>4p4 

39i r<»l*2)ri2*f(»lr(rt*2>*»l 

39* rptr<p24l>-r<pl4|>rp24r<p342> 

37* i»3*f<i»3*l>lp4*rp3lr*» 

39* "R0T*tri»2+H!fto +2 

39* “IR0T**-rp2+H 

49* ri»t‘»Xtr<Pl*imfr<i>l<>2>4Z 

41* r<p24| >♦! I r<p2*2>*JI r<p2+3)-»K 

42* <MH4l I- JJ-KK>X*2< < I J4HK>Y4< IK-HJ>Z>4rp3 

43* <MH-Il4jJ-KK>Y42(<JK4HI>24<!J-MK>X>4r(p34l> 

44* <MH-I I- JJ4KK>2+2< < !K4HJ>X*< JK-Hl >Y>4r<p3*2>* rtt 

49* "Q90T~tl*Ht94Mf r»2/2+Xtr<P2*l>/2<»Ylr<»242>'242S«o 49 

49* "9X9"* j4H4Mlpto 43 

47* "9X9C"*l4Ht*l4MI«io 42 

49* "9CX9"*-14HM4M 

49* Mr<p24|>4XIMr<P242>4YIMr<p243)42lrp24M 
99* Hr<Pl4l)+IfHr<Pl42>*j*Hr<Pl43>4Ktr»:4H 
91* HH-!X-JY-k'Z4rp3 
92* MX4lU4J2-kY*r<p34l> 

93* HY4JM4kX-!2*r<p3+2> 

94* H24KH4lY-JX-»r<o343>*r«t 

99* "0231"*l*M*r<»l4|>+Ilr<e!*3)+Ki »*o 42 

99* "92l3"*-|4MI-r<p|43)*H r<Pl4l)*K 

97* rPl 4 Hlr<Pl* 2 > 4 j 

93* MH-!!4JJ-KK«*XI2<JK-M1)4Y 

99* <2< I J4HI >.r<XX4YY>>4r<e2^|> 
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fTRAJ Cont 


60* »0TM2*<>l.H*K>*XI»8TN2»<J-!iH-IO4V 
61* *ftHG2»<X*Y>4rp2l •fmG2»<M<X-Y>>*r<p2*2>l rtt 
621 '0313’trPl*Hlr<PlM>Mtr<Pl*2>«Jlr<Pl*3>*K 
63* 2» ftTN2» <r< I !♦ JJ> » r<H»H<K> >*r<P2*l > 

641 'ATN2*<KtH>*XI’f)TN2*<JtI>+Y 

63* »flNG2» CX*Y)*ri»2l »ANC2» <X-Y)*r<p2*2)l r*i 

661 ’231Q*i l+HI«o *2 

6?i "213Q" l -1*M 

68* rpl/2+Xtcos(X)+Usin<X)+X 

631 r<Pl^l>/2*Ylco*<Y>*JUln<Y>4Y 

70i r<Pl*2>/2*2lco*<2>*KI*in<2>*2 

71* IJK-MXY2*rp2 

72* 1 J2+MXYK*f < p2^2~M> 

73* XJK+MIY2*r<p2*2> 

74* IYK-MXJ2*r<p2+2+M>lr*t 

73* "3130"»r<Pl+l>/2*Ylcc»<O+Jliin<Y>*Y 

76* <rPl+r<Pl+2>>/2+Xlco»<X>*n»ip<X>*X 

77* <rpl-r<pU2>)/2+21co»<2>+KUirt<2>*2 

78* JI*rp2IYK*r<P2M>IY2'»r<p2*2>l JX*r<p2*3>Jr#t 

73* ’lMAT0’*r<p!*3)*nr<Pl*6>+Jlr(pltl)*K 

88* r(pl*7)-»X»r(Pl*2>*Ytr(pl^3)-*2 

81* l*rpl*r<»l*4) + r<Pl*8>^HI If *2 

82* l 4W*rp2f <I-X »M*r*.P2+l >1 < J-Y >\4+r> i>2+2 > ; •> -Z >W-»r<p2*3>* r«t 

83* 2-H+M(U*2rPl+Htif H<lS«to *2 

84* l/2fH-»UI<l-X>M+rp2l l/4M<*r'P2M >1 <.r.*2)\4*rU2*2 >1< J«Y>U*r<p2*3>l r«t 
83* M*2r<Pl*4>*HI If H<l»»t© ♦ 2 

86* l/2fH+WI < J-Y>W*rp2l <K*2>M*r < p 2*1 >1 l/4M+r<p2*2>l < !*;:*ll*r<p20>l r«« 
87* W42r<Pl*8>*H 

88* l/2rH-»MI<K-2)M*rp2*OY)U-»r*.P2^1>l<I*X>M*r«»2*2>l|. 4U>»r<»2*3‘<* »«t 
lift “DER1VS"*23p1-0 
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#TRAJ Cont 


90« “GRAV-tell 'QCXQ» <0+25. 0+29. 0+8?>l cl 1 »QXQ* <0+87. 0+80. 0+129) 

9H 0+rl+r2ll+r3tctl *R0T* <1.0+25.0+133>l-l+r3lcll *ROT* <1.0+129. 1) 

92* r<0433>4Rlr24/RR4Clr<04l33)4Ulr<0+134>4Tlr<04!3S)48 

93 * . 883248 1 < r23/R> 1 2+B I C < 0< 388- 1 >/2- 1 > +fl I CBS+B 

94* BU+r7*8T+r81BS-A+r9 

95* r<0+7?>+Lt r<0+?8)+MI r<0+79)+NI3C'R+C 

96* Cr2r3<N-M>+rl61Crlr3<L-N>+rl?ICrlr2<l1-L>+rl8 

97* rl24r<0+75)+PI if P*8lfor 1*18 to 1518+rIlnext Il<no “ANGflC" 

98 * "ATM" * . 08669342+E I TT+UU+C I -ESfC+DI EC+C 1 D/’ < 1 -C > +D 
99* 28855591 < 1 +< 1 +3C/4>C/2> +B1 BD/R+AI <R-B> < 1 -Dfl/2) +H 
188* .802377/«xP<r<<H-394134>^12665>+17.449>4Q 
181* 7.292U-5+WlR<r<0+35>+TW>+rll-RUW+r2l-r<0+34>+r3 
182* ell * ROT* ( 1*0+87* 1 > 1 if elMJltO "PL AERO" 

103* "AERO"*t*lrl+r3r3+A?f<a+r2r2>+Vlr2/'V+$lrA/V+C 
104* 8+Yll+Xlif C#01r3/VC+Ylrl/VC+X 
105* ’ATN2* <S*C)+B1 * ATN2* <Y.X)+Aiotb$<$)+Sl sin<2B)+E 
106* < . 786+2. 413«bs< Y)tl . 316X 1-S>+1 . 928S+. 646abs<EYY>+D 
107* -1345PDQV+D! cl 1 »SXV» <D. 1 . 10) 

108* .831<sin<3B)-Ecos<2A>>-.036B+I 
109* <.045co*<3fl)-.162Y)CC-.045Xcos<3B>+J 
110* .067sin(2.096B)-.01<Y+Xrabs<E)>E+K 
111* PQVV+AI53218A+6J 104986A+A1 AI +rl31 BJ+r 145 AK+r 151 9to +4 
112* "PLAERO”* rl22obs<rl >+rl23f < r2r2+r3r3)+A 
113* -APQ+Alcl 1 *SXV*<A.1»18> 

114* -rl21/2+rl50+r2+r3lcll ’OftSP* <1. 18. 13> 

115* •ENDERO"*cll »CRSP* <0+84. 18.1) 

116* cli ’ VSUB* <13. 1 * 13) I cl 1 » IROT’ <10.0+87. 18) 

117* 32.174/r<0+39)+G5cll * SXV' <G. 18. 10) 

118* ell ’ROT* <13.0+80. 1 >lcl 1 * VADD» < 16* 1 . 16) 

119* ”ANGAC“*cll ’ROT’ <0+36.0+80. DILrl + r45Mr2+r55Nr3+r6 


«2>n 
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#TRAJ Cont. 

120* ell 'CRSP'<l»4»4)lcll *VSUB'<16»4»1> 

121* rl/L+r 1 I r2/H+r2l r3/N+r3lcl 1 • IROT* < 1.0+80. 0*136) 

122* if f>l*llno “PLTHR" 

123* “THR“* r94+Jl Jnod2+Il (J-D-^+K 
124* 19. 16K*rl411 19. 16I+rl4010+rl39 
125* 6000G<I+K)+Alcl 1 *SXV*<A»91.13> 

126* r95+Jlif J=81for 1=142 to 14710+rIlnext I5<*to +11 

127* for 1=13 to 151 rl+RC 1-12* J)+r Jlnext I 

123* for 1 = 136 to 1381 rl+RC 1-132* Jl+rll next I 

129* for 1 = 142 to 1471RH-135. Jl+rIJnext Il-jto +8 

130* “PLTHR"* l+rl7410+V+rl3*rl4+rl5+rl64lif rll9=0l9to +7 

131* r?4*TJ if T>AC2 31*to +6 

132* 1+V1 T-AC 1 3+T1 if T<=fll 7 Jl^to *3 

133* AC 7 J+X 1 AC 8 3+Y5 s read 3. fit 73.AC83* if T>AC731no +0 

134* <AC 8 3-Y)/<AC 7 1-X)+AC 6 31 Y-XAC63+AC 5) 

135* -< AC 5 3+TAC 6 3)AC 3 )+r 164+Gi -32. 174GAC 4 )/r64+G 
136* for 1=13 to l5JGr(I+103)+rltnext I 
137* “ENDTHR “sell > IROT* < 13.0+87. 13) 

138* for 1=7 to 95 rl + r< l+3) + r< I +6)*r< I +0*123)1 next I 
139* “AMC“* r<0+76)+AJ if A=01*to "ENDBER" 

140* if A=llfor 1=16 to 181 r< I+O+120)+rI Intxt I*9to +9 
141 * 0+rll-r<O+35)+r21 -r<0+131 )/Rr2+r3 

142* ell ’ I ROT’ <0+36 .0+87. 7>»cll * VSUB* < 7. 1 . 4 > * c 1 1 * CRSP’ • 4. 13. 4 > 

143* 2r5+Alif pl=llif V#0*A-(AC6J'rl64+rl64 r64)*A 

144* 'cl 1 *CRSP* <0+133. l»4>lr<0+34>*Slr<0+35>*W 

145* . 0832481 r24(r23/RR>t2<r5r<0+135)+r6r<0+134))-4Sr8. R+A+A 

146* <r(O+130)-2SM)/R+X10+r4l-X+r55A'RM-<S/R+X/M>r3+r6 

147* ell * CRSP’ <1. 7.7)1 cl 1 ’ VADD’ <4. 7.4) 1 cl 1 ’ROT* (4.0+37.4) 

143* for 1 = 16 to 181 r<I+0+120)-r<I-12)*rIlne?.t I 
149* if r-l = l 1 no "PL AM I “ 
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I50» I*? to 910+rIlnext 1 

1311 for K»16 to 181K-12+HI2H+JI J-l+Gt if rK<01G+JI J+l+G 

1 52* abs<rK>+fltabs<RtH» JD+EI if A>EIE+A1 l+r285 

1531 .0l4*Flr99+Ctr<H+92>+DtF+4r<ftD>/E+TICint<T/C>+T 

154: if T>0tT-F+SS0*B*Y+ZI*to +2 

155: C+TlT-F+$tEESS/16IHV>BtB''oibs<RCH»Gl>+ZtZT/S+Y 

156: <A+BVE+X1XT'S+U 

157: for 1*7 to 91 I-6*HtrI+XRCH» J3+ZRCH»G J+rllnext I 

158: for 1=136 to 1381I-132+H1 rI+XRCH> J1+2RCH>G J+rI?r.ext I 

159: for 1=142 to 1471 1-135+H1 rl+WRCH. J3+YRtH>G3+rIlnext I 

160: next Klcll • IROT’ <?i 8?i7)lcll *VADD* <130t?» 130>l*to *5 

16t: "PLANrtfor 1=161 to 1631 rI-rU-143>+rIlnext I 

162: ell ’ROT* <16*t05*l)l-Lrl-»rli -Mr2+r21 -Nr3*r3 

163: ell ’ IROT* <1* 105* l)lfor 1=1 to 31rl/2+Y 

164: «bs<Y)+X!X+Y+r< I +164)1 X-Y+r< 1+167)1 next 1 

165: “EHDDER“tf or 1=7 to 91 r< I+0+123)+rI {next I 

166: r<0+34)*S+r<0+133)tr<0+35>+W{RNW-r9+r<0+134> 

167: <r7-2$W)/R+r<O+135)10+rll-W+r2f r8/RW+r3 
168: ell ’©DOT* <0+25* 1> 0+125) tell ’QDOT* <0+29>0+36»0+129>t ret 
169: "RUN”: ?et "\Rf1ftT"> 169 
170: end 
“’£018 
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C.2 Response Matrix Computation Link (35RMAT) 

0* "*RNAT<0403'24SEP79>“*<no “RUM" 

i* wtb e.nuist #6 

it lit “V.THIT“*0»1 

3* “RUN“*lkd I rad) if fl*lll<9to +5 

4i din A[ 33.BC 323.FC12. 12 3. IC 3.33. JC 8. 12 3 

5* din L( 33.U1 12. 12 3.RC 12. 12 3.TC 7.443 

6i din BSC IS It CSC 12* 4 1* DSC 25 3. JSC 3.4 3.KII 53 

?« din NSC 44. 3 3* P$C 3 3. V$C 80 3. W$( 80 31 sf * 11 

8* "* WARNING* SS "*Bf 

9* " CONTROL IMPOSSIBLE WITH “->M 

10* " V"*J#C 1 31 " P"-»J«23r PZI"4J$C3J 

11* " JETS-4K* 

12* fnt 

13* osen “SCMDID*. 1 1 sread l.C* 

14* asm “SJFTM“» llfor 1*1 to 441sread l*Mf[I31next I 
15* asen "*JFr. 11 sread 1 » TC * 3 

IS* assn “1*”» 11 sread 1.81*31 if BC323*91*to “RMEND” 

17* for 1*1 to 31B[I + 13-»ICI»I31next I 

18* -BC5 3+IC 2*3 3* IC 3*21 

19* -BC 6 3*K 1 *3 3* It 3* 1 3 

20* -BI73*II 1.2 3*112.13 

21* im> 1*1 

22* for T*1 to 3id*P "JjRHAT“* JSt T31 in* T 

23* files #JSV.*V.*VRleto *3 

24* files #JSP.*P**PR>*PF$sto +2 

23* files fJSPZI»*PZI.»PZIR.*PZIF 

26* ina U.R.Fl" M *V#-.W$lsreod l.JC*3 

27* for C*1 to 12lfor N*1 to 8 

28* JIN.C3-.Mlif M*05 sto +13 

29* for 1*1 to 61 T[ I »M3*rItnexl I 
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30t for 1*4 to 61<rI-BCH-43>/l2+rIlnext I 
311 -r4+r4l-r6+r6lcll *CRSP*(4il»4> 

32* for 1*4 to 6»rI+Tt?iMlr<I-3>+rIfn*xt I 
33* for 1*1 to 61UI I >C 3+rI-MJt I.C JJnext I 
34* M*CM3+P» 

35* if P#C1»13*“R“ 19+11 *to +3 
36 * if P«1.13*"L“l8+Il*to +2 
37* 7+1 

38* if 061 1+3+1 

39* 3. 1071 -Ml if M>381 ,0923+W 

40* UEIfCl+H+UCItCl 

41* next N 

42* 32.174/Bt 13+fi*for 1*1 to 31AUt l*C3+UtI»C31next I 
43* for 1*4 to 61 Ut I *C 3+Lt 1-3 31 next I 
44* not IL+A 

45* for 1*4 to 61 AC I -3 3+Ut I » C 3! next I 
46* if T*l! if C<75 9to +4 

47* int<<C+l> '2>+P54P-2C-l+Slif SUtP»C3>05*to +3 
48* Bf&Cft C 3&Df&J$C T 3&K$+V#Jwrt 6. 1 * V*J " "+V* 

49* if C>6!B$$c‘'ATT”&D$$.J$[ T 3&Kf+U$ 

50* next CS**rt 2>U[ *3>Wtlara U+R 
51* if W## ,,, 'l9to ••REND" 

52* for C*7 to 12* int < (C+l >/2>+P 
53* for H*1 to 10 
54* for H*0 to l5 4+*<P+H>nod3+E 
55* -RC Ei C 3+X1 2E+K? i f X,0}K-1+K 

56* X/Ut E»K 3+Ylf'or 1*1 to 121 Rt I*C 3+YUt I *K3+Rt I • C 31 next I 
57* next H18+X 

58* for H=0 to 1 1 4+<P+H)nod3+El X+Rt E* C 3t2*X 
59* next HlX'Rt P*C 3I2+X1 if X<le-8!*to +2 


XRMAT Cont 


60* next Kioto *2 

61* 4P-2C-HSII# SRCP.C3>01*to +2 

62* B$i"ATT H *'D*U*m&l<**M*lurt 6.1iM*taro U*Rloto “REND* 

63* for H*8 to 1 1 4+<P+H>nod3+EI 6*RC E»C 31 next H 

64* next C 

63* for C=1 to 6 

66* for EM to 6 

67* -RtE»C3*XI2E*KJif X>0JK-1*K 

63* XRCE>K3*Ylfor 1=1 to 12IRC I.CJfYRCI*K3+R[ I»C3Jnext I 
63* 3+RCE* CIS next E 
73* if T=ll*to +3 

71* ir»t<<C+l)/2)-»P54P-2C-l+SI if SRCP»C3>0Jeto +2 
72* BffROT COMP *«<C«C3tD*8<J*me<K#+V$J wrt 6 . 1 *V$ 

73* next C 

74* ”REND”*SPrt 3»RC*3>W#Jif T=l»*to "HEXTT“ 

73* ora R+F 

76* if W*# ,,M leto “FEND” 

77* if V##““*9to “ABORTF" 

78* for C*1 to 6* int < <0+1 V2)*P 
79* for N*1 to 10 
80* for H=0 to lU + <P+H)nod3*E 
81* -Ft E*C 3*X* 2E*K» l f X>0<K>1*K 

82* X/R[E*K3*Ylfor 1 = 1 to 125 Ft I . C 3+YRtI . 1 3*FC 1 . C 35 next 1 
83* next HJ0*X 

84* for H=t to i;i+<P*H>nod3*EJX+FC£»C3t2+X 
S3* next HJX/FC P*C 1124X1 if X'le-Slato +2 
86* next Nieto "ABORTF" 

87* 4P-2C-l*8*if SFCP»C3< a 0Jeto ‘ABORTF” 

88* for H=0 to 1 5 1 ♦ ( P+H)rtod3^E 5 0*Ft E» C 31 next H 
89* next C 
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90* for C»7 to 12 

91* for E*1 to 3 

92* -FtE»Cl*XI2E*Klif X>0IK-1*K 

93* X/FCE.Kl+YJfor 1 = 1 to 12JFI I#C]+YFU»K]*F[ I>C]lnext I 
94* 0*FIE»C3»next E 

95* int<<C+l>/2>+PJ4P-2C-l+S»if SFCP*Cl<=01*to “ABORTF" 
96* next Cleto “FEND“ 

97* “ABORTF" * B$i"FULL COMP’tDftJfCTHcKf-HI* 

98* wrt 6.1*W$lara R*F 
99* "FEND"*SPrt 4iFC*]»M$ 

100* “NEXTT“* next T 
101* 9+BC321 

102* asen “l*".15sf>rt liBC*] 

103* “RMEND”* eet "%TN1 T"i169 

104* end 

*24706 
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C.3 Trajectory Initialization Link (%TNIT) 

A. 

•* *3TNIT( 1230/16OCT79) " * eto “RUN" 

ll wtb 6. 111 list K 

2i nt "5iS8VU“.0»l 

3« "JD“* int<<p2-l4)/l2>->P4 

41 P3-32073+ i nt < 1 46 1 <p1+4800+p4>/'4>*p5 

51 p5+int<3670»2-2-12p4>/12>+i>5 

6* ret p5-int<3int<<Pl+490O+p4>/100)/4> 

7* “PLOTYP"* 

8t if PM“ “fret 0 
9i if PM"N0NE“J ret 0 

10t if P$*“RS8Y" I ret 1 

v - 

11* if PM“CPLV“lret 2 
12* ret -9 
13* “UNIT"* 

14* if PM“ FT"! ret 1 
15* if PM“ KFT" J ret 1O00 
16* if PM" NMI"? ret 6076.115 
17* if PM" M’ Jret 1/.3048 

18* if PM" KM" I ret 1 000/. 3048 


19* 

ret -9 




20* 

"MONTH"! 

[ 



21* 

if 

PM" 

JAN' 

1 ret 

1 

22* 

if 

PM" 

FEB‘ 

1 ret 

2 

23* 

if 

PM" 

MAR’ 

1 ret 

3 

24** 

if 

PM" 

APR’ 

1 ret 

4 

25* 

if 

PM" 

MAY’ 

I ret 

5 

26* 

if 

PM" 

JUN’ 

J ret 

6 

27* 

if 

”0 

II 

JUl* 

1 ret 

7 

28* 

if 

PM" 

AUG’ 

1 ret 

8 

29* 

if 

PM" 

SEP* 

1 ret 

9 
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%TNIT Cont 


301 if P«»" OCT" tret 10 
31* if Pff NOV" tret 11 
32* if P$»" DEC" I ret 12 
33* ret -9 

34* "C"* ret p3<p2-1>+p1-1 

35* "T"*«oix<p1>p2Xp 3I ret p3<p3-l>/2*nin<pl»p2>-l 

36* "DIflG“*0-»M*X 

37* for J*1 to Pllfor l*i to J 

38* » T» < 1 1 J>*Kt P2+K4HI P3+K+KI rH*rK*Y 

39* if I a JIWm*Wll*r<p4+»C»<J»J.pl>>l9to +2 

40* X+2YY*XI0*r<P4*»C'<IiJtPimr<p4+*C»<J»I»Pl>> 

41* next Iinext J 

42* if p0<5ir(M+X)/pl/le6'»P5 

43* rX/pl+p6J0+p7*p8 

44* for J=2 to Pllfor I»1 to J-l 

45* p 3+*T»<I. JXZlif abs<rZX=P6l*to *12 

46* 1->p8Ip3+»T»<I.I>*HJp3+»T»<J.J>*W 

47* rH-rW+XI2rZ+YJ Y/f <XX+YY>-»KI if X<0l-K+K 

48* K/r < 2 < 1 +f < 1 -KK ) ) ) +Y I r < 1 -YY> +XI rH*K 

49* KXX+2XYr Z+ YY rW* r H I KYY-2XY rZ+XX rW+ rW 1 0* r Z 

58* for K»1 to Pllif K a I I eto +4 

51* if K a JI «to +3 

52* pS+’T* <K*I>-*H5p3+’T , <K*J)+W 

53* rH+ZlXZ+YrW*rHIXrW-YZ-»rW 

54* p4+»C’<K*I»p1)->HIp4+ , C'<K.J»p1>->W 

55* rH*ZIXZ+YrW*rHIXrW-YZ-»rW 

56* next K 

57* next Iinext Jlif p8°1 I 0+p 8I eto *3 
58* if P6<p5lret 
59* p6/p!+p6 
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XTNIT Cont. 


€01 P?M*P?lif p7< 1001 «to -t€ 

€li fnt “DIAG ERROR -"i«10.2»2/lurt 6»Plp6lfnt trM 

€21 "RUtTtftb Still if fhlllito +3 

63* din RtdltBC 321tG( 121tR( l2t 12 It T( 73 

64* din fl*C9.21iB«401iC«2t41tD#C3i43tE«3t41iF«4»41 

63* din Gf C 3» 4 It Htt 8* 3 It ISC €» 3 It Jfl 6i i It K*C 60 1* L« 4 It H#t € 1* Ntt € ] 

66* din 0#I6.31iP*C41tQ$C6i33»R«6»31*$$£6t31tUfC 103»M»£801»2f£8011sf* 11 

67* 1 11 “PL"+A*£2 31 ”TP"*fl#£31l "DX"+fl*£ 4 31 "DY^fift 5 31 ”DZ"->fl#C63 

68* “DR“-»R$C71I “DS"->fl#C8 JI"DT“ , >fltC9l 

69* "OI150"*C4£ 131” IMLD"*C$t 2 1 

7B* w M50“+D4£ 1 31 " SLV"*D«23I" PLV"-»D«33 

71 * "RSBY“*E« 1 II "RSLV"-»Ef[ 2 It “C$LV“*E« 3 II "ESLV”-»E« 4 It "$M50"*E« 3 1 

72* "RPBY"+F4C 1 11 “RPLV"+F« 2 II “CPLV->F« 3 II “EPLV"+FtC 4 3 

73* “SENS“*G« 1 II “PC0N"*G$f 2 31 "RIf1P"*G« 3 3 

74* "SNfT+H*! 1 11 ”ECC"->H«2 3I "INC“-»H$£33I “RflN”-»H« 4 3 

75* ,, ARG*+Hf£5 3l "TRA”*HIE6 3I "HA "♦H*£73I"HP ”+Hf[8 3 

76* “PCH” + I4t 1 II ”YAW"-tI$£211 "ROL ,,, »I$t3 3l "RXB”*It£ 4 3rRYB”*I#£5 3l “RZB’ -tltEei 

77* ”X“+J$E 1 II “Y"+J#E21I "2"*Jtt33l "R"*J#C411 "S"*J#C53I "T""»J4t63 

78* *>X “-» 0 « 13 r+Y " *08(211 "+Z "*0$£33I"-X “♦0«43t"-Y "•>0«53l“-2 “*0*E63 

79* ”RFX"*ClSE 1 31 “RLX" *Q$£ 2 31 "RRX"*Q8t 331 "RAX” +Q8E 4 31 "RX ” *Q8£ 5 31 "0ML"*Qt£ 6 1 

80* ”RFR"*R$C 1 31 ”RLR”+R« 2 31 “RRR“-»R$f 3 31 "RAR"*R$£ 4 31 ”RR “-»R« 531 "0MR”*Rt£ 6 ) 

8t * ”RFT“*S$C 1 31 ' RIT ”*S*E 2 31 ”RRT“*S#£ 3 31 “RAT“*S*I 4 31 “RT 331 “OMT”*St£ 6 3 

82* o**n "#ID"t ll*r««d 1»Z#I "#TRAJ”*Z8£ 1.53 

83* beeolem ’ENTER TRAJECTORY ID TEXT*'»K#I if flfl3l«to +8 

84* 80-Wn<K*>*Kl if K«0l*io +2 

83* for H»1 to KIKtfc" “-tKflntxt H 

86* fmt 2 ■ / »c80»2/*c80#/lwrt 6tZf*Kflfnt 

87* Ikd lrad*0-»r2851if/180-»Q 

88* 2O925722*r23»1.4076469«16*r24 

89* oun "l*" 1 1* stead l.Bf * 31 BC 1 3*r39 
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*TNIT Cont 


90* Bt 2 3*r 1 1 Bt 3 3*r3l BC 4 3*»r6t -Bl 5 ]*rSI -Bt € 3*r4l -BC 7 3*r2 

91* ell *DIAG’(3iltli7>lrl*r77lr3*r78lr6+r79lcU * IMflTQ* <?»8B> 

92* <10?6.?-Bt81>/12*r84lBt91/i2*»-85l<3?5-BtlB3>/l2->r86 
93* < 1826 . ?- 8 m 3>xl2*r273IBC12 3/12->r2?6l<375-B£13)>/l2-»r277 
94* for 1=1 to 3IQBC l+13 3-»rllnext Hell »231Q* <1(2?8> 

95* Btl?)/l*3+r99lfor 1-40 to 47l0*rlln*xt ! 

96* o**n ”2*“» ltsroad l(Bt*3JBt 1 3*r64 

97* Bt 2 3*r II Bt 3 l*r3l Bt 4 3->r6i -Bt 5 3*r5l -Bt 6 )+r4l -Bt 7 3*r2 

98* ell *018G'(3» 1» 1(7)1 r Hr 1021 r3+rl03l r6+r 1041 ell * MATO* <7. 105) 

99* -Bt83/12-»rl09IBl93/12-»rll0l~BC103/12*rlll 

1 00 * -Bt 11 V 1 2 ♦ r 282 1 8(1 2 3/ 1 2* r 283 1 -Bt 1 3 )/ 1 2* r 284 

101* Bt 14 J/24-»R+rl20IBt 153s*2+L+rl21 

102* «RR+rl22»2Rl+rl23 

103* for 1=65 to 70l0+rlln»xt I 

104* os«n "3*“»Hfor 1=1 to lll*r«ad lfGtlllntxt ! 

105* ft» (Gt 1 3>+P#l * PLOTYP' +Gt 1 Jl if Gt 1 >01 *to +13 
108* ft* <Gt2)HP*»»UNIT»-»«23 

107* b< client "READY TO PLOT rilltf fl*13=0lno +0 
108* PClrfwrt 785( " i*540( 878# 14?78( 18838" I *>tn# 2lfxd 0 
109* tel 8( 14(0(9lp*nlolt 0(-.3l)bl Kf 
118* GC31+CIGI 9 3-CGl 83*YI if Gt53»0l*to +4 
111* id 0* - 1 4C»8* -9CI of * -Gt 8 1( -Gt 9 )l Gt 8 J-CGt 5 3*X 
112* xo* 0.GtlCO.0.Gt8J.Gtllllx«x 8(-Gt 1 0 3* 0 * 1 ©* « 1 1 3 
113* >ox 8(Gt 183(0(Gt93(Gt 11 31 vox 8i -Gt 18 3.0( YiGt 1 1 3 
114* if Gt 4 3=01 ft© +4 

115* jel 8* -14C (8» -9Clof $ -CGt 5 l-Gt 7 3» -Gt 9 31 Gt 7 3-CGt 4 3+X 
118* ;;ox O(Gtl03(0*GC7 3-C 1 8 * Cl 1 1 31 xox 0. -Gt 103(0iX(GC 1 1 3 
117* .ox 0(Gt 103(8(GC93(Gt 11 31 vox 0.-GI 183(0(Y.Gt 11 3»*>«n4 1 
118* oi*n "4*“#llir«ad 1(BC*3 
119* BMl-Ylfti * Bt 2 3>+Ptl ‘MONTH' +I1IBC 33+D 
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XTNIT Cont 


1201 (• JO* (Yi H« D) -2433282V36524 » 224T 

1211 <2384.948*<.302*.018T>T>T*X 

1221 <2884. 2S6-<. 426*. 042T>T>T*Y 

1231 <2304. 948*<1. 093*. 018T>T)T+2 

124t «/648000*Cl «/2-CX*rl ICY*r21 -«/2-C2*r3 

I25t ell * 3130* <1* 19) 

1201 'SECS* <8C53>+r741 8C63*rl24 

1271 «stn "?*"»llon «nd lt”G02*lsr«ed 1»BC*3 

1281 "G02"*8t5 J*r751BC 183*rl0010+r76*r94+r95+rl01+rll9 

1291 «s«n ”5#”»lUr«ad l,B[*)lft* <8tt3>*P* 

1301 6078. 1 15*11 r24*GI if P»»”OH50"IT8t23*fllBC33*EI*t© *2 
131 1 TBC 2 3*r23*Rl TBt 3 3+r23*PI <«-PV<ft*P>*EI <«*P>/2*A 
1321 8<l-EE)-»PIQBl7J-»FIP/<l*Eco*<F>>'»R»r<C/P>E*in<F)4$ 

1331 QBt 5 !♦ r 1 1 QBt 4 3* r2l QBt 6 3+F*r3 

1341 If Pf«“OM50"lcl 1 *313Q*<l»25>tcll *QCXQ»<19»25»25>lfto *8 

1351 rl*HI r2*ll r3*Ulco»<I>*CI»in<I>*8 

1361 co»<U>*»:i*in<U>*UKKM.l*M2KL*M 

1371 1.62405*-3*JIJr23T2*VlV/P*B 

1 381 R*BDDM/6*RI fi+8< l -3®®/2>/3*AI r<G'A>/A*M 

1391 S-BWD»H/3’»Slft+2flflV<l-3ODLL>/3RRR*0 

1401 U-B<l-7D®/6)N/2P+r3ll+VC0M/2RR*r2IH+BCN/4P^rl 

1411 ell *3130* <i*25> 

1421 rf/2*r41-r4*r5ie+r6»cll '313Q* <4.0>lcll 'OXQ* <25*0*25) 

1431 r<<2/R-l/fl>G-S$>''R*»‘35!$+r341R-*r33 
144t for l»l to 3IOBt I*83*rllnoxt I 
1451 ell '2310' < !*87)lcll *GXQ» <25»8?»29> 

1461 for I»36 to 3810Bt l-233*rfln«xt S 
1471 cl l 'DER1VS' • 0>l*'i80*G 

1481 ft* • 8t 12 if P#»" SlVIcll *R0T*<l*87»4)lcll ' VADO* <4* 36» 36) 
149* on* “6*“*lltr#od 1*81*3 
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XTNIT Cont 


ISBifor !•? to 12fB(l-9)-»r!ln«xt I 
1911 ft* <Btn>4P*fif Pf«”RSLV"f MO M 
192* ell *R0T»<1.8?.4>lcll *V$UB* 06.4.4) 

1931 Cll *CRSP , <4.7»4>lcll ’ VADD’ < 10.4, 10 ) 

1541 ell * 1 ROT* <7.87. 7)tcl 1 * 1R0T* <10.87* 1B> 

1551 r9-r33*r9l r!2-r34*rl2 

1561 Cll *CRSP*U.7.4>lcll * VADD* < 10*4. 10) 

1971 cll * I ROT* <7.25* 7)tcll * 1R0T* < 10.25. I©> 

158* f *DOTP* <7. 7>*R*r58lcl 1 »CRSP*<7.10 4> 

159« r*D0TP*<4.4)*H»H/RR*W*r60 

!<A» for 1»4 to 61 -rl/'Horll-r<I*3>''R*r<I*3>lnext I 

I6H -*DOTP»<10.?>*.591cll *CRSP»<4.7.1> 

I62t cll * IMATQ* < 1 .30) 

1631 for 1*1 to 310BC I ♦Sl+r II n*xt Hell *2310* <1*1 12> 

1641 ft* <Bm)*P#lif Pf«" PLV-lno *2 

1651 cll *0CXO* <50. 29*8) tell *GXCI*<0.112.112> 

166* cll *QXQ* <58. 112*54) 

1671 for 1*61 to 631GBC I-483*rlln«xt 1 
168* ft* <8( 12)>*Pfl >f P#*“ HSe'tMO *6 
169* if P#*~ PlV”l*»o *3 
1701 cll *0CXQ* , 29*?4*0>»ell ’ROT* < 36.0.4) 

1711 Cll ♦VAIi&‘ 4.$l.*p»fto *3 
172* cll *I*EP.IVS*-P 

1731 cll ’ROT* <■ 1 * 112* 4 *1 cl 1 'VW> 4.61*61 > 

1741 0 *«n “7#“»1 

175* if GC 1 ]«llchoin “*/.$SVU“. 169 
176« choir. ’.SHIT-. 169 
177* «n<J 
*23100 
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€.4 OrMter Geometry Link (XSSVU) 

0* "k88VU<U08/168EP?9> a i Vio "RUN" 

It Mb 6* III list M 
2i eet ”*SNIT"»0il 

3t N SYH”iobs<int<p8>>‘»HlP5/H+Jtcos<J>*It sln<J)+J 

4t eos<P4>*Xtsin<P4>*YI0+K 

9t Pit p14Xp3*p2*Yp 3I if K-HI rtt 

01 JX*IY*MIIX-JY*XIN*YIK*l*KI«o -l 

71 "RUM"«d*«l 12G(2)+MDG(3)+C 

8t Ptni 21 if GC9)*0t*to N 8VU N 

9t tel 0.-l4C.0.9C»of* -DGC81.DGC91 

Set lin DG( 8 3* DGt 8 J-CGt 5 )» -DGt 9 1* CGC 6 J-DGC 9 ) 

111 oft -12r*85t l2r80-3?5 
12i -!♦« 


131 

PtnlPlt 9*201 

141 

Pit 

$50* 202 

191 

Pit 

S1O8'203 

101 

Pit 

$408* 380 

171 

Pit 

$409*312 

191 

Pit 

$298*327 

19i 

Pit 

$130*330 

29i 

Pit 

$125*403 

2H 

Pit 

$150*399 

22» 

Pit 

$175*339 

231 

pH 

$200,403 

24i 

pH 

$220,412 

251 

pit 

$241,424 

201 

Pit 

$232*424 

27» 

Pit 

$220,422 

28i 

Pit 

$200*413 

291 

pit 

$175*409 
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XSSVU Cont 


301 Pit 8150*400 
311 Pit 8125*413 
32s pit 8120*415 
331 Pit SI 30*430 
341 Pit $139*450 
35s Pit S143»462 
36s pit S1454^5 
37s pit S144* 490 
38s Pit S140*509 
39s pit S130*515 
40s Pit S115*524 
41s pit S100*526 
42* Pit S80»526 
43* Pit S60*524 
44* Pit $40*518 
45* Pit S30*512 
46* Pit S28*513 
47* Pit S4*816 
48* if S<0Jl*Sl9to -35 
49* Pit -S4*816 
50* line 2* llS15*X5 484.8*Yi240Z 
51* PenlPlt X.YHlt X+$Zcos< 52> * Y+Zsin<52> 

52* PenlPlt X*YlPlt X+$Zcos<l 19.5)* Y+ZsinU 19.5)11 ine 
53* "SVU"* if Gt43=0J 9to "END" 

54* scl 0» 14C* 0* 90S of s CGt 5 3+DGt 7 3* DGt 9 1 

55* li« -DGt 7 3 * CGC 4 3 -DGt 7 3 * -DGt 9 3 * CGC 6 J-DGC 9 3 

56* of s 12r84-1076.?*12r36-3?5 

57* PenlPlt 1613*291 

58* Pit 1534*298 

59* Pit 1526*336 
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%SSVU Cont 


601 

Pit 

15501335 

6lt 

Pit 

1575* 334 

62t 

Pit 

1600*334 

63 i 

Pit 

1623*336 

64t 

Pit 

1609*412 

651 

Pit 

1600*410 

66< 

Pit 

1575*402 

6?* 

Pit 

1550*392 

68 « 

Pit 

1525*382 

69: 

Pit 

1517.376 

70t 

Pit 

1502.441 

71* 

Pit 

1525.441 

72i 

Pit 

1550.442 

73: 

Pit 

1575.445 

741 

Pit 

1598.453 

75i 

Pit 

1578.526 

76* 

Pit 

1550.517 

771 

Pit 

1525.506 

781 

Pit 

1500.489 

79: 

Pit 

1492.483 

88: 

pH 

1482.519 

811 

pH 

1576.562 

82: 

pH 

1700.816 

83: 

Pit 

1593-816 

84: 

Pit 

1325.548 

85: 

Pit 

1317.540 

86: 

Pit 

1311.533 

87: 

Pit 

1308.524 

88: 

Pit 

1307.516 

89: 

Pit 

1 !07.420 
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%SSVU Cont 


90 » Pit 9?0»420 
911 Pit 970(900 
92* pH 909(900 
93* Pit 491(499 
94* Pit 480(493 
99* Pit 472(489 
90* Pit 434(490 
97* Pit 300(394 
98* Pit 283(389 
99* Pit 270(378 
100* Pit 299(372 
101* Pit 250(364 

102* ell * SYM* (262.5(338(27.5(135(102(10) 

103* Pit 263(308 

104* Pit 286(300 

105* Pit 312(294 

106* Pit 337(289 

107* Pit 362(286 

188* Pit 387(283 

109* Pit 414(280 

110* Pit 450(278 

111* Pit 492(275 

112* Pit 535(273 

113* Pit 575(272 

114* Pit 625(271 

115* Pit 675(269 

116* Pit 725(268 

117* Pit 775(267 

118* Pit ?00«265 


119* Pit 1000(264 



*SSVU Cont 


120! 

Pit 

110Oi263 

121! 

Pit 

1150*262 

1221 

Pit 

1250*261 

123! 

Pit 

1300*261 

124! 

Pit 

1350*264 

1251 

Pit 

1400*267 

126! 

Pit 

1450*271 

127! 

Pit 

1500*275 

128! 

Pit 

1550*281 

129! 

Pit 

1600*288 

1301 

Pit 

1613*290 


131 ! line 2. t» 545. 2-»Xl 484. 8*Y; 2402 

132* PenlPlt X*YJpH X+Zcos<100>> Y*2sin< 100) 

1331 PenJplt X.YJplt X*Zcos<38>. Y+Zsin<38> 

1341 560+X5480+Y 

135! PenJplt XiYJplt X*Zcos<30>* Y+Zsin<30) 

136« PeniPlt X>YJplt X+395cos<-32>. Y+395sin(-32>U i 

137* “END"s pen 1 pen# l 

1381 chain ”*SNIT”*169 

139t end 

*22773 


C.5 Flight Segment Initialization Link (XSNIT) 

01 "*SNIT<i340/27OCT79>"t*to “RUH* 

It wtb 6'UIUst »6 
2t ut "5{PROP**0» 1 
3t "KftH M s 

4t if P*»“ 0”} ret 0 

3$ if Pf=“ IRH"J ret I 
6 i if P$=”LVRH"I ret 2 
?s ret -9 
8t "JSEL"$ 

9s if Pfa- V“l ret 0 


10S 

if 

Pfa" 

P* 

'1 ret 

1 

1 IS 

if 

Pf*" 

pzr 

'» ret 

2 

12s 

res 

. -9 




13s 

“COMP-s 




14S 

if 

Pfa” 

• 

‘5 ret 

0 

15s 

if 

P$=" 

NONE’ 

'5 ret 

0 

18 s 

if 

Pfa” 

ROT' 

' J ret 

1 

17s 

if 

Pfa” 

FULL’ 

‘1 ret 

2 

18s 

ret 

-9 




19s 

"KTOMS" 

5 



29 s 

if 

Pfa" 

• 

' 5 ret 

0 

21$ 

if 

Pf=" 

HONE' 

'* ret 

0 

22s 

if 

Pfa" 

L' 

'* ret 

1 

23$ 

if 

Pfa" 

R' 

'1 ret 

2 

245 

if 

Pfa” 

L+R' 

1 ret 

3 

235 

ret 

-9 




26: 

"KTRCS” 

5 



275 

if 

Pfa" 

■ 

’ i ret 

0 

285 

if 

Pfa" 

NONE' 

' J ret 

0 

295 

if 

Pf a " 

+X' 

’ J ret 

1 
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mmssssk 


r 


O 


i ( 




43 

44 

45 
48 
4? 

48 

49 
58 

51 

52 

53 

54 

55 
58 
5? 

58 

59 


%SNIT Cont. 


set 

it 

PM* 

-X* 

tr«t 

2 

31« 

it 

Pl»" 

♦Y* 

lr*t 

3 

321 

it 

Pt«“ 

-V 

‘Iret 

4 

33 * 

it 

p#a« 

♦2* 

'1 ret 

5 

341 

it 

p#a“ 

-2‘ 

1 ret 

8 

35: 

it 

P#=" 

+R0L' 

‘fret 

7 

38: 

if 

r*= M 

-ROL' 

‘ 1 ret 

8 

3?: 

it 

p*=“ 

+PCH* 

‘1 ret 

9 

38: 

it 

Pfa“ 

-PCH* 

‘1 ret 

10 

39: 

it 

P*=" 

♦YAW 

*1 ret 

11 

40: 

it 

pfa" 

-YAW ' 

• * ret 

12 

41: 

ret -9 




42: 

“DflTYP” 

: jmp 

el 



*B*r 

"■»U#lret 

3 

■»B$1 “HHMM.SS 

"Hit! ret 

1 

*B$1 "HHMM.SS 

“->Uflret 

1 

■»B41 "HHMM.SS 

ret 

1 

+B*r 

"♦Ut Iret 

1 

♦B*r 

“+U#I ret 

3 

■»B*mEG'$EC 

"+Utl ret 

1 

+B#5“DEG"SEC 

"*Uf 1 ret 

1 

*B*rHEG - SEC 

"*UII ret 

1 

+B*P 

"■»Ut 5 ret 

3 

■»Bf 1 "DEG 

“-►Utlret 

1 

♦BfPDEG 

“■MJtlret 

1 

^Bfl'HEG 

“ Hit! ret 

1 

■*B#1 " 

"•Hltlret 

3 

•>B*!“ 

“•►Uf! ret 

3 

+Btl" 

” *U1 J ret 

3 

*6f» ” 

ret 

3 
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» 




m a i l 


%SNIT Cont. 


60) "PL AERO FACTOR "+Btl" 

611 "PL RATE OPT (INCR* IRiLVR) "♦8*1" 

62t "PL X8 RATE OR INCR "♦B*l "BEG/SEC 

63t "PL YB RATE OR INCR “♦B*l "DEG/SEC 

64) "PL ZB RATE OR INCR "♦B*! "DEG/SEC 

63) "PL ATT NAINT OPT <D»IRH*LVRH> "♦B*l" 

66) "PL THR TABLE (NONE* IL*$A»$D). . . "-»BfI" 

6?) "SLIST") 


63) tnx l»2/»"FLT PROFILE SEGMENT” »f 3*0 

69) tn t 2*2't "ITEN*»2x» “DESCRIPTION* »38x» "VALUE" *2/ 

70) fnt 3*f4.0*c42*f I4.4*cil 

71) fnt 4*f4.0*c42*f 14.0> cl 1 

72) fnt 5*f4.0*c42*cl4*cll 

73) fnt 9*2/lff>t 

74) if Pl>llwtb 6*11 
73< wrt 6.1*Pllwrt 6.2 

76) for 1=1 to p25JnP ’DATYP* (!) 

77) wrt 6.3* I * B$*BC 1 1»U*I eto +3 
73* wrt 6.4* I *6liBC I l*U$l?to +2 

79) fts <BC I 3>*Pf Jwrt 6.5* I *B#*Pf *(J# 

80) next Ilwrt 6.91 ret 

81) "RUN") rod?«/180*Q 

82) TCS3M+T[6110*GC 12 It if Tt63*ll 1*GI 1211 "”+L# 

83) 24*Llon end l. 'FIN" 

34) for 1 = 1 to Llftread l»BC!llnext Hell ’SLIST’ <TC6l*L> 

85) for 1 = 1 to 31 ’SECS* (Bt I + l 1>*TC I 31 next 1 1 Tm+r74”T[ 1 3 

86) Bt53*r?S 

87) ft » (BC103>*P*1’KAM’+Ker76 

88) ft* (Bt 14 3>+P#J * JSEL* 

89) if J=0i"*V"+mi9to +3 


"♦Utlret 1 
"♦U*l r*t 3 
"♦USIrtt l 
“♦Utlret 1 
"♦Utlret 1 
“♦Utlret 3 
"♦Utlret 3 
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XSNIT Cont 


90* if j»ir*p"+n*!«to *z 

9i* "#P2r«>M$ 

92» fts <BU51>*P$I'C0MP»*C 

93t if K>0Snax<l»C)*C 

94* if J=0?nin<l»C)+C 

95* if C»05M*i. M R"*N#lsto +3 

96* if CMIHttlT+Ht+Ntlito +2 

97* 

98* ts*n M*»2lsre«d 2>R[*J*W$ 

99* if Mf # ""If «t />c80»/5urt 6iW*5fnt 

100* fts <BU7])*P#5*KTRCS**r93 

101* fts <Bt 163)*P$5 , KTOMS , *r945 if r94=05*to +6 

102* <1518-1076. 7>/12-»rl?0*r25<492-3?5>/12-*r3 

103* if r94=1588/12*r2 

104* if r94=2J-88/12*r2 

105* ell ’VflDD* <l>84»91)5f*D0TP* (91*9l)+X 

106* for 1=91 to 935 rI/X*rI5next I 

107* for 1=96 to 98 * GBC I -85 3* r U next I 

108* ell 'DERIVS'^Jir/lSB+Q 

109* 0+r7*r85 -r33*r95cl 1 »CRSP* <1»7. 10) 

110* rl2-r34+rl2! cl 1 ’ROT’ <10.87» 10) 

111* fts <B£6])*P# 

112* if P#=" INCR" 5 f or 1=4 to 650*rI5next 15-9*0 +3 
113* if Pt=" IR"5for 1=4 to 65 r< I+32)*rI5next I5?to +2 
114* ell ’POT* <1 » 87» 4>5cl 1 ’ VSUB* < 36i 4» 4) 

113* for 1=4 to 65 GBt I+33-rI + rHnext I 
116* oi*n Hf * 25 sread 2»R[*1.W# 

117* if Mf*“"5fnt /.c80*/5wrt 6»Wf 5 fht 
118* for K=4 to 65 if rK=05 jto +6 
119* 2K+C5 if r\ 35C-1+C 
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XSNIT Cont 


D« 



I 



♦ 





120« n*x<0«rK/R[KiC]>+T 
12U for 1-10 to 121 rl+TRC I-9»C)+rltn*xt 1 
1221 for I»36 to 381 rI*TRC I-32.C3+rI Inext I 
1231 for 1=42 to 471 rl+TRI I~35»C]+rllnoxt I 
124* noxt K 

125* ell *IR0T’(?»25'7>lcU * IROT* <18»29> 10) 

128* ell »CRSP»<7* 10.4)1 r33*R 
127* r»D0TP'<4i4)-»HlH/RR4M-»r35 

128* for 1=4 to 61-rI«'H*rIl-r<I+3>/R-»r<I*3>lnext I 
129* -*DOTP»<10»7Hr34!cll ’CR$P» <4.7. 1> 

138* ell * IMATQ’<1»25> 

131* as*r» Mt»21sreod 2.RC*3 

132* 6C 183+rl001fts <Bt 23 3>->P#l * KAM'+rlOl 

133* n* <BC 24 3>+P# 

134* if P$="N0NE"5 no +4 
135* if P#=” "1 no +3 

138* if Lf="“lno +4 

137* fot "SRM ” » c4 # ” IGNITION COMMAND IGNORED" » /I urt 6*PfIfp»t 
138* if rl 19401 f nt "SRM ".c4.“ BURN CONTI HUES"..' Iwrt 6.L#lfwt 
139* no +8 
140* Pf+L# 

141* for 1 = 1 to 45 if P$t I » I ]='' " 1 next I 
142* “00”t.Pll I.4 3*N4iN$[3»6 3*P# 

143* os ?n Nf»31&read 3» AC 2 ]♦ At 3 3* AC 4 ]* AC 8 3 

144* r74+AC 1 31 r?4+R[ 2 3+At 2 35 At 8 3’At 5 31 8*At 8 3*At 7 3 

145* l-»rll9-»rll610-»rll7->rll8 

148* os*n "8*"»51 iPrt S.Pt.KS 

147* ell * DERI VS* ''1>1 tf/180+Q 

148* 0*r7+r81-r58*r91ell 'CRSP* < 1 » 7. 18) 

149* rl2-r59*rl2!cll * ROT* <10. 1 12* 10> 
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XSNIT Cont 


iset rts <bci9]>+p* 

1511 if Pt*"INCR' , Ifor 1*4 to 6IB*rIln*xt lino *3 
1521 if PM" IR"lfor 1*4 to 61 r<l*5?>->rllnext lino +2 
1531 ell *R0T*(l»112»4>lcll *VSUB'<61.4i4> 

1541 for 1*4 to 6IQBC I*16]-rl+rlln#xt I 
155« for 1*61 to 631 rl+r<l-57>*rlln*xt I 
156* ell »R0T»<4. 105.1) 

157* for 1*1 to 31 rlr(l+101)*rlln*xt I 
158* ell *IROT»<1»105>1) 

159* for 1*1 to 3JrI/2*Ylob»<Y>-»X 
160* r<I+64)+Xmr<I+64>lr<I+6?>+X-Y*r<I*67>lntxt I 
161* chain "V.PROP”. 169 
162* "FIN"*if rl 19*01 fto +3 

163* fnt "SRM ".C4." STILL BURNING AT TRAJECTORY TERMINATION". 2/»wrt 6.L* 
164* tort 5* "end" 

165* fnt 2/.-:SO»2/»c80#2/lwrt 6.2f.Kt 
166* for 1*1 to 3lbeeplw«it 5O0ln*xt I 
167* if L##""*9et "4PIDI" 

168* *nt "START NEW RUN ?".I*if fU13-0lno *2 
169* ««t "ftDBED" 

1 70* *tP 
171* end 
* 2'?819 
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C.6 State Propagation Link (%PR0P) 





t 


t 

l 


j 


ei - kPROP < 1 1 03/ 1 60CT 79 > " » «t © "RUN* 

It wtb 6*111 list §6 
2l ««t "#PIDr»0»l 

3* "REUP"*«bs<int<p5>>*HJp4/H-»Jlco»<J><MI»in<J>*J 

41 cos<p3>*Xtsin<P3>*YI0+K 

5* Pit Xp1p6-Yp2p7*Xp1p7*Yp2p 6I if K*Hfr*t 

6t JX+IY+MI IX-JY*XlW+YIKM+Kloto -I 

7* , 'RK4‘ , *p14p2Ip2/2'*p3I P2/6+P4 1 1 *p3 

8* for 1*23 to 741rI+r<I+150>‘»r(l+200>ln«xt I 

9* ell ’DERIVS* (0)fcl 1 »DERIVS*<1> 

10t for 1*25 to 74lr<l+150>+P3r<l+100)*rlln?xt I 

11 1 for K*0 to 25 by 251 for J*25 to 29 by 41 J*K-»HI0«*M 

12* for 1*0 to 3iW+r<H>H>t2*Mln#xt llfM-HI 

13* for 1*0 to 31 r< I+H>/W+r<I*H>lnoxt lln*xt Jlnext K 

14* if p5*4lr«t 

15* for 1*225 to 2741 rI+p4r<I-10O>*rIln#xt 1 
16* if p5*11p2/3->p41 sto +3 
17* if p5»2lP2*p3lfto *2 

18* p2''6-»p31 f or 1*175 to 2241 r< l+50>*rlln*xt I 

19* p5+1*p51«o -10 

20* "RUN"* rodlO+Tt 71 

21* fnt 1*2 '* "GET *"*fl0.3*“ HHnM.S$“*/ 

22* f «t 2* e2* “ * “ * c4* ♦ 9. I*c4*f8.1*c4*f8.2*c4*t8.2*c4«t8.2*c4«t8.2*c4 
23* fnt 3* c2* " * “ * c4» f 9. 1 * c4* f 8.5*c4*f 8. 2* c4* ♦ 8. 2* c4* t 8. 2* c4* f 8. 2*c4 
24* frit 4*c2*“*“ *c4tf9.2*c4*f8.2*t4»f8.2*c4*f8.3*c4*f8.3*c4»f8.3*c4 
25* frit 5*c2* "* "*c4# < 10.2* c2* f 9.2*c2*f 9.2*e2*f 10.2*c3*fl0.2*c3*fl8.2*c3 
26* fiit 6*c2» ”* "»c4* f 10.0*e2»f9.0*c2*f 9.0*c2*t 10.2*c3*f 10.2* c3*f 10.2* c3 
27* fnt 7*c2* “* “ *c4* f 10.0*e2*f 9.2»c2*f 9.2*c2»f 10.2*c3*f 10. 3*c3*f 10. 3*c3 
28* fnt 8*c2*c3*f9.0*c4ff8.0*c4*f8.0*c4ff8.0*c4*f8.O*c4*f8.0*c4 
29* fnt 
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XPROP Cont 


301 "PRINT"* if TC7J«ed2-0lwtb 8*11 
91* wrt 6.1.'HN.S*<r?4> 

92* ell * DERI VS’ <0)1 l80/«+Q 

33* if r 285001 f«t "•** RC8 CAN’T CONTROL 88 RTT"#/lwrt 6lf*t |0*r28S 

34* for I-t to 3lrl*r<I-tl74>ln*xt I 

35* for 1-7 to 0fr!*r<!M?8>ln«xt I 

38* "SSTRRN"* r24*GI r33*RI r34*Sl r3S*M 

37* GRr <2G-RRRMN-R88 > *ft I RR< 2/R- 1 /R-88/G > *P 

38*’ 8f<P/C>*YIP/R-t*Xlr<XXW>'»EI*ATNt»<Y»X>*F 

30* #/2+r4i -r4-»r5f 0+ r8l cl 1 *313Q’<4'7)tcll *QX0C* <25»7» II) 

40* ell »0XQ*<19»ll.7>lcll *03l3»<7i4> 

4t* Q»ANGt*<r4>*NIQr5-HIQ»flNGl*<r0-F><U 
42* 8078. 115+TIQF+F 

43* wrt 0. 3» Af tl)* C# 1 1 1. ft ' 1 , HOC I )• E t Hf 1 2 h I ♦ Hf I 3 h Nt Hit 4 )♦ J» H*t 5 1* F# H*C 0 J 

44* 1.82405t-3*Jl Jr23t2-»VIV/P-*B 

45* ell *0313* <ll»4)lc©*<r5>*CI*in<r5>*D 

40* co»<r0>*Ktsin<r0>+LIKK-LL+HI2KL*N 

47* A-2ARV< 1 -3DDLL > / 3RRR+AI f <G/A > /A*H 

48* A-B< 1 -3DD/2 > /3+fl I R-BDDM/0+RI S+BUDDN/3+8 

40* RR<2/R-l/«-SS/G)’»PI$r<P/G>*YIF/R-l4X 

50* r<XXm>*EI*ATNl’<Y.X>*F 

51 * r8*B< 1 -7DD- 8>M/2P*»UI r5-VCDM/2RR* I ^ r4-BCN/4P*H 
52* GI*llG*ANGl* <H)-*NIO*flNGl* <U-F>*JIQF+F 
53* <A<l-E>-r23VT*PI<A<l*E)-r23>/T*A 

54* wrt €.2tfttm»CiC21#B»HtC7J#P.H*C81t I.H*m»NtH#C4J# J.Hfl5)»F.H«0) 
55* wrt 8 

58* “SSR0T“*dl *0X0* < 10* 20(7) I cl 1 *Q23t*<7»4> 

57* Or4*nOr5‘*J«C»*ANGl*<r6>H'lOr30*XIQr3?*Y»Qr38'*Z 

58* wrt 8. 4* Aft 1 3*Dft 1 ]» It If( 1 )*J« IK 2 )«Kf IK 31* X* If( 4 )»Yt lit 5 )«Zt lit 8) 
59* ell *R0T* <1*87* t > I c 1 1 *VSUB» <38» 1* 7)lel 1 *Q231*<8?.4> 


XPROP Cont. 


68* Qr4+ I IQr3*JIQ» UNCI * < r6>4KI0r74XI QrBoYI OrB* 

6H vrt 6.4*At( 1 )»89(2)»I* Iff 1 )*J*I9t2)#K*l8t3)»X»I4(4)»Y*I$tS)»Z*t8t6] 
62* wrt 6 

63< M PLROT*iell ’DERIV8* <1 >1180/00 
641 for !■! to 3lrl4r<l4l77>fn*xt I 
651 for 1*7 to 9frl4r<l4l81)fntxt I 
661 oil ’OXO* <19*54* 7) I oil ’Q231’<7*4> 

671 Qr44U0r54,J| O' RNGl • < r6>4KI Qr6l 4X1 0r62+Yf Qr63+Z 

68t wrt 6«4*AIC2)*Dft 1 )# t.Ift 1 )* J* tf(2)*K*!8(3)*X*18(4)*Y*It(S)*Z»I4C6) 

69* Oil ’R0T*<1»112* l>lell »VSUB*<61*t*7>fell 'Q231» (112.4) 

70* Or44llOr54jlO»ANGl»<r6>4KIOr74XIQr84YI8r94Z 

71* wrt 6*4*AS(2)*Df(3)«I»I4( 1 )*J*I*(2)*K*!*(3)*X»!$(4)»Y*It(5)*Z»ISt6) 
72* wrt 6 

73* ell *OCXO' (58*25*8)1 if 0(1)421 oto 49 

74* if G( 12)*0I14C( 12)1 oto 48 

75* "CPIV" * GC 2 )C( 3 )♦$! GC 2 )G( 9 )4T 

76* ell »O2l3*<0*4>l-r4r584r7lr5r584r8lr58-r334r9 

77* for L«8 to If If GC5-U»8l»to 44 

78* sol e*-14$*e*-9$IGC2)GCe-U+X4Mtlf U»1IM+SG(5)4M 

79* of* -M. -Til in X* X-SGC 5-L )» T • T-SGt 6 ) 

88* of s r(8-L>.r9»*>lt B.eietnlcll *REIIP* <$/28. S/28.8.2o 38. 1*0> 

81 < o«nlr»*xt L 

82* 84r7*r8f-r33*r9lell *CRSP' <175.7. 18>I rl2-r34*rl2 
83* 04rl34rl4l-r58+rl5fcll *CR$P» <178* 13* 16>l rl8-r59*rl8 
841 oil *R0T'<13*8*l3>lcll *AOT» < 16.8. 16> 

85* Oil ' VSU6*» 13*7* 13) I ell * VSUB* <16* 18* 16>l if rll9«8ltt© 412 
86* “PL8M58“*cll ’ROT* <188*8* 188>l cl 1 ' V8U8' < 1 88* 1 85* !89> 

87* ell *IR0T»<ll«.54.1>lcll ' 1R0T* < 1* 19. 1 > 

88* sort 5*r?4~ACl)»rS64*rl*r2*r3 

89* ell ' I POT* < 13*25* 1>I ell »!R0T»fl*19*l> 
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XPROP Cont 


Ml ell 'IftOT'(14»25*4>tdt • IROT* <4* 19*4) 

91* i»rl 9trl»r2»r3»r4»r3»r4 

Mi wrt 4.6.A*C21.E*t5J.rl.J*Ul.r2»J«2).r3.Jtt31.r4.ft*t41*rS.AM51.r6.fl*C41 
931 ell * IROT* <I89«29» Did 1 ' IROT* < I. I9» l> 

941 t»n 9»rl»r2tr3 

95i if r?4<AC21l «to +2 

94t %pr% 9»"*ncri0*rll9 

97i “PLTRAN"»f or !«! to 31 r<l*l2>*rlln«xt I 

981 cl! *CRSP*<i?9»tt4>lcll *V$UB*<l6i4»4> 

99t wrt 4*9»Af(2]»Ef(2)»rl«J*( 1 ]• r2t J#( 2 It r3«Jf(3)»r4 »Af( 4 It r9*Af(9Ji r4* ATI 4) 
1881 dl *IR0T* <34t87t7)ldl *V$U8» (175.7.7) 
tell dl 'CR$P*<?#lt7>lcll 'VADD* <4t7t4> 

1821 6ll *R0T'<Ii8?tl>ldl *ROT*<4t87.4> 

183t wrt 6.4»A*C21*Eft 1 lirlt JfC 1 Itr2t JtC2I#r3t JfC 3)tr4tA*C4)»r5tA$C9)ti*4'AfC6] 
1841 wrt 4 

1891 ell •8CXe»(94.29»131>in CC 1 If II ft o M3 
1841 If CC 121*81 net 12 31 fto ♦ 14 

1871 ‘ , R88Y N iC(2]CC3]<»8ICC2)C(9)<»Tlrl28*RI>rl2nrl8l8«rlI*rl2 
1881 dl • ROT* <189* 191 1 7)1 ell *R0T*<18tl81tl8> 

1891 for L«8 to II if CC3-t.]«8lfto *19 

1181 sel 8»*14S#8»”9SICC21CC0*L J*X*MI if 1»1IW*SGC5)*W 

lilt of I -Mt-TIlin X* X-SGt 5-1 1# T# T-SGt 4 ) 

1121 of* r<2-L>* r3l»lt OtOlefnlcll *PEUP* <S 28t S '20*« • 2*2«t 3* l*9>l(-«n 

113* rt 1 1-L>*UI rl2*VI -r(8~L>*XI -r9*Y 

1 14* if <l-2L»r«l8*L> 8IX*U*XIY*V*YI-U*U1-V*V 

115* f<U0*VV)*Elftr.-,«/<9.1-«E rl21>t2)'»BI»f E>81 1*CI0»0I no *Z 

114* V-E+Ci-U/E*D 

U7t of* X.YIcll , REUP*<ft.B»8.2«*38.C*8>l»« £>81*10 *2 
1 18* of* U*V(cll 'PEL IP* <fl*B>8»«*15*C*B>l»olt -U#-V 
1191 Mft*no*t L 
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%PROP Cont 


1201 “TPSENS**cIl * VSUB* (282* 109» ?> I cl 1 ‘ROP <7. 181»7> 

121 > ell • ROT’ <61 i 181 » 13) tell •V8UB»<13»36»13)tcll •CRSP'<13.7»10> 

1221 ell *VADD* <4*10i4)lcll »VA0D»<li?»l> 

123s ell ‘VADD't 1*84*1) tell *VSUB* (1.273.1) 

124i ell , R0T , <l*278*l)fell »ROT» <4.278»4> 

1251 rl-»XI r2*YI r3*Zl r4+Al r5*Bt r6*C 
1261 XX+YY*Qir<Q+ZZ>*RlrQ-»<) 

1271 *ATN1* <Y*X)-»S1 ’ ATN2' <Z*Q)+T 
128* <AX*BY+CZ)/'R+U 

129* if Q/R<le-650*Vlr<AA+BB+CC-UU>/R»M»9to +2 
130* <BX-AY)/'QQ-»VKCR-UZ)/'QR-*M 
131* 180/«+QIQS+S*QWIQV-»V»CIW+W 

132* wrt 6.7* All 3 1. G«C 1 3* R. J$C 4 ). S. J« 5 3.T. J« 6 1. U. Aft 7 1* V* Aft 8 3. U> AM 9 3 
133* wrt 6 

134* , 'RIPC0N'‘*for I®1 to 61 r< l+64>*rltnext I 

135* wrt 6. 8. Aft 2 3* G« 3 3* r 1 . 0#C 1 3. r2. 0« 2 3. r3. Oft 3 3. r4. Oft 4 3. r5. Of t 5 3* r6* Oft 6 3 
136* wrt 6 

137* for 1*0 to 75r<I+40>*rIfnext 1 
138* r3+r4*AJA+r2*B 
129* r6+r7+ClC+r5*D 
140* r2+r5*Il r3+r6*Jl r4+r7+K 
141* A*C*XlB+D-tYl r0+rl*Z 

142* wrt 6.8* Aft 1 3*Gft 2 j* r2*0f[ 1 3. r3.0ft 2 3. r4*Oft 33. A. Oft 4 3. B> Qft 5 3. r0.QtI 63 
143* wrt 6.8* Aft 1 3*6112 3. rS.Rft 1 3*r6.Rft23. r7.Rtt3 3»C.Rft4 3.B,Rft5 3. rl.Rft63 
144* wrt 6.8* Aft 1 3* G#C 2 3. 1 .Sft 1 3. J. Sft 2 3* K. Sft 3 3. X. Sft 4 3* Y. Sft 5 3*2. Sft 63 
145* wrt 6 

146* if obs<r74-Tn 3 X. 06011 chair. "5{$MIT" . 169 

147* “ABVANS" * Tt 2 3(1 -f rc< r74/Tt 2 3> ) *D! if ob*< D>< . 0081 ITt 2 3*D 

148* D*r74-»E1 if 7U 3-E< . ©001 f TC 1 3*EIE-r74*D 

149* irt .99999D/TI 3 3) + l-»Tt 5 3ID/Tt 53*Tt 4 3 
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i <ii mmm 



XPROP Cont. 

150* if Tt 5 1*11 if abs(TC 1 J-Tt 4 l-r?4>< . 0001 1 Tt 1 J-r74*Tt 4 1 
151* fxd 2ld»P "V.PR0P T * *VHN.S»<r?4>lcll ‘RK4»<TC41> 
152* Tt 5 l-l *11 5 II if Tt 5 3*01 Tt 7 J*l*Tt ? II dSP "“leto "PRINT* 
153* 9to -3 
154* end 
*290? 





D-l 






























REG. SYMBOL UNITS 



DESCRIPTION/NOTES 


SS local -vertical system 


orientation versor, with respect 


to MLD ECI system. 



SS body orientation versor. 


with respect to MLD ECI system. 



SS CG distance from center of earth. 


SS Keplerlan angular velocity magnitude. 


SS body angular velocity with 


respect to Inertial space. 


resolved onto SS body axes. 


SS gross weight (constant) 


OMS propellant consumption 



aft left 


translation 


44 ! W PRRX 



aft right control 


forward 


aft left 


aft right | control 


(RESERVED) 


(RESERVED) 



RCS propellant 


consumption 
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REG. | SYMBOL | UNITS 



PL STATE VARIABLES 


DESCRIPTION/NOTES 


1 PL local vertical system 


orientation versor, with 


respect to MLD ECI system 




PL body orientation versor. 


with respect to MLD ECI 


system. 



PL CG distance from center of earth. 


/Ji 2+ llb ft se 


PL Keplerlan angular veolclty magnitude. 


PL body angular velocity with 


respect to Inertial space, 


resolved onto PL body axes. < 


PL gross weight 


Integrals of control torque applied 


* In the positive direction about each 


/t 3+ lb ft se<| j body axis. 



69 /A 1- 


71 I II 


1b ft sed 



Integrals of the control torgue magnitude 


applied In the negative direction about 


each body axis. 


(RESERVED) 


(RESERVED) 


(RESERVED) 


r 
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REG. SYMBOL 
75 Ka 

~TZ 


I SS CONSTANTS 

UNITS 1 DESCRIPTION/NOTES 

SS aero factor. 

SS att (Mint flag (0:D, 1:IRH, 2:LVRH) 


SS principal moments of Inertia 




SS principal axes of Inertia 


1 

orientation versor. with respect 


to SS body axes. 



Position of SS actual CG with 


respect to nominal CG, In SS 


body axes system. 


o q TD (not constant 



value computed In 'DERIVS') 



Ug (QMS thrust unit vector. In SS body axes system. 


QMS thrust flag (Ornone, 1:L, 2:R, 3:L+R) 


RCS thrust command flag (0,1,2,**** 


SS attitude deadbands 



sec | SS DAP cycle time 


TV •*. 









































































































































































































































































































































E.l SHUTTLE DATA FILE 


SHUTTLE DATA FILE t 


I TEH DESCRIPTION 


VALUE 


1 HEIGHT 

2 IXX 

3 IYV 1 

4 1 22 

5 IY2. 

6 IZX 

7 IXY 

8 CG STA 

9 CG BL 

18 CG WL 

11 SENSOR STA 

12 SENSOR BL 

13 SENSOR WL 

14 SENSOR AXES PITCH 

15 SENSOR AXES YAW.. 

IS SENSOR AXES ROLL. 
17 DAP CYCLE 


208817.0888 LB 
887382.8888 $LUG-FTt2 
6386877.8888 SLUG-FTt2 
6694367.0088 $LUG-FTt2 
-971.0000 SLUG-FTt2 
247376.O008 SLUC-FTt2 
5622.0080 SLUG-FTt2 
1095.3800 IH 
0.3008 IN 
377.4080 IN 
1095.3008 IN 
8.3008 IN 
377.4O00 IN 
0.0000 DEG 
0.0000 DEG 
0.0000 DEG 
80.0000 MILLISEC 


'«i»;"NAL PAGE *8 
Of POOR QUALITY 
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E.2 PAYLOAD DATA FILE 


PAYLOAD DATA FILE * 


ITEM DESCRIPTION 


VALUE 


1 WEIGHT 

2 I XX 

3 IYY 

4 IZZ 

5 IYZ 

6 IZX..... 

7 IXY 

8 CG STA 

9 CG BL 

10 CG ML 

11 TARGETPOINT STA 

12 TARGETPOINT BL. 

13 TARGETPOINT ML. 

14 DIAMETER 

15 LENGTH 


372S3.00O8 LB 
7439.0068 SLUG-FTT2 
45468.0060 SL«JG-FTt2 
45433.0600 SLUG-FTt2 
0.0000 SLUG-FTt2 
0.8600 SLUG-FT12 
0.0000 $LUG-FTt2 
300.6360 IN 
0.6080 IN 
0.2600 IN 
360.6300 IN 
0.0000 IN 
0.2060 IN 
112.0660 IN 
428.6600 IN 
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E.3 GRAPHICS DATA FILE 


GRAPHICS DATA FILE $ 


ITEM DESCRIPTION VALUE 


1 PLOT TYPE < NONE » RSBY » CPL V > RSBY 

2 PLOT UNIT. PU <FT»KFT » NKI . M. KM> FT 

3 PLOT SCALE 120 PU/IN 

4 X PANEL WIDTH 14.0060 IN 

5 Y PANEL NIDTH 0.0000 IN 

6 PANEL HEIGHT 9.0000 IN 

7 X MAX 15O0.0000 PU 

8 Y MAX 0.0000 PU 

9 Z MAX 00.0000 PU 

10 TIC INTERVAL 20 PU 

11 LABEL INTERVAL 5 TICS 



E.4 TIME/ATM DATA FILE 


TIME'ATM DATA FILE ♦ 

ITEM DESCRIPTION VALUE 


1 LAUNCH YEAR 1931 

2 LAUNCH MONTH ( JAN' FEB* MAR* FEB 

3 LAUNCH DAY 2? 

4 LAUNCH GMT 1935.0000 HHMM.SS 

5 INITIAL GET 1005.0030 HHMM.SS 

6 ATMOSPHERIC DENSITY FACTOR I.0O00 


E-5 






E.5 SHUTTLE ISTATE DATA FILE 


SHUTTLE ISTATE * 


ITEM DESCRIPTION 


VALUE 


1 STATE TYPE <ON50» IMLD) 

2 SNA 

3 ECC 

4 INC 

5 RAN 

6 ARC 

7 TRA 

8 ATT REF <SLV> 

9 PITCH 

10 YAH 

11 ROLL 

12 RATE REF <H50» SLV> . . . . 

13 XB RATE 

14 YB RATE 

15 ZB RATE 


ON50 

3598.1726 NNI 
0.0012 
28.3014 DEG 
-80.5443 DEG 
-29.4152 DEG 
3.3863 DEG 
SLV 

-100.9080 DEG 
-33.9000 DEG 
-168.0000 DEG 
N50 

0.0000 DEG 'SEC 
0.0000 DEG 'SEC 
0.0000 DEG/SEC 


E.6 PAYLOAD I STATE DATA FILE 



PRYLOAD I STATE $ 



ITEM DESCRIPTION 


VALUE 



1 STATE TYPE <RSLViRSBY> 

2 X 

3 Y .. 

4 Z 

9 XDOT 

6 YDOT 

7 ZDOT 

8 ATT REF <PLV»SBY> 

9 PITCH 

10 YAH 

11 ROLL 

12 RATE REF <M50. PLV»SBY> 

13 XB RATE 

14 YB RATE 

19 ZB RATE 



RSBY 

6.3779 FT 
-O.0250 FT 
-8.1292 FT 
0.2120 FT/SEC 
0.0000 FT/SEC 
-0.3392 FT/SEC 
SBY 

58.0000 DEG 
0.0000 DEG 
0.0000 DEG 
1150 

0.0000 DEG/SEC 
0.0000 DEG/SEC 
0.0000 DEG' SEC 



APPENDIX F 


USER-SUPPLIED FLIGHT PROFILE DEFINITION 
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F.1 FLIGHT PROFILE SEGMENT 


FIT PROFILE SEGMENT 20 


ITEM DESCRIPTION VALUE 


1 FLT SEGMENT TYPE <PROP> PROP 

2 FLT SEGMENT LENGTH 2.4000 HHMM.S8 

3 PRINT/PLOT INTERVAL 0.0500 HHMM.SS 

4 MAX INTEG STEPSIZE 10.0000 HHMM.SS 

5 SS AERO FACTOR 1.0000 

6 SS RATE OPT < INCR. IR.LVR) IR 

7 SS XB RATE OR IHCR 0.0000 DEG'SEC 

8 SS YB RATE OR INCR 0.0000 DEG'SEC 

9 SS ZB RATE OR INCR 0.0000 DEG'SEC 

10 SS ATT MAINT OPT (D. IRH.LVRH) IRH 

11 SS XB DEADBAND 2.0000 DEG 

12 SS YB DEADBAND 2.0000 DEG 

13 SS ZB DEADBAND 2.0000 DEG 

14 SS RCS OPT (V.P.PZI) P 

15 SS XC COMPENSATION (NONE* ROT* FULL) ROT 

16 SS OMS CMD (NONE# L*R*L+R> 

17 SS RCS CMD (NONE* +X* -X* +Y* ... * ♦YAM. -YAM) 

18 PL AERO FACTOR 1.0000 

19 PL RATE OPT ( INCR* IR.LVR) IR 

20 PL XB RATE OR INCR 0.0000 DEG'SEC 

21 PL YB RATE OR INCR 0.0000 DEG' SEC 

22 PL ZB RATE OR INCR 0.0000 DEG/SEC 

23 PL ATT MAINT OPT <D. IRH.LVRH) IRH 

24 PL THR TABLE (NONE* IL*SA*$D> IL 


APPENDIX G 



DIGITAL OUTPUT DATA 


Vr ” * 






J V f fL*W .* w i . ■ 


6.1 DIGITAL OUTPUT DATA 


GET • 1056.350 HHNN.88 


S0ION50 
SSI INLD 

3E17.3 SNR 0.00407 ECC 
107. S HR 149.0 HP 

20.31 INC 
29.46 INC 

279.22 RAN 
279.57 RAN 

26.00 RRG 
19.99 RRG 

153.42 TRR 
159.50 TRR 

SS« M50 
SSI SLV 

•21.13 PCH 
-114.70 PCH 

-57.40 YAM 
0.00 YRH 

. 90. 13 ROL 
190.00 ROL 

-0.090 RX8 
0.000 RXB 

0.000 RY0 
-0.065 RY0 

0.000 R2t 
0.000 R20 

PU N50 
PLI PLV 

101.70 PCH 
-3.51 PCH 

-12.72 YRM 
-10.14 YRM 

87.96 ROL 
0.93 ROL 

0.009 RX0 
-0.011 RX0 

9.000 RY0 
0.068 RYB 

0.000 R2B 
0.013 R20 

PLISH50 

PLiRSLV 

PHRSBY 

-325994 X 
413200 X 
32454 X 

-277410 Y -196164 2 -869.89 OX 

-992 Y 225774 2 1273.34 DX 

1011 Y 469739 2 101.56 DX 

-716.98 DY 
-155.27 DY 
155.31 DY 

-347.56 D2 
106.44 DZ 
1164.21 DZ 

TP l SENS 

47O800 R 

1.70 S 

; 06.05 T 1169.77 DR 

0.269 DS 

-0.003 DT 

PL RINP 

173 *H 

205 *Y 

272 42 

173 -X 

331 -Y 

261 -2 

$$ PCON 
$$ PCON 
SS PCON 

62 RFX 
40 RPR 
100 RPT 

0 RLX 
30 RLR 
30 RLT 

0 RRX 
31 RRR 
31 RRT 

0 RRX 
61 RRR 
61 RRT 

62 RX 
107 RR 
169 RT 

671 OHL 
671 ONR 
1341 ONT 

GET • 1056.400 HHNN. 

SS 





SSION50 
SSI INLD 

3617.3 SNR 
187,5 HR 

0.00487 ECC 
149.8 HP 

29.31 INC 
28.46 INC 

279.22 RAN 
279.57 RAN 

25.90 ARC 

19.90 ARC 

153.85 TRR 
159.91 TRR 

SSI M50 
SSt SLV 

-21.13 PCH 
-114.37 PCH 

-57.48 YRM 
0.00 YRH 

98.13 ROL 
180.80 ROL 

•0.000 RX0 
0.000 RX8 

0.000 RYD 
-0.O65 RYB 

8.000 RZB 
0.000 PZ0 

PLI M50 
PLI PLV 

161.78 PCH 
-3.16 PCH 

-12.72 YRM 
-10.07 YRM 

97.80 ROL 
O.03 ROL 

0.000 RX0 
-0.011 RXB 

0.000 RYB 
0.069 RYB 

0.000 RZB 
0.014 RZB 

PLISN50 

PLiRSLV 

PLtRSBY 

-330744 X 
420806 X 
32821 X 

-281099 Y - 
-1849 Y 
1868 Y 

198034 2 -1033.87 DX 
226314 2 1452.11 DX 
475968 2 44.17 DX 

-759.41 DY 
-188.22 DY 
188.27 DY 

-401.51 D2 
189.08 D2 
1338.19 DZ 

TPl SENS 

477182 R 

3.26 S 

86.05 T 1330.80 DR 

0.323 DS 

8.004 DT 

PL RIMP 

173 *X 

205 +Y 

272 *2 

173 -X 

331 -Y 

261 -Z 

SS PCON 
SS PCON 
SS PCON 

62 RFX 
46 RFR 
188 RFT 

0 RLX 
30 RLR 
30 RLT 

0 RRX 
31 RRR 
31 RRT 

0 RAX 
61 PAR 
61 RAT 

62 RX 
107 RR 
169 RT 

671 ONL 
67! ONR 
1341 ONT 
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APPENDIX H 

GRAPHICAL OUTPUT DATA 







Figure HZ. Illustration of Graphical Plot for RSBY Plot Opti 



